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ABSTRACT 

The WMAP 5-year data provide stringent limits on deviations from the minimal, 6-parameter 
ACDM model. We report these limits and use them to constrain the physics of cosmic inflation 
via Gaussianity, adiabaticity, the power spectrum of primordial fluctuations, gravitational waves, 
and spatial curvature. We also constrain models of dark energy via its equation of state, parity- 
violating interaction, and neutrino properties such as mass and the number of species. We detect 
no convincing deviations from the minimal model. The 6 parameters and the corresponding 68% 
uncertainties, derived from the WMAP data combined with the distance measurements from the 
Type la supernovae (SN) and the Baryon Acoustic Oscillations (BAO) in the distribution of galaxies, 
are: VL b h 2 = 0-02267±g;ggg||, VL c h 2 = 0.1131 ± 0.0034, 0\ = 0.726 ± 0.015, n s = 0.960 ± 0.013, 
r = 0.084 ± 0.016, and Afj = (2.445 ± 0.096) x 10~ 9 at k = 0.002 Mpc" 1 . From these we derive 
cr 8 = 0.812 ± 0.026, H a = 70.5 ± 1.3 km s" 1 Mpc" 1 , n b = 0.0456 ± 0.0015, ft c = 0.228 ± 0.013, 
Q m h 2 = Q.1358±g:gg|^, z ro ion = 10.9 ± 1.4, and t = 13.72 ± 0.12 Gyr. With the WMAP data 
combined with BAO and SN, we find the limit on the tensor-to-scalar ratio of r < 0.22 (95% CL), 
and that n s > 1 is disfavored even when gravitational waves are included, which constrains the 
models of inflation that can produce significant gravitational waves, such as chaotic or power-law 
inflation models, or a blue spectrum, such as hybrid inflation models. We obtain tight, simultaneous 
limits on the (constant) equation of state of dark energy and the spatial curvature of the universe: 
-0.14 < 1 + w < 0.12 (95% CL) and -0.0179 < n k < 0.0081 (95% CL). We provide a set of ll WMAP 
distance priors," to test a variety of dark energy models with spatial curvature. We test a time- 
dependent w with a present value constrained as —0.33 < 1 + u>n < 0.21 (95% CL). Temperature 
and dark matter fluctuations are found to obey the adiabatic relation to within 8.9% and 2.1% 
for the axion-type and curvaton-type dark matter, respectively. The power spectra of TB and EB 
correlations constrain a parity- violating interaction, which rotates the polarization angle and converts 
E to B. The polarization angle could not be rotated more than —5.9° < Aa < 2.4° (95% CL) between 
the decoupling and the present epoch. We find the limit on the total mass of massive neutrinos of 
^m„ < 0.67 eV (95% CL), which is free from the uncertainty in the normalization of the large- 
scale structure data. The number of relativistic degrees of freedom, expressed in units of the effective 
number of neutrino species, is constrained as N e # = 4.4±1.5 (68%), consistent with the standard value 
of 3.04. Finally, quantitative limits on physically motivated primordial non-Gaussianity parameters 
are -9 < /^£ al < 111 (95% CL) and -151 < f^ 1 < 253 (95% CL) for the local and equilateral 
models, respectively. 

Subject headings: cosmic microwave background, cosmology: observations, early universe, dark matter, 
space vehicles, space vehicles: instruments, instrumentation: detectors, telescopes 
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the questions that were the focus of cosmology for the 
past 50 years: How old is the universe? How fast is it 
expanding? What is the size and shape of the universe? 
What is the composition of the universe? What seeded 
the formation of galaxies and large-scale structure? 

By accurately measuring the statistical properties of 
the microwave background fluctuations, WMAP has 
helped establish a standard cosmology: a flat ACDM 
model composed of atoms, dark matter and dark en- 
ergy, with nearly scale-invariant adiabatic Gaussian fluc- 
tuations. With our most recent measurements, WMAP 
has measured the basic parameters of this cosmology to 
high precision: with the WMAP 5-year data alone, we 
find the density of dark matter (21.4%), the density of 
atoms (4.4%), the expansion rate of the universe, the 
amplitude of density fluctuations, and their scale depen- 
dence, as well as the optical depth due to reionization 
(jDunklev et alJl2008f ). 

Cosmologists are now focused on a new set of ques- 
tions: What is the nature of the dark energy? What 
is the dark matter? Did inflation seed the primordial 
fluctuations? If so, what is the class of the inflationary 
model? How did the first stars form? Microwave back- 
ground observations from WMAP, Planck, and from the 
upcoming generation of CMB ground and balloon-borne 
experiments will play an important role in addressing 
these questions. 

This paper will discuss how the WMAP results, partic- 
ularly when combined with other astronomical observa- 
tions (mainly the distance measurements), are now pro- 
viding new insights into these questions through con- 
straints on gravitational waves and non-adiabatic (en- 
tropic) fluctuations, measurements of primordial non- 
Gaussianity, accurate determination of the primordial 
spectral index and the geometry of the universe, and lim- 
its on parity- violating interactions. 

This paper is one o f 7 papers on t h e ana lysis of the 
WMAP 5-year data: iHinshaw et~aTl (|2008l ) report on 
the da ta processi n g, ma p-making, and systematic error 
limits, iHill et ail (2008) on the physical optics model- 



ing of beams and the 5-year window functions (beam 
transfer functions), I Gold et abl (|2008h on the modeling, 
understanding, and subtraction of the tempe r ature and 
polarized foreground emission, IWright et alj (|2008l ) on 
the catalogue of p oint sources detec ted in the 5-year 
temperature data, N ofta et all (|2008f ) on the measure- 
ments o f the temperature an d polarization power spec- 
tra, and lDunklev et ail (|2008l ) on the parameter estima- 
tion methodology, the cosmological parameters inferred 
from the WMAP data alone, and comparison between 
different cosmological data sets. 

This paper is organized as follows. In § [2] we briefly 
summarize new aspects of our analysis of the WMAP 5- 
year temperature and polarization data. In § [3] we con- 
strain the spatial curvature of the observable universe, 



Gaussianity/adiabaticity/scale-invariance of the primor- 
dial fluctuations, and the amplitude of primordial grav- 
itational waves. We discuss their implications for the 
physics of the early, primordial universe. In § [4] we 
demonstrate that the power spectra of TB and EB cor- 
relations 16 which are usually ignored in the cosmologi- 
cal analysis, can be used to constrain a certain parity- 
violating interaction that couples to photons. In §[5] we 
explore the nature of dark energy, and in § |6] we study 
the properties of neutrinos in cosmology. We conclude in 

§m 

2. SUMMARY OF 5- YEAR ANALYSIS 
2.1. WMAP 5-year data: temperature and polarization 

With 5 years of observations of Jupit er and an exten - 
sive physical optics modeling of beams (jHill et al.ll2008h . 
our understanding of the beam transfer function, bi, 
has improved significantly: the fractional beam errors, 
Abi/bi, have been nearly halved in most Differencing As- 
semblies (DAs). In some cases, e.g., W4, the errors have 
been reduced by as much as a factor of 4. 

Many of the small-scale CMB experiments have been 
calibrated to the WMAP 3-year data at high multi- 
poles. Since the new beam model raises the 5-year 
power spectrum almost uniformly by ~ 2.5 % relative 
to th e 3-year power spectrum over I > 200 ( Hill et aTl 
2008), those small-scale CMB experiments have been 
under-calibrated by the same amount, i.e., ~ 2.5% in 
power, and 1.2% in temperature. For ex ample, the 
latest ACBAR data (jReichardt et al.1 l2008h report on 
the calibration error of 2.23% in temperature (4.5% in 
power), which is twice as large as the magnitude of mis- 
calibration; thus, we expect the effect of mis-calibration 
to be sub-dominant in the error budget. Note that the 
change in the beam is fully consistent with the 1-a er- 
ror of the previous WMAP beam reported in lPage et al.l 
(l2003h . Since the ACBAR calibration error includes the 
previous WMAP beam error, the change in the beam 
should have a minimal impact on the current ACBAR 
calibration. 

While we use only V and W bands for the cosmological 
analysis of the temperature data, the treatment of bi in 
Q band affects our determination of the point source con- 
tamination to the angular power spectrum, A ps 17 . The 
5-year estimate of the point source correction, A ps — 
0.011 ± 0.001 fiK 2 sv (jNolta et al.ll2008h . is slightly lower 
than the 3-year estim ate, A ps — 0.014 ± 0.003 [iK 2 sr 
(Hins haw et al.ll2007h . partly because more sources have 
been detected and masked by the 5-year source mask 
(390 sources have been detected in the 5-year tempera- 
ture data, whereas 323 sourc es were detected in the 3- 
year data (| Wright et alJl2008ft ). 

Note that the uncertainty in A ps has been reduced 
by a factor of 3. The uncertainty in the previous esti- 
mate was inflated t o inc lude the lower value found by 
iHuffenberger et all (|2006f ) (0.011) and higher value from 

16 Here, "TB" refers to the power spectrum of a cross-correlation 
between the temperature and B-mode polarization, while "EB" 
refers to a correlation between the E-mode and B-mode polariza- 
tion. 

17 This quantity, A pa , is the value of the power spectrum, C;, 
from unresolved point sources in Q band, in units of the antenna 
temperature. To con vert this value to the thermodynamic units, 
use C ps = 1.089 A pa jNolta et al.ll2008h . 



WMAP 5-year Cosmological Interpretation 



3 



TABLE 1 

Summary of the cosmological parameters of ACDM model and the corresponding 68% intervals 



Class Parameter WMAP 5-year ML" WMAP+BAO+SN ML WMAP 5-year Mean 6 WMAP+BAO+SN Mean 



Primary 


100tt b h 2 


2.268 


2.262 


2.273 ± 0.062 


9 9R7+0.058 
— 0.059 




n c h 2 


0.1081 


0.1138 


0.1099 ±0.0062 


0.1131 ± 0.0034 




n A 


0.751 


0.723 


0.742 ± 0.030 


0.726 ± 0.015 




n s 


0.961 


0.962 


Q63+ 014 


0.960 ± 0.013 






0.089 


0.088 


0.087 ±0.017 


0.084 ± 0.016 




A|(fco c ) 


2.41 x 10~ 9 


2.46 x 10" 9 


(2.41 ±0.11) x 10" 9 


(2.445 ± 0.096) x 10" 9 


Derived 




0.787 


0.817 


0.796 ±0.036 


0.812 ± 0.026 




Ho 


72.4 km/s/Mpc 


70.2 km/s/Mpc 


71.9±f% km/s/Mpc 


70.5 ± 1.3 km/s/Mpc 




n b 


0.0432 


0.0459 


0.0441 ± 0.0030 


0.0456 ± 0.0015 




n c 


0.206 


0.231 


0.214 ±0.027 


0.228 ± 0.013 






0.1308 


0.1364 


0.1326 ±0.0063 


n iqco+0.0037 

u.iooo_ 0036 




z ■ f 

^•rcion 


11.2 


11.3 


11.0 ± 1.4 


10.9 ± 1.4 




to 9 


13.69 Gyr 


13.72 Gyr 


13.69 ±0.13 Gyr 


13.72 ± 0.12 Gyr 
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(2008). "ML" refers to the Maximum Likelihood parameters 

(2008). "Mean" refers to the mean of the posterior distribution of each parameter 
(200g). "ML" refers to the Maximum Likelihood parameters 
Dunklev et a l. (2008). "Mean" refers to the mean of the posterior distribution of each parameter 
e k = 0.002 Mpc" 1 . A^(fc) = k' i PT Z (k)/(2n 2 ) (Eq. [15]) 

■'"Redshift of reionization," if the universe was reionized instantaneously from the neutral state to the fully ionized state at 

^reion 

9 The present-day age of the universe 



our original estimate (0.017). Much of the discrepancy 
between these estimates is due to the multipole range 
over which A ps is fit. With the improved beam model 
from the 5-year analysis, the dependence on the multi- 
pole range has disappeared, and thus we no l o nger need 
to inflate our uncertainty. See iNolta et alj (|2008| ) for 
more details. 

The method for cleaning foreground emission in both 
temperature and polarization data is the same as we 
used for the 3-year data, i.e., the template-based clean- 
ing method described in § 5.3 o f Hinsh aw et al.l (|2007f ) 
for tem perature and § 4.3 of lPage et aLl (|2007D for polar- 
ization. I Gold et all (|2008| ) describe the results from the 
template cleaning of the 5-year data with the new coeffi- 
cients . In addition. iGold et aLl (|2008f ) and lDunklev et all 
explore alternative modelings of the foreground 



polarization power spectra. 

We have improved our treatment of the temperature 
and polarization power spectra at lower multipoles, as 
described below. 

Low-l temperature - We use the Gibbs sampling tech- 
nique and the Blackwell-Rao (BR) estimator to evalu- 
ate the likelihood of the t e mperature power spectrum at 
I < 3 2 dJewell et al.ll200l I Wandeltl 120031: jWandelt elaTl 



emission. All of these methods give consistent results. 
I Gold et aLl (|2008t ) also describe definitions of the new 
masks, KQ75 and KQ85, that replace the previous 
masks, KpO and Kp2, that are recommended for the anal- 
ysis of Gaussianity tests and the power spectrum, respec- 
tively. 

The method for measuring the TT and TE spectra 
at higher multipoles, i.e., I > 33 for TT and / > 24 
for TE, is also the sa me as we used for the 3-year data 
(|Hinshaw et al.ll2007h . As for the estimation of the cos- 
mological parameters from these spectra, we now include 
the weak gravitational lensing effect of CMB due to the 
interv ening matter fluctuations (see iLewis fc Challinorl 
I2006L for a review), which was not included in the 3- 
year analysis. We continue to marginalize over a po- 
tential contribution from the Sunyaev-Zel'dovich effect 
(SZE), using exactly the same template SZE power spec- 
trum that we used for t he 3-year analysis: Cf ZE from 
iKomatsu fc Seliakj (|2002t ) with fl m = 0.26, fl b = 0.044, 
h = 0.72, n, = 0. 97, and a s = 0.80 (see also § 2.1 of 
ISpergel et~al1l2007t >. We continue to use the V and W 
band data for estimating the high-/ temperature power 
spectrum, and the Q and V band data for the high-/ 



200i lO'Dwver et a.l.ll2004J: lEriksen et aLl l2004 12007c, 
Chu et al.ll2005tlLarson et al.ll2007l ). For the 3-year anal- 
ysis we used the resolution 4 Internal Linear Combination 
(ILC) temperature map (AWc = 16) with a Gaussian 
smoothing of 9.183° (FWHM). Since the ILC map has 
an intrinsic Gaussian smoothing of 1°, we have added 
an extra smoothing of 9.1285°. We then evaluated the 
likelihood directly in the pixel space for a given C;. For 
the 5-year analysis we use a higher resolution map, the 
resolution 5 ILC map (N s ia c = 32) with a smaller Gaus- 
sian smoothing of 5° (FWHM). The potential foreground 
leakage due to smoothing is therefore reduced. The BR 
estimator has an advantage of being much faster to com- 
pute, which is why we have adopted the Gibbs sampling 
and the BR estimator for the 5-year data release. We 
have confirmed that both the resolution 4 pixel-based 
likelihood and the resoluti on 5 Gibbs-based lik elihood 
yield consistent results (see lDunklev et all 120081 . for de- 
tails). Both options are made publicly available in the 
released likelihood code. 

Low-l polarization - While we continue to use the direct 
evaluation of the likelihood of polarization power spec- 
tra in pixel space from coadded resolution 3 (-/V s id = 8) 
polarization maps (Stokes Q and U maps), we now add 
the Ka band data to the coadded maps; we used only 
Q and V band data for the 3-year analysis. We believe 
that we understand the polarized foreground emission 
(dominated by synchrotron, traced well by the K-band 
data) in the Ka band da ta well enough to justify the in- 
clusion of the Ka band (|Gold et all [2008). This, with 2 
years more integration, has led to a significant reduction 
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TABLE 2 

Summary of the 95% confidence limits on deviations from the simple (flat, Gaussian, adiabatic, 

power-law) ACDM model 



Section 



Name 



Type 



WMAP 5-year 



WMAP+BAO+SN 



3X31 



3,1 



m 

i a 



Gravitational Wave" 
Running Index 
Curvature' 1 
Curvature Radius 9 

Gaussianity 

Adiabaticity 

Parity Violation 
Dark Energy 

Neutrino Mass s 
Neutrino Species 



No Running Ind. 
No Grav. Wave 

Positive Curv. 
Negative Curv. 
Local 

Equilateral 
Axion 
Curvaton 
Chern-Simons" 
Constant w° 
Evolving w{z) q 



r < 0A3 b 
-0.090 < dn s /d\nk < 0.019 c 
-0.063 <U k < 0.017 e 
R CUIV > 12 h^Gpc 
i?,cmv > 23 h~ 1 Gpc 
-9 < /^£ al < 
-151 < /^£ U < 253 s 
a < 0.16 J 
a-i < 0.011' 
-5.9° < Aa< 2.4° 
-1.37 < 1 + w < 0.32P 
N/A 
Em„< 1.3 eV 
iVeff > 2.3" 



r < 0.22 
-0.068 < dn s /d\nk < 0.012 
-0.0179 <fl k < 0.0081/ 
iJcurv > 22 /i~ 1 Gpc 
-Rcmv > 33 /i _1 Gpc 
N/A 
N/A 
ao < 0.072* 
ot-i < 0.0041™ 
N/A 

-0.14 < 1 + w < 0.12 
-0.33 < 1 + wq < 0.21 r 

E%< 0.67 eV" 
Af cff =4.4=fc 1.5 m (68%) 



"In the form of the tensor-to-scalar ratio, r, at k = 0.002 Mpc 
punklev et~aT1 |2008]1 
TOunklev etall {2008} 

d (Constant) dark energy equation of state allowed to vary (w ^ — 1) 

e With the HST prior, H = 72 ± 8 km/s/Mpc. For w = -1, -0.052 < U k < 0.013 (95% CL) 
/ For w = -1, -0.0178 < Q fc < 0.0066 (95% CL) 

9 Rcurv = (c/tf )/vW = 3/vW /l-^pc 

^Cleaned V+W map with Z ma x = 500 and the KQ75 mask, after the point source correction 
'Cleaned V+W map with l max = 700 and the KQ75 mask, after the point source correction 
TOunklev etaD lf2008) 



k ln terms of the adiabaticity deviation parameter, 
are found to obey the adiabatic relation (Eq. 1361 ) 



«5 (c : 7) = 

a at 

to 8.9% 



Y^a/3 (Eq. 1391 ). the axion-like dark matter and photons 



IDunklev et"aH (;2008) 
"In terms of the adiabaticity deviation parameter, 



\/ct/3 (Eq. 1391 ). the curvaton-like dark matter and 



photons are found to obey the adiabatic relation (Eq. 1361 ) to 2.1%. 
"For an interaction of the form given by [c/>(t)/M]F a pF a P , the polarization rotation angle is Aa : 
°For spatially curved universes (f2 fc ^ 0) 
PWith the HST prior, H = 72 ± 8 km/s/Mpc 
? For a flat universe (Q k = 0) 
r wo = w(z = 0) 
'J]m, = 94(n^fe 2 ) eV 



M~ 1 ff<]> 



1 Dun klev etaLI (2008) 
"For w = —1. For w ^ 
IDunklev etaLI (|2008) 
"With the HST prior, H = 72 ± 8 km/s/Mpc. The 95% limit is 1. 



-1, £m„ < 0.80 eV (95% CL) 



< N efi < 7.6 



of the noise power spectra (averaged over I = 2 — 7) in 
the polarization EE and BB power spectra by a factor 
of as much as 2.3 compared to the 3-year analysis. As a 
result, the EE power spectrum averaged over I — 2 — 7 
exceeds the noise by a factor of 10, i.e., our measurement 
of the EE power spectrum averaged over I = 2 — 7 is now 
limited by cosmic variance and the possibility of residual 
foreground emission and/or systematic errors 18 , rather 
than by noise. In addition, we have added a capability 
of computing the likelihood of TB and EB power spectra 
to the released likelihood code. This allows us to test 
models in which non-zero TB and EB correlations can 
be generated. We discuss this further in § 01 

We continue to use the Markov Chain Monte Carlo 
(MCMC) technique to explore the posterior distribu- 
tion of cosmological parameters given the measured tem- 
perature and polarization power spectra. For details 
on the implementat ion and convergence criteria, see 
IDunklev et al] (|2008T ). 

18 For o ur limits on the residual polarized foreground contami- 
nation, see IDunklev etU] |[2008T) . 



2.2. Comments on systematic errors in the cosmological 
parameters derived from the WMAP 5-year data 

iHinshaw et al] (|2008l ) give extensive descriptions of our 
limits on the systematic errors in the WMAP 5-year 
temperature and polarization data. With the improved 
treatment of beams and gain calibrations, we are confi- 
dent that the instrumental systematic errors in the cos- 
mological results derived from the temperature data are 
negligible compared to the statistical errors. As for the 
polarization data, we find that the W band polarization 
data still contain the systematic errors that we do not 
understand fully, and thus we do not use the W band po- 
larization for the cosmological analysis. We do not find 
any evidence for the unaccounted instrumental system- 
atic errors in the other bands that we use for cosmology. 

The most dominant systematic errors in the cosmolog- 
ical analysis are the foreground emission. Since CMB 
dominates over the foreground emission in the temper- 
ature data in V and W bands outside the galaxy mask, 
and we also reduce the sub-dominant foreground contam- 
ination at high Galactic latitudes further by using the K- 
and Ka-band data (for synchrotron emission), the exter- 
nal Ha map (for free- free emission), and the external dust 
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TABLE 3 

Sound horizon scales determined by the 

wmap 5-year data. cmb: the sound 
horizon scale at the photon decoupling 
epoch, z t , imprinted on the cmb power 
spectrum; Matter: the sound horizon 
scale at the baryon drag epoch, z d , 
imprinted on the matter (galaxy) power 
spectrum. 





Quantity 


Eq. 


5-year WMAP 


CMB 




deal 


1090.51 ± 0.95 


CMB 


r a (z*) 


© 


146.8 ± 1.8 Mpc 


Matter 






1020.5 ± 1.6 


Matter 


r s (z d ) 




153.3 ± 2.0 Mpc 



map, the systematic errors from the foreground emission 
are unimportant for the temperature data, even at the 
lowes t multipoles whe re the foreground is most impor- 
tant jGoiditinilool). 

We, however, find that the uncertainty in our mod- 
eling of the polarized foreground is not negligible com- 
pared to the statistical errors. For the 5-year polariza- 
tion analysis we have used two independent foreground- 
cleaning algorithms: one based upon the t emplate fit- 
ting (as developed for the 3-year analysis; see lPage et al.l 
l2007f) . and the other based upon the Gibbs sampling 
(|Dunklev et all 120081 ). The optical depth, r, is the 
parameter that is most affected by the uncertainty in 
the polarized foreground model. The template fitting 
method gives r = 0.087 ± 0.017. The Gibbs sampling 
method gives a range of values from r = 0.085 ± 0.025 to 
r = 0.103 ± 0.018, depending upon various assumptions 
made about the properties of the polarized synchrotron 
and dust emission. Therefore, the systematic error in r 
is comparable to the statistical error. 

This has an implication for the determination of the 
primordial tilt, n Sl as there is a weak correlation between 
n s and r (see Fig. Q]): for r = 0.087 we find n s = 0.963, 
while for r = 0.105 we find n s = 0.98. Since the statisti- 
cal error of n s is 0.015, the systematic error in n s (from 
the polarized foreground) is comparable to the statisti- 
cal one. The other parameters that are correlated with 
n s , i.e., the baryon density (Fig. QJ, the tensor-to-scalar 
ratio (Fig. [3]), and the amplitude of non-adiabatic fluc- 
tuations (Fig. [5] and O , would be similarly affected. For 
the parameters that are not correlated with n s or r the 
systematic errors are insignificant. 

2.3. External data sets: Hubble constant, luminosity 
and angular diameter distances 

Aside from the CMB data (including the small-scale 
CMB measurements), the main external astrophysical 
results that we shall use in this paper for the joint cos- 
mological analysis are the following distance-scale indi- 
cators: 

• A Gaussian prior on the present-day Hubble's con- 
stant from the Hubble Ke y Project final results , 
H = 72±8km s" 1 Mpc" 1 (jFreedman et al.ll2001f ). 
While the uncertainty is larger than the WMAP's, 
determination of Ho for the minimal ACDM model 
(see Table [1} , this information improves upon lim- 
its on the other models, such as models with non- 
zero spatial curvature. 



• The luminosity distances out to Type la super- 
novae with their absolute magnitudes marginalized 
over uniform priors. We use the "union" super- 
nova samples compiled by Kowals ki et al.l ([2008) . 
The union compilation contains 57 nearby (0.015 < 
z < 0.15) Type la supernovae, and 250 high-;? 
Type la supernovae, after selection cuts. The high- 
z samples contain the Type l a supernovae from 
the Hubble Space Tele scope (|Knop et all 120031 : 
iRiess et al.l 12004 12007ft. the Super Nova Legacy 
Survey (SNLS) (jAstier et all I2006D . the Equa- 
tion of State: SupErNo vae trace Cosmic Expan - 
sion (ESSENCE) survey (jWood-Vasev et al.l[200l . 
as well as those used in the original papers of 
the discovery of t h e acceleration of the u niverse 
(|Riess et al.l 119981: iPerlmutter et all I1999D. and 
the sa mples from Barris et al.l (|2004l) ; iTonrv et al.l 
(120031) The nearby T ype l a supernovae are taken 
from iHamuvetaLl d 19961): IRiess et all (Il999f) : 
Jha et al.l(l2006l ) :lKrisciunas et al.l (|200lLl2004al|b1 ). 
Kowalski et alj (|2008f ) have processed all of these 
sample s using the sam e light curve fitter called 
SALT (|Guv et all 120051 ). which allowed them to 
combine all the data in a self-consistent fashion. 
The union compilation is the largest t o date . The 
previous compilation by iDavis et alJ (|2007h used 
a smaller number of Type la supernovae, and did 
not use the same light curve fitter for all the sam- 
ples. We examine the difference in the derived 
ACDM cosmological parameters between the union 
compilation and Davis et al.'s compilation in Ap- 
pendix [DJ While we ignore the systematic errors 
when we fit the Type la supernova data, we exam- 
ine the effects of systematic errors on the ACDM 
parameters and the dark energy parameters in Ap- 
pendix |Dj 

• Gaussian priors on the distance ratios, r s j Dy(z), 
at z = 0.2 and 0.35 measured from the Baryon 
Acoustic O scillations (BAO) in the distribution 
of galaxies (jPercival et al.ll2007l ). The CMB ob- 
servations measure the acoustic oscillations in the 
photon-baryon plasma, which can be used to mea- 
sure the angular diameter distance to the photon 
decoupling epoch. The same oscillations are im- 
printed on the distribution of matter, traced by the 
distribution of galaxies, which can be used to mea- 
sure the angular distances to the galaxies that one 
observes from galaxy surveys. While both CMB 
and galaxies measure the same oscillations, in this 
paper we shall use the term, BAO, to refer only 
to the oscillations in the distribution of matter, for 
definiteness. 

Here, we describe how we use the BAO data in more 
detail. The BAO can be used to measure not only the 
angular diameter distance, Da{z), through the clustering 
perpendicular to the line of sight, but also the expansion 
rate of the universe, H(z), through the clustering along 
the line of sight. This is a powerful probe of dark en- 
ergy; however, the accuracy of the current data do not 
allow us to extract Da(z) and H(z) separately, as one 
can barely measu re BAO in the spherica lly averaged cor- 
relation function (|Okumura 



() 
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The spherical a verage gives us the fol lowing effective 
distance measure (|Eisenstein et al.ll2005l ): 



Dv(z) 



{l + zYD\{z) 



H{z) 



1/3 



(1) 



where Da(z) is the proper (not comoving) angular diam- 
eter distance: 



D A (z) = 



fk 



dz 



(z') 



o (i + «)vTO 



(2) 



where fk[x] = sinx, x, and sinhx for < (k = 1), 
Slfc = (k = 0), and flk > (k = —1), respectively. 

There is an additional subtlety. The peak positions 
of the (spherically averaged) BAO depend actually on 
the ratio of Dy(z) to the sound horizon size at the so- 
called drag epoch, Zd, at which baryons were released 
from photons. Note that there is no reason why the 
decoupling epoch of photons, z», needs to be the same as 
the drag epoch, za- They would be equal only when the 
energy density of photons and that of baryons were equal 
at the decoupling epoch - more precisely, they would be 
equal only when R(z) = 3p b /(Ap^) = (3S1 & /4J1 T )/(1 + 
z) ~ 0.67(f2 b /i 2 /0.023) [1090/(1 + z)] was unity at z = 
z». Since we happen to live in the universe in which 
fibft 2 ~ 0.023, this ratio is less than unity, and thus the 
drag epoch is slightly later than the photon decoupling 
epoch, Zd < z*. As a result, the sound horizon size at 
the drag epoch happens to be slightly larger than that 
at the photon decoupling epoch. In Table [3] we give the 
CMB decoupling epoch, BAO drag epoch, as well as the 
corresponding sound horizon radii that are determined 
from the WMAP 5-year data. 

We use a fitting formula for Zd proposed by 
lEisenstein k. Hul (|1998l ): 



1291(f2 m /i 2 ) 



2\0.251 



1 + 0.659(f2 m /i 2 ) 



.828 



[l + 6 1 (r2 h / l 2 ) 62 ] , (3) 



where 



&i=0.313(f2 m /i 2 r°' 419 [1 
6 2 = 0.238(r2 m /i 2 ) ' 223 . 



0.607(f2 m /i 2 ) 



2\0.6741 



(4) 
(5) 



In this paper we frequently combine the WMAP data 
with r s (zd)/Dv(z) extracted from the Sloan Digital 
Sky Survey (SDSS) and th e Two Degree Field Galaxy 
Redshift Survey (2dFGRS) dPercival et al.ll2007ri . where 
r s (z) is the comoving sound horizon size given by 



,(*) 



c 

V3 



V(i+2) 



da 



a 2 #(a) v / l + (31V4f2 7 )a' 



(6) 



where Oy = 2.469 x lO" 5 ^" 2 for T cmh = 2.725 K, and 



H{a) = H 



a 3 a 4 a 2 a 3(i+™ctf(a)) 



1/2 



(7) 



The radiation density parameter, Q r , is the sum of pho- 
tons and relativistic neutrinos, 



a r = Oy (1 + 0.2271A cff ) 



(8) 



where N c f[ is the effective number of neutrino species 
(the standard value is 3.04). For more details on N e g, 
see 5 16.21 When neutrinos are non-relativistic at a, one 



needs to reduce the value of N e g accordingly. Also, the 
matter density must contain the neutrino contribution 
when they are non-relativistic, 

On = fic + V b + tt v , (9) 

where £l v is related to the sum of neutrino masses as 



(10) 



94h 2 eV ' 

For more details on the neutrino mass, see § 16.11 

All the density parameters refer to the values at the 
present epoch, and add up to unity: 



- - r ( 



f2. r + Ok + n A = i. 



(11) 



Throughout this paper, we shall use J1a to denote the 
dark energy density parameter at present: 



= &de{z = 0). 



(12) 



Here, w e s(a) is the effective equation of state of dark 
energy given by 

An a 

w ff(a) = i / d\na' w(a'), (13) 

In a J 

and w(a) is the usual dark energy equation of state, i.e., 
the dark energy pressure divided by the dark energy den- 
sity: 

P de {a) 



w(a) = 



Pde{a) 



(14) 



For vacuum energy (cosmological constant), w does not 
d epend on time, and w = — 1. 

iPercival et al.1 (|2007l ) have determined r s (zd)/Dv(z) 
out to two redshifts, z = 0.2 and 0.35, as r s {zd)/Dv{z — 
0.2) = 0.1980 ± 0.0058 and r s (z d )/D v (z = 0.35) = 
0.1094 ± 0.0033, with a correlation coefficient of 0.39. 
We follow th e desc ription given in Appendix A of 
IPercival et alj (|2007t ) to implement these constraints in 
the likelihood code. We have checked that our calcu- 
lation of r s (zd) using the formulae above (including Zd) 
matches the value that they quote 19 , 111.426 h" 1 Mpc, 
to within 0.2 Mpc, for their quoted cosmological pa- 
rameters, fi m = 0.25, n b h 2 = 0.0223, and h = 0.72. 

We have decided to use these results, as they measure 
the distances only, and are not sensitive to the growth of 
structure. This properly enables us to identify the infor- 
mation added by the external astrophysical results more 
clearly. In additio n to these, we shall a lso use the BAO 
measurement by lEisenstein et all (I2005D 20 . and the flux 
power spectrum of Lya forest from ISeliak etTai] (p006T ) 
in the appropriate context. 

For the 3-year data analysis in lSpergel et alj (|2007f ) we 
also used the shape of the galaxy power spectra measured 
from the SDS S main sample as well as the Luminous 
Red Galaxies (iTeemark et al.ll2004E 12001 . and 2dFGRS 
(|Cole et al.ll2005D . We have found some tension between 
these data sets, which could be indicative of the degree 

19 In Pcrcival ct al. (2007) the authors use a different notation 
for the drag redshift, z„, ins tead of z^. We have co nfirmed that 
they have used equation (6) of Eiscnstein & Hu ( 1998) for r s , which 
makes an explicit use of the drag redshift (Will Percival, private 
communication) . 

20 We use a Gaussian prior on A = Dy(z = 

0.35)^/n m /^/(0.35c) = 0.469(n s /0.98)-°- 35 ± 0.017. 
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by which our understanding of non-linearities, such as 
the non-linear matter clustering, non-linear bias, and 
non-linear redshift space distortion, is st i ll limi ted at low 
redshifts, i.e., z < 1. See lDunklev et alj (120081) for more 
detail ed study on this issue. See also iSanchez fe Cold 
(2008) on the related subject. The galaxy power spec- 
tra should provide us with important information on the 
growth of structure (which helps constrain the dark en- 
ergy and neutrino masses) as our understanding of non- 
linearities improves in the future. In this paper we do 
not combine these data sets because of the limited un- 
derstanding of the consequences of non-linearities. 

3. FLAT, GAUSSIAN, ADIABATIC, POWER-LAW ACDM 
MODEL AND ITS ALTERNATIVES 



the idea 
period of 



that 



the uni- 
acceler- 



The theory of inflation, 
verse underwent a brief 
ated expansion jStarobin skv _ 

Kazanasl Il98nt iGuthl Il981h iSatol TTqSH ILindd llQSl i 



period, ol rapid 
^vl [19791 : jStarobinskvl [1981 



Albrecht fc Steinhardtl fl982). has become an indispens 



able building block of the standard model of our universe. 

Models of the early universe must explain the follow- 
ing observations: the universe is nearly flat and the 
fluctuatio ns observed by WMA P appear to be nearly 
Gaussian (Komats u et alJ l2003t). scale-invariant, super- 
horizon, and adiabatic (ISpergel fc Zaldarriagal 119971 : 
ISpergel et all l2003t iPeiris et alj l2003h . Inflation mod- 
els have been able to explain these properties suc- 
cessfully dMukhanov fc Chibisovl pL98lt [ Hawking j!982t 
IStarobinskylll982tlGuth fc Pilll982tlBardeen et alJll983f ). 

Although many models have been ruled out observa- 
tionally (see iKinney et "all 120061: lAlabidi fc Lythl l2006at 
iMartin fc R .ingcval 200i| for recent surveys), there are 
more tha n one hundred candidate inflation model s avail- 
able (see ILTddle fc Lvthll200Clt iBassett et alll2006t ILindd 
120081 for reviews). Therefore, we now focus on the ques- 
tion, "which model is the correct inflation model?^ This 
must be answered by the observational data. 

On the other hand, an inflationary expansion may not 
be the only way to solve cosmological puzzles and cre- 
ate primordial fluctuations. Contraction of the primor- 
dial universe followed by a bounce to expansion can, 
in principle, make a set of predictions th at are qualita- 
tively similar to those of inflat i on mo d els dKhourv et al. | 
it 



200 1 l2002allbl: iKhourv et al l 12001 iBuchbinder et al 
20071 [ 2008: Ko vama fc Wanddl2007tlKovama et alj|2007t 
Creminelli fc Senatord l2007) , although building concrete 



models and making robust p redictions have been chal- 
lengi ng (|Kallosh et al.ll200l|[?l : lLind3[200l iKallosh et all 
2008). 

There is also a fascinating possibility that one can learn 
something about the fundamental physics from cosmo- 
logical observations. For example, a recent progress in 
implementing de Sitter v acua and inflation in the con- 
text of String Theory (see lMcAllister fc Silversteinll2008l 
for a review), makes it possible to connect the cosmolog- 
ical observations to the ingredients of the fundamental 
physics via their predicted properties of inflation such as 
the shape of the power spectrum, spatial curvature of the 
universe, and non-Gaussianity of primordial fluctuations. 

3.1. Power spectrum of primordial fluctuations 
3.1.1. Motivation and Analysis 



The shape of the power spectrum of primordial cur- 
vature perturbations, Pn(k), is one of the most power- 
ful and practical tool for distinguishing among inflation 
models. Inflation models with featureless scalar-field po- 
tentials usually predict tha t Pji(k) is nearly a power-law 
(jKosowskv fc Turnerlll995l ) 



n s (k.Q) — l-\-^dn s / d In k 



(15) 

which gives an ap- 



Here, A^(fe) is a useful quantity, 
proximate contribution of 1Z at a given scale per log- 
arithmic interval in k to the total variance of 7£, as 
(7^ 2 (x)) = J dink A^(fc). It is clear that the special 
case with n s = 1 and dn s /d\nk = results in the "scale- 
invariant spectrum," in which the contributions of 1Z at 
any scales per logarithmic interval in k to the total vari- 
ance are equal. (Hence the term, "scale invariance.") 
Following the usual terminology, we shall call n s and 
dn s /dlnk the tilt of the spectrum and the running in- 
dex, respectively. We shall take fco to be 0.002 Mpc~ . 

The significance of n s and dn s /dhik is that differ- 
ent inflation models motivated by different physics make 
specific, testable predictions for the values of n s and 
dn s /d\nk. For a given shape of the scalar field poten- 
tial, V (</)), of a single-field model, for instance, one finds 
that n s is given by a combination of the first deriva- 
tive and second derivative of the potential, 1 — n s — 
ZM^V'/Vf - 2M* l (V"/V) (where M p 2 ; = 1/(8ttG) 
is the reduced Planck mass), and dn s /d\nk is given 
by a combination of V'/V, V"/V, and V" /V (see 
iLiddle fc Lvtl3l200a for a review). 

This means that one can reconstruct the shape of V(4>) 
up to the first three derivatives in this way. As the expan- 
sion rate squared is proportional to V(<p) via the Fried- 
mann equation, H 2 = V/(3M? l ) ) one can reconstruct 
the expansion history during inflation from measuring 
the shape of the primordial power spectrum. 

How generic are n s and dn s /d\nk? They are moti- 
vated physically by the fact that most inflation models 
satisfy the slow- roll conditions, and thus deviations 
from a pure power-law, scale-invariant spectrum, n s = 1 
and dn s /dlnk = 0, are expected to be small, and the 
higher-order derivative terms such as V"" and higher 
are negligible. On the other hand, there is always 
a danger of missing some important effects, such as 
sharp features, when one relies too much on a simple 
parametrization like this. Therefore, a number of people 
have investigated various, more general ways of recon- 
struc t ing the shape of Pyi(k) d Matsumiya et al. [2002L 
20031: iMukheriee fc Wanel l2003bOct iBridle et all 120031: 



2004 



Kogo et alJ 
2004 lHannestad 



Wang iJULMb.a.c; LSrid le et al.l \Z UU6: 
Kogo et all 120051: iHu fc Okamotol 



120041: iSh aficloo fc Souradeepf 200 



Sealfon et al.1 l200a Tocchini-Valcnti ni et al.l 200E 
2006t ISpergel et all 1200 71: IShafielo o fc Souradeepl 2007 
Verde fc Peirid I2008D . and V(6) (iLidsev et alf Il997 



Grivell fc Liddld 120001: [Kadota et al 



2006!; Ees gourgucs fc Va lkenburg 20071: 
2007). 



20051: iCovi et"aLl 



Powell fc Kinnevl 



These studies have indicated that the parametrized 
form (Eq. [15]) is basically a good fit, and no significant 
features were detected. Therefore, we do not repeat this 
type of analysis in this paper, but focus on constraining 
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TABLE 4 

Primordial tilt n a , running index dn a /dlnk, and tensor-to-scalar ratio r 



Sec. 


Model 


Parameter" 


5-year WMAP 6 


5-year WMAP+CMB C 


5-year WMAP+ACBAR08 d 


5-year WMAP+BAO+SN 


S 13.1.21 


Power-law 


n a 


u.ytM_ 015 


0.960 ±0.014 


0.964 ± 0.014 


0.960 ±0.013 


§ |3.1.3| 


Running 


n a 

dn a /d\nk 


, „„-, +U.U54j 
i.uoi_ 055 

-0.037 ±0.028 


1 059 +uubi 
-0.053 ± 0.025 


1.031 ± 0.049 

-0 n35+ ' 024 
u.u^o_ () 025 


1 -" LI -0.043 

-0.028 ±0.020* 


§|3.2.4| 


Tensor 


n a 
r 


0.986 ± 0.022 
< 0.43 (95% CL) 


0.979 ±0.020 
< 0.36 (95% CL) 


u.y»t>_ 020 
< 0.40 (95% CL) 


0.970 ±0.015 
< 0.22 (95% CL) 


§ |3.2.4| 


Running 
+Tensor 


n a 
r 

dn a /d\nk 


1.08f_ 073 
< 0.58 (95% CL) 
-0.050 ±0.034 


i 1Q7 +U.UVb 
1 ' 1 ' a < -0.071 

< 0.64 (95% CL) 
-0 072+ - 031 

u ' u ' -0.030 


i.uoo_ 062 
< 0.54 (95% CL) 
-0.048 ± 0.027 


l.U89_ 06g 
< 0.55 (95% CL) ; 
-0.053 ±0.028 m 



"Defined at kg = 0.00 2 Mpc" 
tounklev etaLl ll2008l1 



"CMB" includes t he sma ll-scale CMB measur ements fro m CB I (Mason ct al. 2003; Sievcrs ct al. 2003, 2007; Pearson ct al. 2003; 
I Readhead eTaIll2004T), VSA lIDickinson et al.ll200l) . ACBAR IIKuo et al.H2004L 120071) . and BOOMERanG l IRuhl eraLll2003l ; IMontrov et al.l 

12006k IPiacentini et alj|200D 

d "ACBAR08" is the complete ACBAR data set presented in Rcichardt ct al. (2008). We used the ACBAR data in the multipole range of 
900 <l < 2000 
e Defined at fcp = 0.00 2 Mpc" 1 
IDunklev etaLl IpOOl l 

9 "CMB" includes t he sma ll-scale CMB measur ements fro m CBI (IMason et al.l |2003|; ISievers et all |2003|. [20071; IPearson et al.l |2Q03|; 
Rcadhead ct al. 200% VSA. lIDickinson et al .112001) , ACBAR (IKuo et al.H2004L 120071) . and BOOMERanG l IRuhl eraL| [2003; Mon trov etaLl 



2006k IPiacentini et al.ll20Tj6D 



' ll ACBAR08" is the complete ACBAR data set presented in |Rcichardt ct al. (2008). We used the ACBAR data in the multipole range of 
900 < I < 2000 

'At the pivot point for WMAP only, k p i vot = 0.080 Mpc -1 , where n a and dn B /d\nk are uncorrelated, n s (k p i vo t) 
^At the pivot point for WMAP+BAO±SN, k p i vo t = 0.106 Mpc" 1 , where n B and dn s /d Ink are uncorrelated, n s (k p 
*With the Lyct forest data (Scliak ct al. 2006), dn s /d\nk = -0.012 ± 0.012 
'With the Lya forest data (Scliak ct al. 2006), r < 0.28 (95% CL) 
m With the Lya forest data l|Seliak et al.l l2006L dn s /d\nk = -0.0171q 



■■ 0.963 ± 0.014 
kpivot) = 0.961 ± 0.014 




0.92 0.94 0.96 0.98 1.00 



0.92 0.94 0.96 0.98 



1.00 



0.92 0.94 0.96 0.98 

n 



1.00 



Fig. 1. — Constraint on the primordial tilt, n s (§ 13. 1.21 1. No running index or gravitational waves are included in the analysis. (Left) 
One-dimensional marginalized constraint on n a from the WMAP-only analysis. (Middle) Two-dimensional joint marginalized constraint 
(68% and 95% CL), showing a strong correlation between n a and Qi,h 2 . (Right) A mild correlation with t. None of these correlations are 
reduced significantly by including BAO or SN data, as these data sets are not sensitive to Qi,h 2 or r; however, the situation changes when 
the gravitational wave contribution is included (see Fig. [3}. 

WMAP 5-year temperature and polarization data, the 
small-scale CMB data, and/or BAO and SN distance 
measurements. In Table [4] we summarize our results pre- 
sented in § 13X21 13X31 and I3T41 



the parametrized form given by Eq. (|15p . 

Finally, let us comment on the choice of priors. We im- 
pose uniform priors on n s and dn s /d\nk, but there are 
other possibilities for the choice of priors. For example, 
one may impose uniform priors on the slow-roll param- 
eters, e = (M p 2 ; /2)(V> /V) 2 , r, = M^V" /V), and £ = 
Mp^V'V" /V 2 ), as well as on the number of e-foldings, 
N, rat h er than on n s and dn s /d\n.k (jPeiris fc Eastherl 
l2006allbl ; [E~asther fc Peiriil2006l) . It has been found that, 
as long as one imposes a reasonable lower bound on N . 
N > 30, both approaches yield similar results. 

To constrain n s and dn s /d\nk, we shall use the 



3.1.2. Results: Tilt 

First, we test the inflation models with dn s /dhik = 
0, and negligible gravitational waves. The WMAP 
5-year temperature and polarization data yield n s — 
0.9631q'qj5, which is slightly above the 3-year value 
with a smaller unc ertainty, rt s (3yr) = 0.958 ± 0.016 
(|Spergel et al.1 l2007t ). We shall provide the reason for 
this small upward shift in § 13.1.31 
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The scale-invariant, Harrison-Zel'dovich-Peebles spec- 
trum, n s = 1) is at 2.5 standard deviations away from 
the mean of the likelihood for the WMAP-only analysis. 
The significance increases to 3.1 standard deviations for 
WMAP+BAO+SN. Looking at the two-dimensional con- 
straints that include n s , we find that the most dominant 
correlation that is still left is the correlation between n s 
and fli,h 2 (see Fig. [TJ. The larger il^h 2 is, the smaller 
the second peak becomes, and the larger n s is required 
to compensate it. Also, the larger f^/i 2 is, the larger 
the Silk damping (diffusion damping) becomes, and the 
larger n s is required to compensate it. 

This argument suggests that the constraint on n s 
should improve as we add more information from the 
small-scale CMB measurements that probe the Silk 
damping scales; however, the current data do not im- 
prove the constraint very much yet: n s = 0.960 ± 
0.014 from WMAP+CMB, where "CMB" includes the 
small - scale CMB measurements f r om CBI (Maso net al.l 
2003t iSievers et all 120031 120071: iPearson et alj I2003t 
Readhead et alj |2004|)"7~V S A (Dickinson et alj |2004[ ). 
ACBAR (iKuo et alj l200"l 120071) and BOOMERanG 
(iRuhl et al.ll2003t iMontrov et all 120061 : iPiacentini et all 
2006), all of which go well beyond the WMAP angular 
resolution, so that their small-scale data are statistically 
independent of the WMAP data. 

We find that the small-scale CMB data do not im- 
prove the determination of n s because of their relatively 
large statistical errors. We also find that the calibration 
and beam errors are import ant. Let us examine th is us- 
ing the latest ACBAR data (|Reichardt et alj |2008). The 
WMAP+ACBAR yields 0.964 ± 0.014. When the beam 
error of ACBAR is ignored, we find n s — 0.964 ± 0.013. 
When the calibration error is ignored, we find n s = 
0.962 ± 0.013. Therefore, both the beam and calibration 
error arc important in the error budget of the ACBAR 
data. 

The Big Bang Nucleosynthesis (BBN), combined with 
measurements of the deuterium-to- hydrogen ratio, D/H, 
from quasar absorption systems, has been used exten- 
sively for determining Vihh 2 , independent o f any other 
cosmological parameters (see lStei gman 2007, for a recent 
summary). The p recision of the lates t determination of 
fl b h 2 from BBN (jPettini et al.l 12008( 1 is comparable to 
that of the WMAP data-only analysis. More precise 
measurements of D/H will help reduce the correlation 
between n s and flbh 2 , yielding a better determination of 
n s . 

There is still a bit of correlation left between n s and the 
electron-scattering optical depth, r (see Fig. [1}. While 
a better measurement of r from future WMAP observa- 
tions as well as the Planck satellite should help reduce 
the uncertainty in n s via a better measurement of r, the 
effect of ilbh 2 is much larger. 

We find that the other data sets, such as BAO, SN, 
and the shape of galaxy power spectrum from SDSS or 
2dFGRS, do not improve our constraints on n s , as these 
data sets are not sensitive to fl^h 2 or t; however, this will 
change when we include the running index, dn s /d\nk 
(§ I3.1.3[) and/or the tensor-to-scalar ratio, r (§ I3.2.4[) . 

3.1.3. Results: Running index 

Next, we explore more general models in which a siz- 
able running index may be generated. (We still do 



not include gravitational waves; see § 13.21 for the anal- 
ysis that includes gravitational waves.) We find no 
evidence for dn s /d\nk from WMAP only, -0.090 < 
dn s /d\nk < 0.019 (95% CL), or WMAP+BAO+SN 
-0.068 < dn s /d\nk < 0.012 (95% CL). The improve- 
ment from WMAP-only to WMAP+BAO+SN is only 
modest. 

We find a slight upward shift fro m the 3-year result , 
dn s /d\nk = -0.055to;ol? (68% CL: lSpergel elaT1[2003Tl . 
to the 5-year result, dn s /d\nk = -0.037 + 0.028 (68% 
CL; WMAP only). This is caused by a combination of 
three effects: 

(i) The 3-year number for dn s /dlnk was derived from 
an older analysis pipeline for the temperature data, 
namely, the resolution 3 (instead of 4) pixel-based 
low-Z temperature likelihood and a higher point- 
source amplitude, A ps = 0.017 fiK 2 sr. 

(ii) With two years more integration, we have a better 
signal-to-noise near the third acoustic peak, whose 
amplitude is slightly higher than the 3-year deter- 
mination (|Nolta et al.ll2008r ). 

(iii) With the improved beam model (|Hill et al.ll2008h . 
the temperature power spectrum at I > 200 has 
been raised nearly uniformly by ~ 2.5% (Hil l et aTl 
l200llNolta et al.ll2008h . 

All of these effects go in the same direction, i.e., to in- 
crease the power at high multipoles and reduce a negative 
running index. We find that these factors contribute to 
the upward shift in dn s /d\uk at comparable levels. 

Note that an upward shift in n s for a power-law model, 
0.958 to 0.963 (§ l3~l~2]) . is not subject to (i) because 
the 3-year number for n s was derived from an updated 
analysis pipeline using the resolution 4 low-i code and 
A ps — 0.014 fiK 2 sr. We find that (ii) and (iii) contribute 
to the shift in n s at comparable levels. An upward shift 
in as from the 3-year value, 0.761, to the 5-year value, 
0.796, can be explained similarly. 

We do not find any significant evidence for the run- 
ning index when the WMAP data and small-scale CMB 
data (CBI, VSA, ACBAR07, BOOMERanG) are com- 
bined, -0.1002 < dn s /d\nk < -0.0037 (95% CL), or 
the WMAP data and the lates t results from the ana lysis 
of the complete ACBAR data (|Reichardt et al.ll2008T l are 
combined, -0.082 < dn s /d\nk < 0.015 (95% CL). 

Our best 68% CL constraint from WMAP+BAO+SN 
shows no evidence for the running index, dn s /dhik = 
—0.028 ± 0.020. In order to improve upon the limit on 
dn s /dlnk, one needs to determine n s at small scales, as 
dn s /dlnk is simply given by the difference between n s 's 
measured at two different scales, divided by the loga- 
rithmic separation between two scales. The Lya forest 
provides such information (see § [3 also Table |4|) . 

3.2. Primordial gravitational waves 

3.2.1. Motivation 

The presence of primordial gravitational waves is 
a robust prediction of inflation models, as the same 
mechanism that generated primordial density fluctua- 
tions should also generate primordia l gravitational waves 
(|Grishchuld 119751 (Starobinskvl [l979h . The amplitude of 
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gravitational waves relative to that of density fluctua- 
tions is model-dependent. 

The primordial gravitational waves leave their 
signatures i mprinted on t he C MB temperature 
aniso tro pv (iRubakov et all fl98l iFabbri fc Pollock 



amsotrop 
1983; \) 



Abbott 



Crittenden et ajj 



Wisel 
1991. 



1981 iStarobinskvi EMI; 
as w ell as on polariza- 
tion dBasko fc Polnareyl Il980t iBond fe Efstathioul 
l98l iPolnarevI 119851 : iCrittenden et all 119931 119951 : 



Coulson et al. fl 199 Hi.- 1 The spin-2 nature of gravita- 



tional w aves leads to two types of polarization pattern on 
the sk y (jZaldarriaga fc Seliak 1997t lKamionkowski et alJ 
1997): (i) the curl- free mode (E mode), in which the 
polarization directions are either purely radial or purely 
tangential to hot/cold spots in temperature, and (ii) the 
divergence-free mode (B mode), in which the pattern 
formed by polarization directions around hot/cold spots 
possess non-zero vorticity. 

In the usual gravitational instability picture, in the 
linear regime (before shell crossing of fluid elements) ve- 
locity perturbations can be written in terms of solely a 
gradient of a scalar velocity potential u, v = Vit. This 
means that no vorticity would arise, and therefore no 
B mode polarization can be generated from density or 
velocity perturbations in the linear regime. However, 
primordial gravitational waves can generate both E and 
B mode polarization; thus, the B mode polarization of- 
fers a smoking-gun signatu re for the presence of pri- 
mordial gravitational wave s (jSeliak fc Zaldarriagalll997l : 
iKamionkowski et al.|[l997fh This gives us a strong mo- 
tivation to look for signatures of the primordial gravita- 
tional waves in CMB. 

In Table H we summarize our constraints on the am- 
plitude of gravitational waves, expressed in terms of the 
tensor-to-scalar ratio, r, defined by equation (f^0|) . 

3.2.2. Analysis 

We quantify the amplitude and spectrum of primordial 
gravitational waves in the following form: 



_ k 3 P h (k) 



2vr 2 



^ )(# 



(16) 



where we have ignored a possible scale dependence of 
n t (k), as the current data cannot constrain it. Here, by 
Ph{k) we mean 

(Mk)^(k')) = (2n) 3 P h (k)S 3 (k-k'), (17) 

where hij (k) is the Fourier transform of the tensor metric 
perturbations, gij — a 2 (6ij + h^), which can be further 
decomposed into the two polarization states (+ and x ) 

with the appropriate polarization tensor, e\ - ' , as 

Mk) = Mk)e+(k) + fc x (k)e*(k), (18) 



with the normalization that ejl-e -1- '^ 



e?.e x ' ij ' = 2 and 



e^e*' 13 = 0. Unless there was a parity- violating inter- 
action term such as f{(t>)R l _ L vpaR^ l,( " 7 where /(</>) is an 
arbitrary function of a scalar field, Ry, vpa is the Riemann 
tensor, and W v P a is a dual tensor (|Lue et al.| [l999). both 



21 See, e.g., Watanabc & Komatsu (2006) for the spectrum of 
the primordial gravitational waves itself. 



polarization states are statistically independent with the 
same amplitude, meaning 

(\h + \ 2 ) = (\h x \ 2 )^(\h 2 \) : {K~h* + ) = o. (19) 

This implies that parity-violating correlations such as 
the TB and EB correlations must vanish. We shall ex- 
plore such parity- violating correlations in § |4]in a slightly 
different context. For limits on the difference between 
(|ft + | 2 ) and (\h x \ 2 ) fro m the TB and EB spectra of the 
WMAP 3-year data, see lSaito et"all (|2007f ). 

In any case, this definition suggests that Ph(k) is given 
by P h (k) = i(\h\ 2 ). Notice a factor of 4. This is the defi- 
nition of Ph{k) that we have been us i ng co n sistently since 
the first year release ( Peiris et all 120031: ISpergel et al.l 
120031 12007t iPage et al.ll2007f) . We continue to use this 
definition. 

With this definition of A 2 , (A:) (Eq. [16]), we define the 
tensor-to-scalar ratio, r, at k = fco, as 



A^(fco) 
A^fco)' 



(20) 



where lS.-n(k) is the curvature perturbation spectrum 
given by equation (fT5|) . We shall take fco to be 
0.002 Mpc -1 . In this paper we sometimes call this quan- 
tity loosely the "amplitude of gravitational waves." 

What about the tensor spectral tilt, n t ? In single-field 
inflation models, there exists the so-called consis tency 
relation between r and n* (see lLiddle fc Lvthll2000L for a 



nt 



(21) 



In order to reduce the number of parameters, we continue 
to impose this relation at fc = fco = 0.002 Mpc^ 1 . For 
a discussion on how to impose this constraint in a more 
self-consistent way, see iPeiris fc Eastherl (|2006aD . 

To constrain r, we shall use the WMAP 5-year temper- 
ature and polarization data, the small-scale CMB data, 
and/or BAO and SN distance measurements. 

3.2.3. How WMAP constrains the amplitude of 
gravitational waves 

Let us show how the gravitational wave contribution is 
constrained by the WMAP data (see Fig. [2]) . In this ped- 
agogical analysis, we vary only r and r, while adjusting 
the overall amplitude of fluctuations such that the height 
of the first peak of the temperature power spectrum is 
always held fixed. All the other cosmological parameters 
are fixed at n k = 0, n b h 2 = 0.02265, n c h 2 = 0.1143, 
Ho = 70.1 km s- 1 Mpc" 1 , and n s = 0.960. Note that 
the limit on r from this analysis should not be taken 
as our best limit, as this analysis ignores the correlation 
between r and the other cosmological parameters, espe- 
cially n s . The limit on r from the full analysis will be 
given in § 13.2.41 

(1) (The gray contours in the left panel, and the upper 
right of the right panel of Fig. [2]) The low-1 polar- 
ization data (TE/EE/BB) at I < 10 are unable to 
place meaningful limits on r. A large r can be com- 
pensated by a small r, producing nearly the same 
EE power spectrum at I < 10. (Recall that the 
gravitational waves contribute also to EE.) As a 
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result, r that is as large as 10 is still allowed within 
68% CL. 22 

(2) (The red contours in the left panel, and the lower 
left of the right panel of Fig. [2]) Such a high value 
of r, however, produces too negative a TE corre- 
lation between 30 < I < 150. Therefore, we can 
improve the limit on r significantly - by nearly an 
order of magnitude - by simply using the high-1 TE 
data. The 95% upper bound at this point is still as 
large as r ~ 2. 23 

(3) (The blue contours in the left panel, and the upper 
left of the right panel of Fig. |5J) Finally, the low-l 
temperature data at I < 30 severely limit the excess 
low-^ power due to gravitational waves, bringing 
the upper bound down to r ~ 0.2. Note that this 
bound is about a half of what we actually obtain 
from the full Markov Chain Monte Carlo of the 
WMAP-only analysis, r < 0.43 (95% CL). This 
is because we have fixed n s , and thus ignored the 
correlation between r and n s shown by Fig. [3] 



3.2.4. Results 

Having understood which parts of the temperature and 
polarization spectra constrain r, we obtain the upper 
limit on r from the full exploration of the likelihood space 
using the Markov Chain Monte Carlo. Figure [3] shows 
the 1-d constraint on r as well as the 2-d constraint on r 
and n s , assuming a negligible running index. With the 
WMAP 5-year data alone, we find r < 0.43 (95% CL). 
Since the B-mode contributes little here, and most of the 
information essentially comes from TT and TE, our limit 
on r is highly degenerate with n s , and thus we can obtain 
a better limit on r only when we have a better limit on 
n s . 

When we add BAO and SN data, the limit improves 
significantly to r < 0.22 (95% CL). This is because 
the distance information from BAO and SN reduces the 
uncertainty in the matter density, and thus it helps to 
determine n s better because n s is also degenerate with 
the matter density. This "chain of correlations" helped 
us improve our limit on r significantly from the previ- 
ous results. This limit, r < 0.22 (95% CL), is the best 
limit on r to date. 24 With the new data, we are able 
to g et more stringent t he limits than our earlier analy- 
ses ()Spergel et al.ll2007l) that combined the WMAP data 
with measurements of the galaxy power spectrum and 
found r < 0.30 (95% CL). 

A dramatic reduction in the uncertainty in r has an 
important implication for n s as well. Previously n s > 1 

22 We have perfor med a similar, but different, analysis in § 6.2 
of IPage et all 12007"). In this paper we in clude both the sc alar 
and tensor contributions to EE, whereas in Page ct al. ( 2007) we 
have ignored the tensor contribution to EE and found a somewhat 
tighter limit, r < 4.5 (95% CL), from the low-Z polarization data. 
This is because, when the tensor contribution was ignored, the EE 
polarization could still be used to fix t, whereas in our case r and 
t are fully degenerate when r > 1 (see Fig. [2}, as the EE is also 
dominated by the tensor contribution for such a high value of r. 

23 See IPolnarev et~aTI p008h : IMiller et aEI pOOl ) for a way to 
constrain r from the TE power spectrum alone. 

24 This is the one-dimensional marginalized 95% limit. From 
the joint two-dimensional marginalized distribution of n s and r, 
we find r < 0.27 (95% CL) at n 3 = 0.99. See Fig. [3] 



was within the 95% CL when the gravitational wave con- 
tribution was allowed, owing to the correlation between 
n s and r. Now, we are beginning to disfavor n s > 1 
even when r is non-zero: with WMAP+BAO+SN we 
find -0.0022 < 1 - n s < 0.0589 (95% CL). 25 

However, these stringent limits on r and n s weaken 
to -0.246 < 1 - n s < 0.034 (95% CL) and r < 
0.55 (95% CL) when a sizable running index is allowed. 
The BAO and SN data helped reduce the uncertainty 
in dn s /dhik (Fig. 3]), but not enough to improve on 
the other parameters compared to the WMAP-only con- 
straints. The Lya forest data can improve the limit on 
dn s /dlnk even when r is present (see §[7| also Table HJ. 

3.3. Implications for inflation models 

How do the WMAP 5-year limits on n s and r constrain 
inflationary models? 26 In the context of single-field mod- 
els, one can write down n s and r in terms of the deriva- 
tives of potential, V(<f>), as (lLiddle fc Lvthll2000l ): 



1 



9 V" 



-8M, 



pi 



(22) 



(23) 



where M p i = l/V87rG is the reduced Planck mass, and 
the derivatives are evaluated at the mean value of the 
scalar field at the time that a given scale leaves the hori- 
zon. These equations may be combined to give a relation 
between n, and r: 



r = -(l-n s ) + 



16M 2 V" 



(24) 



This equation indicates that it is the curvature of the 
potential that divides models on the n s -r plane; thus, it 
makes sense to classify inflation models on the basis of 
the sign and magn itude of the curvature of the potential 
(jPeiris et alJ[2003l ). 27 

What is the implication of our bound on r for infla- 
tion models? Equation (f2"4"]l suggests that a large r can 
be generated when the curvature of the potential is pos- 
itive, i.e., V" > 0, at the field value that corresponds 
to the scales probed by the WMAP data. Therefore, it 
is a set of positive curvature models that we can con- 
strain from the limit on r. On the other hand, negative 
curvature models are more diffic ult to constrain fr om r, 
as they tend to predict small r (|Peiris et al.ll2003l ). We 
shall not discuss negative curvature models in this paper: 
many of these models, including those based upon the 
Cole man- Weinberg potential, fit the WMAP data (see, 
e.g., iDvali et al J 11991 iShafi fc Seno"gu3[2006l) . 

Here we shall pick three, simple but representative, 
forms of V(4>) that can produce V" > 28 : 

25 This is is the one-dimensional marginalized 95% limit. From 
the joint two-dimensional marginalized distribution of n s and r, 
we find n s < 1.007 (95% CL) at r = 0.2. See Fig. [3] 

26 For recent s urveys of inflatio n model s in light of the WMA P 
3-year data, see Alabi di fc Lvthl |2006al) : IKinnev et all (2006); 
IMartin fc Ringevall H20061 V 

27 This classification scheme is similar to, but different from, the 
most widely used one, w hich is based upon the field value (s mall- 
field, large-field, hybrid) [|Dodelson et a.l.lfl997l: IKinnevlll99lt) . 

28 These choices are used to sample the space of positive cur- 
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Fig. 2. — How the WMAP temperature and polarization data constrain the tensor-to-scalar ratio, r. (Left) The contours show 68% and 
95% CL. The gray region is derived from the low-Z polarization data (TE/EE/BB at Z < 23) only, the red region from the low-Z polarization 
plus the high-Z TE data at Z < 450, and the blue region from the low-Z polarization, the high-1 TE, and the low-i temperature data at 
I < 32. (Right) The gray curves show (r,r) = (10,0.050), the red curves (r, t) = (1.2,0.075), and the blue curves (r, r) = (0.20,0.080), 
which are combinations of r and r that give the upper edge of the 68% CL contours shown on the left panel. The vertical lines indicate 
the maximum multip oles below which the data are used for each color. The data points with 68% CL errors are the WMAP 5-year 
measurements (Nolta ct al. 2008). (Note that the BB power spectrum at Z ~ 130 is consistent with zero within 95% CL.) 
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approximate the Landau- Ginzburg type of sponta- 
neous symmetry breaking potential, V((f>) oc (</> 2 — 
v 2 ) 2 , in the appropriate limits. 



(b) 



V(4>) oc 
unique feature of 



Fig. 3. — Constraint on the tensor-to-scalar ratio, r, at k = 0.002 Mpc 1 (§ 13. 2. 41 1. No running index is assumed. See Fig. [4] for r 
with the running index. In all panels we show the WMAP-only results in blue and WMAP+BAO+SN in red. (Left) One-dimensional 
marginalized distribution of r, showing the WMAP-only limit, r < 0.43 (95% CL), and WMAP+BAO+SN, r < 0.22 (95% CL). (Middle) 
Joint two-dimensional marginalized distribution (68% and 95% CL), showing a strong correlation between n s and r. (Right) Correlation 
between n s and Q, m h? . The BAO and SN data help to break this correlation which, in turn, reduces the correlation between r and n s , 
resulting in a factor of 2.2 better limit on r. 

(a) Monomial (chaotic-type) potential, V(4>) oc 4> a . 
This form of the potential was proposed by, and 
is b est known fo r, Linde's chaotic inflation mod- 
els (jLindd fT983f ). This model also approximates 
a pseudo Nambu-Goldstone bo son potential (natu- 
ral inflation; iFreese et al]ll990t lAdams et aflll993t ) 
with the negative sign, V{<j>) oc 1 — cos(</>//), 
when (/)/ f <C 1, or with the positive sign, V(<f>) oc 
l + cos(</>//), when </>// ~ l. 29 This model can also 

vature models. Realistic potentials may be much more com- 
plicated: see, for example, IDestri et all 12008ft for the WMAP 
3-year limits on trinomial potentials. Also, the classifica- 
tion scheme based upon derivatives of potentials sheds little 
light on th e models with non-canonical kinetic terms such as 
Zc-inflation ( Armcndariz-Picon ct al. 1999; Garriga & Mukhanov 
119991). ghost inflatio n ( Arkani-Hamcd ct al. 2004|), DBI inflation 
dSilverstein fc Ton e 2004; Alisha hiha et al.N2004i ror infrared-DBI 

(IR-DBI) inflation {Chen 2005b|a]), as the tilt, n s , depends also on t t 

the derivative of the effective speed of sound of a scalar field (for (c) (j) plus vacuum energy, V((f>) = Vq + m (j) /2. 

recent co nstraints on this class of models from the WMAP 3-year These models are known as Linde's hybrid infla- 

data, seem^mmmmm^^Mimm- tion (Lindc 1994) . This model is a « h brid « be _ 

a The positive sign case, V(q>) <x 1 + cosfcp/j), belong s to a ' ' 

negative curvature model when <h/f <C 1. See Savage ct al. (20061) r . . . . r . . r „,,,. _ „ , . 

tor constraints on this class ot models trom the WMAP 3-year data. 



Exponential potential, 

exp[-(0/M p; )\/27^]. A 
this potential is that the dynamics of inflation is 
exactly solvable, and the solution is a power-law 
expansion, a(t) oc t p , rather than an exponen- 
tial one. For this reason this type of model 

is ca lled power-law infla t ion (lAbbott fc Wise! 

11984 iLucchin fc Matarresd |1985[) . They often 
appe ar in models of scala r -tensor theories o f grav- 
ity (lAccetta et al.l Il985l lLa fc SteinhardtJ 119891; 



Futamase fc Maedal 119891 : iSteinhardt fc Accettal 



199ft iKalara et al.Hl99n() . 
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Fig. 4 . — Constraint on the tensor-to-scalar ratio, r, the tilt, n si and the running index, dn s /d In fc, when all of them are allowed to vary 
(§ I3.2.4H . In all panels we show the WMAP-only results in blue and WMAP+BAO+SN in red. (Left) Joint two-dimensional marginalized 
distribution of n s and r at k = 0.002 Mpc~ 1 (68% and 95% CL). (Middle) n s and dn s / d\nk. (Right) dn B /dink and r. We find no 
evidence for the running index. While the inclusion of the running index weakens our constraint on n a and r, the data do not support 
any need for treating the running index as a free parameter: changes in \ 2 between the power-law model and the running model are 
X 2 (running) - x 2 (power - law) ~ —1.8 with and without the tensor modes for WMAP5+BAO+SN, and 1.2 for WMAP5. 



cause the potential combines the chaotic-type (with 
a = 2) with a Higgs-like potential for the second 
field (which is not shown here). This model be- 
haves as if it were a single-field model until the sec- 
ond field terminates inflation when <j) reaches some 
critical value. When </> (2V0) 1 ' 2 /m this model 
is the same as the model (a) with a = 2, although 
one of Linde's motivation was to avoid having such 
a large field value that exceeds M p i. 

These potentials 30 make the following predictions for 



and n s 


as a 


function of their parameters 


(a) r = 


8(1 


-O-Z9 

a / a-\-2 


(b) r = 


8(1 


- n s ) 


(c) r = 


8(1 


n ) 



Here, for (c) we have defined a dimcnsionless variable, 
<f> = mtfi/ (2V0) 1 / 2 ■ This model approaches to the model 
(a) with a = 2 for <p 3> 1, and yields the scale-invariant 
spectrum, n s = 1, when <p = \j\[2. 
We summarize our findings below, and in Fig. [5J 

(a) Assuming that the monomial potentials are valid 
to the end of inflation including the reheating of 
the universe, one can relate n s and r to the num- 
ber of e-folds of inflation, N = \ii(a eTl &/ ciwmap), 
between the expansion factors at the end of infla- 
tion, a en d, and the epoch when the wavelength of 
fluctuations that we probe with WMAP left the 
hori zon during inflation, aw map- The relations 
are (|Liddle fe Lvthll200"oT ): 



4a a + 5 

r = — , 1 — rz., = 

N 2N 



(25) 



30 In the language of § 3.4 in Pciris et al. (2003), the models 
(a) and (b) belong to "small positive curvature models," and the 
model (c) to "large positive curvature models" for <j> -C 1, "small 
positive curvature models" for tj> 2> 1, and "intermediate positive 
curvature models" for cj> ~ 1. 



We take N = 50 an d 60 as a reasonable range 
(jLiddle fc Leach! 120031) . For a = 4, i.e., inflation 
by a massless self-interacting scalar field V(4>) — 
(A/4)</> 4 , we find that both N = 50 and 60 are far 
away from the 95% region, and they are excluded 
convincingly at more than 99% CL. For a = 2, 
i.e., inflation by a massive free scalar field V(<f>) — 
(l/2)m 2 o J ) 2 , the model with N — 50 lies outside of 
the 68% region, whereas the model with N = 60 is 
at the boundary of the 68% region. Therefore, both 
of these models are consistent with the data within 
the 95% CL. While this limit applies to a sin - 
gle massive free field, lEasther fc McAllister! (2006) 
show that a mode l with many massive ax ion fields 
(iV-flation model; [Dimopoulos et al.ll2005f ) can shift 
the predicted n s further away from unity, 



1 



N.f. 



(1 



,s.f. 



1 



(26) 



where "N.f" refers to "N fields," and "s.f." to 
"single f ield," and is a free param eter of the 
model. lEasther fc McAllister! (|2006f ) argue that 
f3 ~ 1/2 is favoured, for which 1 — n s is larger 
than the single-field prediction by as much as 25%. 
The prediction for the tensor-to-s calar ratio, r, is 
the sa me as the single-field case (jAlabidi fc Lvthl 
2006b). Therefore, this model lies outside of the 
95% region for N = 50. As usual, however, 
these monomial potentials can be made a better 
fit to the data by invoking a non-minimal cou- 
pling between the inflaton and gravity, as the non- 
minimal coupling can reduce r to negligible levels 
dKomatsu fc Futamaselll999l: iHwang fc Nohl [19981 : 
iTsuiikawa fc Gumiudpail 120041) . iPiaol (|2006l ) has 
shown that A^-flation models with monomial po- 
tentials, V(<fi) oc 4> a , generically predict n s that is 
smaller than the corresponding single-field predic- 
tions. 

(b) For an exponential potential, r and n s are uniquely 
determined by a single parameter, p, that deter- 
mines a power-law index of the expansion factor, 
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Fig. 5. — Constraint on three representative infla tion models 
whose potential is positively curved, V" > (§ 13.30 . The con- 
tours show the 68% and 95% CL derived from WMAP+BA O+SN 
(Top) The monomial, chaotic-type potential, V(tj>) oc <p a lILinde] 
[1983), with o = 4 (solid) and a = 2 (dashed) for single-field 
models, and a = 2 for mu lti-axion field models with /3 = 1/2 
HEasther fc M cAllister 2006) (dotted). The symbols show the pre- 
dictions from each of these models with the number of e-folds of 
inflation equal to 50 and 60. The A0 4 potential is excluded convinc- 
ingly, the m 2 <f> 2 single-field model lies outside of (at the boundary 
of) the 68% region for N = 50 (60), and the rn 2 <j> 2 multi-axion 
model with N = 50 lies outside of the 95% region. (Middle) 
The exponential potential, V(4>) oc exp[— (<j>/M p i)^/2/p\, which 
leads to a power- law inflat ion, a(t) oc V ( Abbott k Wise 1984; 
Lucchin & Matarresc 1985). All models but p ~ 120 are outside 
of the 68% region. The models with p < 60 are excluded at more 
than 99% CL, and those with p < 70 are outside of the 95% region. 
For multi-field models these limits can be translated into the num- 
ber of fields as p — > npi, where pi is the p-parameter of each field 
dLiddle et al.|[T998T) . The data favour n ~ 120/pi fields. (Bottom) 
The hybrid-type potential, V(<j>) = V + (l/2)m 2 </> 2 = Vn(l + <t> 2 ), 
where = mij)/(2V ) 1 / 2 (ILindelll994h . The models with 4> < 2/3 
drive inflation by the vacuum energy term, Vb, and are disfavoured 
at more than 95% CL, while those with <j> > 1 drive inflation by 
the quadratic term, and are similar to the chaotic type (the left 
panel with a = 2). The transition regime, 2/3 < <j> < 1 are outside 
of the 68% region, but still within the 95% region. 
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We find that p < 60 is excluded at more than 99% 
CL, 60 < p < 70 is within the 99% region but 
outside of the 95% region, and p > 70 is within the 
95% region. The models with p ~ 120 lie on the 
boundary of the 68% region, but other parameters 
are not within the 68% CL. This model can be 
thought of as a single-field inflation with p ^ 1 , or 
multi-field inflation with n fields, eac h having pj ~ 
1 or even pi < 1 (assisted inflation; iLiddle et al.l 
1998). In this context, therefore, one can translate 
the above limits on p into the limits on the number 
of fields. The data favour n ~ 120/pi fields. 

(c) For this model we can divide the parameter space 
into 3 regions, depending upon the value of 4> that 
corresponds to the field value when the wavelength 
of fluctuations that we probe with WMAP left the 
horizon. When ^Cl, the potential is dominated 
by a constant term, which we call "Flat Potential 
Regime." When (j> ^> 1, the potential is indistin- 
guishable from the chaotic- type (model (a)) with 
a — 2. We call this region "Chaotic Inflation- like 
Regime." When <fi ~ 1, the model shows a tran- 
sitional behaviour, and thus we call it "Transition 
Regime." We find that the flat potential regime 
with cf> < 2/3 lies outside of the 95% region. The 
transition regime with 2/3 < cf> < 1 is within the 
95% region, but outside of the 68% region. Finally, 
the chaotic-like regime contains the 68% region. 
Since inflation in this model ends by the second 
field whose dynamics depends on other parameters, 
there is no constraint from the number of e-folds. 

These examples show that the WMAP 5-year data, 
combined with the distance information from BAO and 
SN, begin to disfavour a number of popular inflation 
models. 

3.4. Curvature of the observable universe 
3.4.1. Motivation 

The flatness of the observable universe is one of the 
predictions of conventional inflation models. How much 
curvature can we expect from inflation? The common 
view is that inflation naturally produces the spatial cur- 
vature parameter, Qki on the order of the magnitude of 
quantum fluctuations, i.e., flk ~ 10~ 5 . On the other 
hand, the current limit on flk is of order 10 -2 ; thus, the 
current data are not capable of reaching the level of flk 
that is predicted by the common view. 

Would a detection of fi^ rule out inflation? It is possi- 
ble that the value of Qk is just below our current detec- 
tion limit, even within the context of inflation: inflation 
may not have lasted for so long, and the curvature radius 
of our universe may just be large enough for us not to 
see the evidence for curvature within our measurement 
accuracy, yet. While this sounds like fine-tuning, it is a 
possibility. 

This is something we can test by constraining bet- 
ter. There is also a revived (and growing) interest in 
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measurements of f2fc, as Q k is degenerate with the equa- 
tion of state of dark energy, w. Therefore, a better de- 
termination of f2fc has an important implication for our 
ability to constrain the nature of dark energy. 

3.4.2. Analysis 

Measurements of the CMB power spectrum alone do 
not strongly constrain Q k . More precisely, any experi- 
ments that measure the angular diameter or luminosity 
distance to a single redshift are not able to constrain 
Ofe uniquely, as the distance depends not only on Q k , 
but also on the expansion history of the universe. For 
a universe containing matter and vacuum energy, it is 
essential to combine at least two absolute distance indi- 
cators, or the expansion rates, out to different redshifts, 
in order to constrain the spatial curvature well. Note 
that CMB is also sensitive to J7a via the late-time in- 
tegrated Sachs- Wolfe (ISW) effect, as well as to il m via 
the signatures of gravitational lensing in the CMB power 
spectrum. These properties can be used to reduce the 
corre lation b etween Q k and Q m (IStompor fc Efstathioul 
11999( 1 or fl A (|Ho et all 120081 iGiannantonio et al.ll2008l ). 

It has been pointed out by a number of people (e.g., 
lEisenstein et al.l |2005[ ) that a combination of distance 
measurements from BAO and CMB is a powerful way 
to constrain Q k . One needs more distances, if dark en- 
ergy is not a constant but dynamical. 

In this section, we shall make one important assump- 
tion that the dark energy component is vacuum energy, 
i.e., a cosmological constant. We shall study the case in 
which the equation of state, w, and tt k are varied simul- 
taneously, in § [5j 

3.4.3. Results 

Figure [6] shows the limits on £1a and £l k . While 
the WMAP data alone cannot constrain Qk (see the 
left panel), the WMAP data combined with the HST's 
constraint on Hq tighten the constraint significantly, to 
-0.052 < n k < 0.013 (95% CL). The WMAP data com- 
bined with SN yield ~ 50% better limits, -0.0316 < 
n k < 0.0078 (95% CL), compared to WMAP+HST. Fi- 
nally, the WMAP+BAO yields the smallest statistical 
uncertainty, -0.0170 < il k < 0.0068 (95% CL), which 
is a factor of 2.6 and 1.7 better than WMAP+HST and 
WMAP+SN, respectively. This shows how powerful the 
BAO is in terms of constraining the spatial curvature 
of the universe; however, this statement needs to be 
re-evaluated when dynamical dark energy is considered, 
e.g., w + 1 — 1. We shall come back to this point in §[5] 

Finally, when WMAP, BAO, and SN are combined, 
we find -0.0178 < fl k < 0.0066 (95% CL). As one can 
see from the right panel of Fig. [51 the constraint on £l k 
is totally dominated by that from WMAP+BAO; thus, 
the size of the uncertainty does not change very much 
from WMAP+BAO to WMAP+BAO+SN. Note that 
the above result indicates that we have reached 1.3% 
accuracy (95% CL) in determining Cl k , which is rather 
good. The future BAO surveys at z ~ 3 are expected to 
yield an order of magnitude better determination, i.e., 
0.1% level, of n k (|Knoxll2006H . 

It is instructive to convert our limit on £l k to the lim- 
its on the curvature radius of the universe. As Q k is 
defined as Q k = — fcc 2 /(_ff 2 i? 2 urv ), where i? curv is the 



present-day curvature radius, one can convert the upper 
bounds on fl k into the lower bounds on R curv , as R CUIV — 
(e/ifoVVPTl = i/VW\ h^Gpc. For negatively 
curved universes, we find R CUTV > 37 /i -1 Gpc, whereas 
for positively curved universes, i? CU rv > 22 h~ 1 Gpc. Inci- 
dentally these values are greater than the particle horizon 
at present, 9.7 /i _1 Gpc (computed for the same model). 

Th e 68% limits from the 3-year data (|Spergel et al.l 
l2007h were n k = -0.012 + 0.010 from WMAP-3yr+BAO 
wher e BAO is from the SDSS LRG of lEisenstein et all 
.2005)), and fl k = -0.011 + 0.011 f rom WMAP-3yr+S N 
(where SN is from the SNLS data of lAstier et ail 0006)). 
The 68% limit from WMAP-5yr+BAO+SN (where both 
BAO and SN have more data than for the 3-year analy- 
sis) is fl k = —0.0050^0 0060- A significant improvement 
in the constraint is due to a combination of the better 
WMAP, BAO, and SN data. 

We conclude that, if dark energy is vacuum energy 
(cosmological constant) with w = — 1, we do not find 
any deviation from a spatially flat universe. 

3.4.4. Implications for the duration of inflation 

What does this imply for inflation? Since we do not 
detect any finite curvature radius, inflation had to last 
for a long enough period in order to make the observ- 
able universe sufficiently flat within the observational 
limits. This argument allows us to find a lower bound 
on the total number of e-foldings of the expansion fac- 
tor during inflation, from the begin ning to the end , 
At t = ln(a cn d/abcgm) (see also § 4.1 of IWeinberd feOOS) . 

When the curvature parameter, £l k , is much smaller 
than unity, it evolves with the expansion factor, a, as 
fl k tx a -2 , a 2 , and a during inflation, radiation, and 
matter era, respectively. Therefore, the observed fl k is 
related to fl k at the beginning of inflation as 31 
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(29) 



where <?* is the effective number of relativistic degrees 
of freedom contributing to entropy, z oq is the matter- 
radiation equality redshift, and T en( j and T cq are the re- 
heating temperature of the universe at the end of infla- 
tion 32 and the temperature at the equality epoch, respec- 



tively. To within 10% accuracy, we take 1 + z c 



3200, 



T eq = 0.75 eV, and g* 



3.9. We find 



N t< 
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-*end 1 end 

10 8 GeV 3 200 ' 



(30) 

Here, it is plausible that <?*, e nd ~ 100 in the Stan- 
dard Model of elementary particles, and ~ 200 when 
the supersymmetric partners are included. The differ- 
ence between these two cases gives the error of only 
AA^tot = —0.2, and thus can be ignored. 

31 To simplify our discussion, we ignore the dark energy contri- 
bution, and assume that the universe is dominated matter at the 
present epoch. This leads to a small error in the estimated lower 
bound on Ntot- 

32 For simplicity we assume that reheating occurred as soon as 
inflation ended. 
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Fig. 6. — Joint two-dimensional marginalized constraint on the vacuum energy density, S7a> and the spatial curvature parameter, Qj. 
(5 I3.4.3I I. The contours show the 68% and 95% CL. (Left) The WMAP-only constraint (light blue) compared with WMAP+BAO+SN (pur- 
ple). Note that we have a prior on Q^, Qjy > 0. This figure shows how powerful the extra distance information is for constraining Q k . (Mid- 
dle) A blow-up of the region within the dashed lines in the left panel, showing WMAP-only (light blue), WMAP+HST (gray), WMAP+SN 
(dark blue), and WMAP+BAO (red). The BAO provides the most stringent constraint on Q^. (Right) One-dimensional marginalized 
constraint on O fc from WMAP+HST, WMAP+SN, and WMAP+BAO. We find the best limit, -0.0178 < U k < 0.0066 (95% CL), from 
WMAP+BAO+SN, which is essentially the same as WMAP+BAO. See Fig. I12l for the constraints on Cj, when dark energy is dynamical, 
i.e., w ^ — 1, with time-independent w. Note that neither BAO nor SN alone is able to constrain Q^: they need the WMAP data for lifting 
the degeneracy. Note also that BAO+SN is unable to lift the degeneracy either, as BAO needs the sound horizon size measured by the 
WMAP data. 



The curvature parameter at the beginning of inflation 
must be below of order unity, as inflation would not begin 
otherwise. However, it is plausible that ri^ cgm was not 
too much smaller than 1; otherwise, we have to explain 
why it was so small before inflation, and probably we 
would have to explain it by inflation before inflation. In 
that case iVtot would refer to the sum of the number 
of e-foldings from two periods of inflation. From this 
argument we shall take 0^ ogm ~ 1. 

The reheating temperature can be anywhere between 
1 MeV and 10 16 GeV. It is more likely that it is between 
1 TeV and 10 s GeV for various reasons, but the allowed 
region is still large. If we scale the result to a reasonably 
conservative lower limit on the reheating temperature, 
T cn d ~ 1 TeV, then we find, from our limit on the cur- 
vature of the universe, 



turn fluctuations, P(<p), of free scalar fields in the 
ground state of the Bunch-Davies vacuum, ip, is a 
Gaussian distribution; thus, the probability distribu- 
tion of primordial curvature perturbations (in the co- 
moving gauge), TZ, generated from ip (in the flat 

gauge) as TZ = — [ iif (<ft)/</>n1y> (Mukhanov & C hibisovl 

19811 lHawkind[l98l lStarobinskvl[T98l iGuth fc Pilll982t 
Bardeen et alj[l983h . would also be a Gaussian distribu- 



Here, H(cj)) is the expansion rate during inflation, 
t>o is the mean field, i.e., <fi = <fio + tp. 



iVtot > 36 + m 



1 TeV 



(31) 



A factor of 10 improvement in the upper limit on |r^ cgm | 
will raise this limit by AAT tot =1.2. 

Again, jV to t here refers to the total number of e-foldings 
of inflation. In § 13. 31 we use N = hi(a en< i/awMAp), which 
is the number of e-foldings between the end of inflation 
and the epoch when the wavelength of fluctuations that 
we probe with WMAP left the horizon during inflation. 
Therefore, by definition N is less than N to t- 

3.5. Primordial non-Gaussianity 
3.5.1. Motivation and Background 

In the simplest model of inflation, the distribution of 
primordial fluctuations is close to a Gaussian with ran- 
dom phases. The level of deviation from a Gaussian 
distribution and random phases, called non-Gaussianity, 
predicted by the simplest model of inflation is well below 
the current limit of measurement. Thus, any detection of 
non-Gaussianity would be a significant challenge to the 
currently favored models of the early universe. 

The assumption of Gaussianity is motivated by the 
following view: the probability distribution of quan- 



tion. 

and (pa is me mean neiu, i.e., cp = (po 

This argument suggests that non-Gaussianity can be 
generated when (a) scalar fields are not free, but have 
some interactions, (b) there are non-linear corrections to 
the relation between TZ and ip, and (c) the initial state is 
not in the Bunch-Davies vacuum. 

For (a) one can think of expanding a general scalar 
field potential V(<p) to the cubic order or higher, V(cp) — 
V + Vip + (1/2) V"ip 2 + (1/6) V"'ip 3 + . . .. The cu- 
bic (or higher-ord er) interaction te rms can yield non- 
Gaussianity in ip (|Falk et al.lll993h . When perturba- 
tions in gravitational fields are included, there are many 
more interaction terms that arise from expanding the 
Ricci scalar to the cubic order, with coefficients contain- 

• 3 

ing derivatives of V and <j>o, such as 4>qV' ', <po /H, etc. 
(|Maldacenall2003l ). 

For (b) one can think of this relation, 1Z = 
— [H ((/>)/ <po](p , as the leading-order term of a Taylor 
series expansion of the underlying n on-linear (gauge) 
transf ormation law between 1Z and ip. iSalopek fc Bondl 
(1990) show that, in the single-field models, TZ = 

AttG f*° +v d<j> (d In H/d<j>)~ 1 . Therefore, even if tp is pre- 
cisely Gaussian, TZ can be non-Gaussian due to non-linear 
terms such as ip 2 in a Taylor series expansion of this re- 
lation. One can write this relation in the following form, 
up to second order in TZ, 



TZ = TZ h 



1 fd 2 \nH 



8nG 



d<j) 2 



TZ 



hi 



(32) 



where TZ^ is a linear part of the curvature perturba- 
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tion. We thus find that the second term makes 1Z non- 
Gaussian, even when 7?.l is precisely Gaussian. This 
formula has also been found independently by other re- 
searchers, extended to multi-fie ld cases, and often re- 
ferred to as the U SN formalism" ( Sasaki fc Stewartlll996t 
ILvth et al.H2005HLvth fc Rodriguezll2005h . 

The observers like us, however, do not measure the 
primordial curvature perturbations, 1Z, directly. A more 
observationally relevant quantity is the curvature per- 
turbation during the matter era, At the linear or- 
der these quantities are related by $ = (3/5)7£l (e.g., 
iKodama fe SasakHll984h . but the actual re lation is more 
comp licated at the non-linear order (see iBartolo et al.l 
12004 for a review). In any case, this argument has moti- 
vated our defini ng the 11 local non-linear co upling param- 
eter," /^£ al , as (jKomatsu fe Spergel|[200ll ) 33 

* - *l + f^T^l- (33) 
If we take equation (|32[> . for example, we find /j^£ al = 
-(5/247rG)(9 2 ln J ff/a0 2 ) (|Komatsul [200l Here, we 
have fo llowed the terminology proposed bv lBabich et al.l 
(2004) and called f§ c ^ 1 the "local" parameter, as both 



of order unity to 



.fjg? al by the time we observe it 



sides of equation (|33|) are evaluated at the same location 
in space. (Hence the term, "local".) 

Let us comment on the magnitude of the second term 
in equation (f3"3"| . Since $ ~ 10~ 5 , the second term is 
smaller than the first term by 10 _5 /]y£ al ; thus, the sec- 
ond term is only 0.1% of the first term for /]^£ al — 10 2 . 
As we shall see below, the existing limits on /]^£ al have 
reached this level of Gaussianity already, and thus it is 
clear that we are already talking about a tiny deviation 
from Gaussian fluctuations. This limit is actually better 
than the current limit on the spatial curvature, which is 
only on the order of 1%. Therefore, Gaussianity tests 
offer a stringent test of early universe models. 

In the context of single-field inflation in which the 
scalar field is rolling down the potential slowly, the quan- 
tities, H, V, and 0, are changing slowly. Therefore, one 
generically expects that /^£ al is small, on the order of 
the so-called slow-roll parameters e and r], which are typ- 
ically of order 10 -2 or smaller. In this sense, the single- 
field, slow-roll inflation models are e xpected to result 
in a t iny amoun t of n o n-Gaussianity (ISalopek fc Bond 
19901: l|4lk_et all Il993t iGangui et all 11994 iMaldacena 



inCMB ()Pvne fc Carrolll Il996t iMollerach fc Matarresel 
Il997f ). On smaller angular scales one must include the 
effects of acoustic oscillations of photon-baryon plasma 
by solving the Boltzmann equations. The second-order 
corrections to the B oltzmann equations can also yield 
/^£ al of order unity (jBartolo et al.ll2006t I2007D . 

Any detection of /]^£ al at the level that is currently 
accessible would have a profound implication for the 
physics of inflation. How can a large /]^£ al be gener- 
ated? We essentially need to break either (a) single field, 
or (b) slow-roll. For example, a multi-field model known 
as the curva ton scenario can result in much larger va l- 
ues of f l £f l (|Linde fc Mukhanovlll997t ILvth et al.|[2003l) , 
so can the models with field-dependent (variable) de- 
cay widths for reh eating of the universe after inflation 
(|Dvali et al.ll2004bH ). A more violent, non- linear reheat- 
ing pr ocess called "prehea t ing" can give rise to a large 
f£? al (lEnqvist et al.l l200l I Jokinen fc Mazumdail 120061 : 
iChambers fc Raia ntic 200f). 

Although breaking of slow-roll usually results in 
a premature termination of inflation, it is possible 
to break it temporarily for a brief period, with- 
out terminating inflation, by some features (steps, 
dips, etc) in the shape of the potential. In such 
a scenario, a large non-Gaussianity may be gener- 
ated at a certain limited scale at which the fea- 
ture e xists dKofman et al.lll991b IWang fc Kamionkowskil 
120001: IKomatsu et al.l I2003D . The structure of non- 
Gaussianity from features is much more complex a nd 
model-dependent than f§f l (jChen et al.ll2007l [2008h . 

There is also a possibility that non-Gaussianity can 
be used to test alternatives to inflation. In a collaps- 
ing universe followed by a bounce (e.g., New Ekpyrotic 
scenario), /]^£ al is given by the inverse (as well as inverse- 
squared) of slow-roll parameters; thus, /^£ al as large as 
of order 10 to 1 2 is a fairly generic prediction of this 
class of models dKovama et al.l l2007t iBuchbinder et al. 
1 2008 : ILehners fc Steinhardtl 120081 : iLehners fc Steinhardt 
2008 . 



Using the angular bispectrum, 34 the harmonic trans- 
form of the angula r 3-point correlation function, 
IKomatsu et al.l (|2002f ) have obtained the first obser- 
vational limit on /]^£ al from the COBE 4-year data 
20031; lAcquaviva et aljllli). These contributions from (|Bennett et al.l [199% ). finding -3500 < /^£ al < 2000 



the epoch of inflation are much smaller than those from 
the ubiquitous, second-order cosmological perturbations, 
i.e., the non- linear corrections to the relation between $ 



and TL, which result in / y f of or der unity jlLiguori et al 



2006: Smith fc Zaldarriagall2006f) . See also lBartolo et al 
( 2004]) for a review on this subject 



One can use the cosmological obs ervations, such as 
the CMB data, to constr ain f$f l (jVerde et all 120001 : 
IKomatsu fc Spergel I2001D . While the temperature 
anisotropy, AT/T, is related to $ via the Sachs-Wolfe 
formula as AT/T = — <I> /3 at the linear or der on 
very large angular scales (|Sachs fc Wolfe 19671), there 
are non-linear corrections (non-linear Sachs-Wolfe ef- 
fect, non-linear Integrated Sachs-Wolfe effect, gravita- 
tional lensing, etc) to this relation, which add terms 

33 Note that /^£ al can be related to the quantities discussed 
in earlier, pioneering wo rk: -$ 3 /2 (IGaneui et al.llT99l). - A lr , t ,/2 
iWnne fc Kamionkowski 2003), and -a JVerde et al.H2000h . 



(95% CL). The uncertainty was large due to a relatively 
large beam size of COBE, which allowed us to go only to 
the maximum multipole of Z max = 20. Since the signal- 
to-noise ratio of f l § c ^ y is proportional to Z max , it was ex- 
pected that the WMAP data would yield a factor of ~ 50 
impr ovement over the COBE data (jKomatsu fc Spergel 
I200H) . 

The full bispectrum analysis was not feasible with the 
WMAP data, as the computational cost scales as 
where ./Vpi x is the number of pixels, which is on the order 
of millions for the W MAP data. The "KSW" estimator 
(|Komatsu et alll2005f ) has solved this problem by invent- 
ing a cubic statistic that combines the triangle configura- 
tions of the bispectrum optimally so that it is maximally 

34 For a pedagogical introduction to the bispectrum (3-point 
function) and trispe ctrum (4-point function) and var i ous to pics on 
non-Gaussianity, see IKomatsu (2001); Bart olo et al.l (120041) . 
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sensitive to /^£ al - 35 The computational cost of the KSW 



3 /2 

estimator scales as N ' 1X . We give a detailed description 
of the method that we use in this paper in Appendix [AJ 

We have applied this technique to the WMAP 1- 
year and 3-year data, and found —58 < /Jvff 1 < 134 
(7 max = 265; IKomatsu et "all [200l an d -54 < f l ° c L al < 
114 (Z max = 350; ISpergel et al.l 12007ft . respectively, at 
95% CL. Creminelli et al. performed an independent 
analysis of the WMAP data an d found similar limits : 
-27 < f l £f l < 121 (l max = 335; ICreminelli et~all[2n0^ 
and -36 < /^£ al < 100 (7 max = 370; ICreminelli etai] 
2007) for the 1 -year and 3-year data, respectively. These 
constraints are slightly better than the WMAP team's, 
as their estimator for /]^£ al was improved from the orig- 
inal KSW estimator. 

While these constraints are obtained from the KSW- 
like fast bispectrum statistics, many groups have 
used the WMAP data to measure /]y£ a ' using vari- 
ous other statistics, such as the Minkowski function- 
al (IKomatsu et al.ll200l lSpergel et alJl2007t iGott et all 
2007; Hikagc et al.l 12008ft . real -space 3-point fu nction 
( Gaztana ga fc Wagd|2003t IChen fc Szapudll2005ft , inte- 
grated bispectrum (jCabella et al.l l2006ft. 2-1 cum ulant 
correlator power spectrum (IChen fc Szap udi 2006), local 
curvature (ICabella et al.ll2004ft. and s pherical Mexican 
hat wavele ufMukheriee fc Wang||200l) . The sub-orbital 
CMB experiment s have also yield ed co nstraints on 
f}^ al : MAXIMA (ISantos et al.1120 03 7 ) . VSA (jSmith et all 
2004). Archeops (ICurto et al.ll2007ft . and BOOMERanG 
(|De Troia et al.ll2007ft . 

We stress that it is important to use different sta- 
tistical tools to measure /jyff if any signal is found, 
as different tools are sensitive to different systemat- 
ics. Th e analytical predict ions for the Minkowski func- 
tional (Hika ge et al.ll2006ft and the angular trispectrum 
(the harmonic transform of the angular 4-point corre- 
lation function; lOkamoto fc Hu 2002; iKogo fc Komatsul 
2006) as a function of /]y£ al are available now. Studies 
on the forms of the trispectrum from inflation models 
have just begun, and some important insights have been 
obtained (iBoubekeur fc Ijtgi^l [Huang fc Shh]|2006l 
Bvrnes et all 120061: ISeerv fc Lidsevl l2007t ISeerv et al 
20071 lArroia fc Kovamal 12008ft . It is now understood 



that the trispectrum is at least as important as the bis- 
pectrum in discriminating inflation models: some mod- 
els do not produce any bispectra but produce significant 
trispectra, and other models produce similar amplitudes 
of the bispectra but produce very different trispectra 
(|Huang fc Shidl27MlBuchbinder et al1l2008ft . 

In this pa per we shall use the estim ator that further im- 
proves upon ICreminelli et alJ (|2006ft by correcting an in- 
adve rtent numerical err or of a factor of 2 in t heir d eriva- 
tion (|Yadav et al.ll2008ft . lYadav fc Wandeltl (|2008ft used 
this estimator to measure /]^£ al from the WMAP 3-year 
data. We shall also use the Minkowski functionals to find 
a limit on /J^£ al . 

In addition to /]^£ al , we shall also estimate the "equilat- 

35 Since the angular bispectrum is the harmonic transform of the 
angular three-point function, it forms a triangle in the harmonic 
space. While there are many possible triangles, the "squeezed tri- 
angles," in which the two wave vectors are long and one is short, 
are most sensitive to /J^£ al (Babich et "aT1l2004h . 



eral non-linear coupling parameter," /^j* 1 , which charac- 
terizes the amplitude of the three-point function (i.e., the 
bispectrum) of the equilateral configurations, in which 
the lengths of all the three wave vectors forming a trian- 
gle in Fourier space are equal. This parameter is useful 
and highly complementary to the local one: while /J^£ al 
characterizes mainly the amplitude of the bispectrum of 
the squeezed configurations, in which two wave vectors 
are large and nearly equal and the other wave vector is 
small, and thus it is fairly insensitive to the equilateral 
configurations, is mainly sensitive to the equilat- 

eral configurations with little sensitivity to the squeezed 
configurations. In other words, it is possible that one 
may detect /]^£ al without any detection of and 
vice versa. 

These two parameters cover a fairly large class of 
models. For example, iv?* can be generated from in- 
flation models in which the scalar field takes on the 
non-standard (non-canonical) kinetic form, such as C — 
P(X,<fi), where X = (dtp) 2 . In this class of models, 
the effective sound speed of <fi can be smaller than the 
speed of light, c 2 s = [l + 2X(d 2 P/dX 2 )/{dP/dX)}- 1 < 1. 

While the sign of is negative for the DBI infla- 

tion, /J^j" ~ — l/c 2 < in the limit of c s <C 1, it 
can be positive or n e gative for more general models 
(ISeerv fc Lidsevl 120051 fChen et all I2007T iCheung et all 
120081: iLi et al.ll2008ft . Such models can be realized in 
the context of String Theory via the non-canonical ki- 
netic action called the Dirac-Born-Infeld (DBI) form 
(jAlishahiha et all [2004) , and in the context of an in- 
frared modification of gravity cal led the ghost conden- 
sation (|Arkani-Hamed et aLll2004ft . 

The observational limits on f^ 11 have been obtained 
from the WMAP 1-year and 3-year data as —366 < 
f^f < 238 (Z max = 405; ICreminelli et~aT1 l2006ft and 
-256 < f^f < 332 (Z max = 475: ICreminelli et al.ll2007ft . 
respectively. 

There are other forms, too. Warm inflation might 
produce a different fo rm of /nl|Moss fc Xiond 2007; 
IMoss fc Grahaml 12007ft . Also, the presence of particles 
at the beginning of inflation, i.e., a departure of the ini- 
tial state of quantum fluctuations from the Bunch-Davics 
vacuum, can result in an enhanced n on-Gaussianity in 
the "flattened" triangle configurations (jChen et al.ll2007l : 
iHolman fc Toll cv 2008) . We do not consider these forms 
of non-Gaussianity in this paper. 

In this paper we do not discuss the non-Gaussian 
signatures that cannot be characterized by /]^£ al , 



cquil 
NL ' 



(the point source bispectrum ampli- 



tude). There have been many studies on non-Gaussian 
signatures in the WMAP data in various forms 
dChiang et alJ I2003L [20071: iNaselskv et~afl l2007t IParl 



2004 1 : Tde Oliveira-Costa et al.l 120041: iTegmark et al 



2004. 



20031: ILarson fc Wandeltl 120041: lEriksen et all 12004 
Erik sen et al.1 120041: lEriksen et al 
Copi et all 120041: ISchwarz et al.l 120041 
2006, 2007:~lGordon et al 



2007a 
Copi et al 



20051: iBielewicz et all 



2005t l.laffe et al.l 120051 120061: iVielva et al 



Cruz et alJ 120051 I2006T l2007blal iCavon et al.l 



200' 
200f 



Bridges et aij|2007t IWiaux et al l 20081: iRath et al.l 200' 
Land fc Magueiidl2005bllaL l2007l:lRakic fc Schwargl^ 



IPark et al.ll2007HBernui et al.ll2007HrIaiian fc Souradeepl 
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2003t lHaiian et al l 120051: lHaiian fe Souradeeul 120061: 
Prunet et alJ 120051 : lHansen et all l2004t lHansen et all 
2004 ) . many of which are related to the large-scale 
features at I < 20. We expect these features to be 
present in the WMAP 5-year temperature map, as the 
structure of CMB anisotropy in the WMAP data on 
such large angular scales has not changed very much 
since the 3-year data. 

3.5.2. Analysis 

The largest concern in measuring primordial non- 
Gaussianity from the CMB data is the potential contami- 
nation from the Galactic diffuse foreground emission. To 
test how much the results would be affected by this, we 
measure /n l parameters from the raw temperature maps 
as well as from the foreground-reduced maps. 

We shall mainly use the KQ75 mask, the new mask 
that is recommended for tests of Gaussianity ( Go ld et al.l 
l2008f ). The important diffe rence between the ne w mask 
and the previous KpO mask (jBennett et al.l[20~03cl ) is that 
the new mask is defined by the difference between the K 
band map and the Internal Linear Combination (ILC) 
map, and that between the Q band and ILC. Therefore, 
the CMB signal was absent when the mask was defined, 
which removes any concerns regarding a potential bias in 
the distribution of CMB on the masked sky. 36 

To carry out tests of Gaussianity, one should use the 
KQ75 mask, which is slightly more conservative than 
KpO, as the KQ75 mask cuts slightly more sky: we retain 
71.8% of the sky with KQ75, while 76.5% with KpO. To 
see how sensitive we are to the details of the mask, we 
also tried KpO as well as the new mask that is recom- 
mended for the power spectrum analysis, KQ85, which 
retains 81.7% of the sky. The previous mask that corre- 
sponds to KQ85 is the Kp2 mask, which retains 84.6% 
of the sky. 

In addition, we use the KQ75pl mask, which replaces 
the point source mask of KQ75 with the one that does 
not mask the sources identified in the WMAP K-band 
data. Our point source selection at K band removes more 
sources and sky in regions with higher CMB flux. We 
estimate the amplitude of this bias by using the KQ75pl 
mask which does not use any WMAP data for the point 
source identification. The small change in /]^£ al shows 
that this is a small bias. 

The unresolved extra-galacti c point source s also 
contribute to the bis pectr um (iRefreeier et alJ 20001 : 
Komatsu fe Spergell 120011 : lArgueso et alJ 120031 : 



Serra fe Cooravl 20081 ). and they can bias our esti 



mates of primordial non-Gaussianity parameters such 
as /]££ al and /jy^ 1 ■ We estimate the bias by measuring 
/]v£ al an d /tvl'' horn Monte Carlo simulations of point 
sources, and list them as A/^£ al and A/^ n in Table E 

36 Previously, the KpO mask was defined by the K band map, 
which contains CMB as well as the foreground emission. By cutting 
bright pixels in the K band map, it could be possible to cut also 
the bright CMB pixels, introducing the negative skewness in the 
distribution of CMB. Since we did not include isolated "islands" 
on the high Galactic latitudes, some of which could be bright CMB 
spots, in the final mask when we defined the KpO mask, the skew- 
ness bias mentioned above should not be as large as one would 
expect, if any. Nevertheless, with the new definition of mask, the 
masked maps are free from th is type of b i as. Fo r more details on 
the definition of the mask, see Gold et al. (200f|). 



TABLE 5 

Clean-map estimates and the corresponding 68% 
intervals of the local form of primordial 
non-Gaussianity, /]^£ al , THE point source bispectrum 

AMPLITUDE, b arc (iN UNITS OF 10 -5 /xK 3 sr 2 ), AND 

Monte- Carlo estimates of bias due to point sources, 

A flocal 
^■JNL 



Band 


Mask 


^raax 


f local 
J NL 


A f local 
J NL 


&src 


V+W 


KQ85 


400 


50 + 29 


1 + 2 


0.26+ 1.5 


v+w 


KQ85 


500 


61 ± 26 


2.5 + 1.5 


0.05 ± 0.50 


V+W 


KQ85 


600 


68 + 31 


3 + 2 


0.53 + 0.28 


v+w 


KQ85 


700 


67 + 31 


3.5 + 2 


0.34 + 0.20 


v+w 


KpO 


500 


61 ± 26 


2.5 + 1.5 




v+w 


KQ75pl a 


500 


53 + 28 


4 + 2 




v+w 


KQ75 


400 


47 + 32 


3 + 2 


-0.50+ 1.7 


v+w 


KQ75 


500 


55 ± 30 


4 + 2 


0.15 + 0.51 


v+w 


KQ75 


600 


61 ± 36 


4 + 2 


0.53 + 0.30 


v+w 


KQ75 


700 


58 ± 36 


5 + 2 


0.38 + 0.21 



"This mask replaces the point-source mask in KQ75 with the 
one that does not mask the sources identified in the WMAP 
K-band data 



TABLE 6 

Null tests, frequency dependence, and 
raw-map estimates of the local form of 

PRIMORDIAL NON-GAUSSIANITY, /J v ° cal 
/max = 500 



I NL 



FOR 



Band Foreground Mask 



flocal 
■I NL 



Q-W 


Raw 


KQ75 


-0.53 + 0.22 


v-w 


Raw 


KQ75 


-0.31 + 0.23 


Q-W 


Clean 


KQ75 


0.10 + 0.22 


V-W 


Clean 


KQ75 


0.06 + 0.23 


Q 


Raw 


KQ75pl a 


-42 ± 45 


V 


Raw 


KQ75pl 


38 + 34 


w 


Raw 


KQ75pl 


43 + 33 


Q 


Raw 


KQ75 


-42 ± 48 


V 


Raw 


KQ75 


41 ± 35 


w 


Raw 


KQ75 


46 + 35 


Q 


Clean 


KQ75pl 


9 + 45 


V 


Clean 


KQ75pl 


47+34 


w 


Clean 


KQ75pl 


60 + 33 


Q 


Clean 


KQ75 


10 ± 48 


V 


Clean 


KQ75 


50 + 35 


w 


Clean 


KQ75 


62 + 35 


v+w 


Raw 


KQ85 


9 + 26 


v+w 


Raw 


KpO 


48 + 26 


v+w 


Raw 


KQ75pl 


41 ± 28 


v+w 


Raw 


KQ75 


43 + 30 



"This mask replaces the point-source mask in 
KQ75 with the one that does not mask the sources 
identified in the WMAP K-band data 

and [7] As the errors in these estimates of the bias 
are limited by the number of Monte Carlo realizations 
(which is 300) , one may obtain a better estimate of the 
bias using more realizations. 

We give a detailed description of our estimators for 
/]v£ al > /tvl' 1 , and b src , the amplitude of the point source 
bispectrum, as well as of Monte Carlo simulations in Ap- 
pendix [A] 



3.5.3. Results: Bispectrum 

In Table [5] we show our meas urement of /w f al from 
the template-cleaned V+W map (|Gold et al.lfe008h with 
4 different masks, KQ85, KpO, KQ75pl, and KQ75, in 
the increasing order of the size of the mask. For KQ85 
and KQ75 we show the results from different maximum 
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TABLE 7 
Clean-map estimates and the 
corresponding 68% intervals of the 
equilateral form of primordial 

NON-GAUSSIANITY, fjf^", AND 
MONTE-CARLO ESTIMATES OF BIAS DUE TO 



POINT SOURCES. 



A fcquil 
•> NL 



Band 


Mask 


^raax 


rcquil 
JNL 


A requil 


V+W 


KQ75 


400 


77 ± 146 


9 + 7 


v+w 


KQ75 


500 


78 ± 125 


14 + 6 


V+W 


KQ75 


600 


71 ± 108 


27 + 5 


v+w 


KQ75 


700 


73 ± 101 


22 + 4 



TABLE 8 

Point source bispectrum amplitude, b s 

^max — 900 



Band 


Foreground 


Mask 


b src [10- 5 fiK 3 sr 2 ] 


Q 


Raw 


KQ75pl a 


11.1 ± 1.3 


V 


Raw 


KQ75pl 


0.83 + 0.31 


w 


Raw 


KQ75pl 


0.16 + 0.24 


v+w 


Raw 


KQ75pl 


0.28 + 0.16 


Q 


Raw 


KQ75 


6.0 + 1.3 


V 


Raw 


KQ75 


0.43 + 0.31 


w 


Raw 


KQ75 


0.12 + 0.24 


v+w 


Raw 


KQ75 


0.14 + 0.16 


v+w 


Raw 


KQ85 


0.20 + 0.15 


Q 


Clean 


KQ75pl 


8.7+ 1.3 


V 


Clean 


KQ75pl 


0.75 + 0.31 


w 


Clean 


KQ75pl 


0.16 + 0.24 


v+w 


Clean 


KQ75pl 


0.28 + 0.16 


Q 


Clean 


KQ75 


4.3 + 1.3 


V 


Clean 


KQ75 


0.36 + 0.31 


w 


Clean 


KQ75 


0.13 + 0.24 


v+w 


Clean 


KQ75 


0.14 + 0.16 


v+w 


Clean 


KQ85 


0.13 + 0.15 



"This mask replaces the point-source mask in KQ75 
with the one that does not mask the sources identified 
in the WMAP K-band data 



TABLE 9 



X 



FUNCTIONALS FOR THE 
TEMPLATE-CLEANED V+W MAP. THE 

RESULTS FROM THE AREA, CONTOUR 
LENGTH, AND EULER CHARACTERISTICS 
ARE COMBINED 



N side 


XwMAp/ doi 


F(> XwMAp) 


256 


51.5/45 


0.241 


128 


40.0/45 


0.660 


64 


54.2/45 


0.167 


32 


46.8/45 


0.361 


16 


44.7/45 


0.396 


8 


61.3/45 


0.104 



multipoles used in the analysis, l m&K — 400, 500, 600, and 
700. The WMAP 5-year temperature data are limited by 
cosmic variance to I ~ 500. 

We find that both KQ85 and KpO for l max = 500 show 
evidence for fftf 1 > at more than 95% CL, 9 < ffif 1 < 
113 (95% CL), before the point source bias correction, 
and 6.5 < /^£ al < 110.5 (95% CL) after the correction. 
For a higher Z max , i max = 700, we still find evidence for 

flocal 
JNL 
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Fig. 7. — Minkowski functionals from the WMAP 5-year data, 
measured from the template-cleaned V+W map at N s i^ c = 128 
(28' pixels) outside of the KQ75 mask. From the top to bottom 
panels, we show the Euler characteristics (also known as the genus), 
the contour length, and the cumulative surface area, as a function 
of the threshold (the number of cr's of hot and cold spots), v = 
AT/o-q. (Left) The data (symbols) are fully consistent with the 
mean and dispersion of Gaussian realizations that include CMB 
and noise. The gray bands show the 68% intervals of Gaussian 
realizations. (Right) The residuals between the WMAP data and 
the mean of the Gaussian realizations. Note that the residuals arc 

flocal 



-57 + 60 (68% CL). From N siic 
(68% CL). 



64 we find / 



local 
NL 



; ± 69 



> 0, 1.5 < f£ c L al < 125.5 (95% CL), after the 



correction. 37 

This evidence is, however, reduced when we use larger 
masks, KQ75pl and KQ75. For the latter we find 
~ 5 < /jvL al < 115 ( 95% CL ) before the source bias 
correction, and -9 < /^£ al < 111 (95% CL) after the 
correction, which we take as our best estimate. This 
estimate improves upon our previous esti mate from the 
3-yea r data, -54 < f l ° c L al < 114 (95% CL: iSpergel elall 
120071 for Z max = 350), by cutting much of the allowed 
region for < 0. To test whether the evidence for 

f^L > can a l so be seen with the KQ75 mask, we need 
more years of WMAP observations. 

Let us study the effect of mask further. We find that 
the central value of /j^£ al (without the source correction) 
changes from 61 for KQ85 to 55 for KQ75 at l max = 500. 
Is this change expected? To study this we have com- 

37 The uncertainty for i ma x > 500 is slightly larger than that for 
(max = 500 due to a small sub-optimality of the estimator of /J^£ al 
llYadav et al.ll2008h . 
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puted /^£ al from each of the Monte Carlo realizations 
using KQ85 and KQ75. We find the r.m.s. scatter of 
_ = 13j 12> 15, and 15, for 

; ma:E = 400, 500, 600, and 700. Therefore, the change 
in /]^ al measured from the WMAP data is consistent 
with a statistical fluctuation. For the other masks at 



It. 



500 we find {{f§f yKp() - /^£ alifQ75 ) 



2\l/2 
/MC 



2\l/2 
/MC 



= 9.7 



= 4.0. 



and ((f£i alKQ75pl - f l £T lKQ75 ) 

In Table[6]we summarize the results from various tests. 
As a null test, we have measured /]^£ al from the differ- 
ence maps such as Q— W and V— W, which are sensitive 
to non-Gaussianity in noise and the residual foreground 
emission. Since the difference maps do not contain the 
CMB signal, which is a source of a large cosmic variance 
in the estimation of /]^£ al , the errors in the estimated 
/Jv£ al are much smaller. Before the foreground cleaning 
( "Raw" in the second column) we see negative values of 
/]v£ a1 ' w hi cn is consistent with the foreground emission 
having positively skewed temperature distribution and 
/Jv£ al > mainly generating negative skewness. We do 
not find any significant signal of /]^£ al at more than 99% 
CL for raw maps, or at more than 68% CL for cleaned 
maps, which indicates that the temperature maps are 
quite clean outside of the KQ75 mask. 

From the results presented in Tabled we find that the 
raw-map results yield more scatter in /]^£ al estimated 
from various data combinations than the clean-map re- 
sults. From these studies we conclude that the clean-map 
results are robust against the data combinations, as long 
as we use only the V and W band data. 

In Table [7] we show the equilateral bispectrum, , 
from the template-cleaned V+W map with the KQ75 
mask. We find that the point source bias is much more 
significant for /jy^' : we detect the bias in f^^ 11 at more 
than the 5a level for / max = 600 and 700. After correcting 
for the bias, we find -151 < ffif < 253 (95% CL; 
'max = 700) as our best estimate. Our estimate improves 
upon the previous on e, —256 < f^ 1 < 332 (95% CL; 
ICreminelli et al.l [20061 for l max = 475), by reducing the 
allowed region from both above and below by a similar 
amount. 

Finally, the bispectrum from very high multipoles, e.g., 
'max = 900, can be used to estimate the amplitude of 
residual point source contamination. One can use this in- 
formation to check for a consistency between the estimate 
of the residual point sources from the power spectrum 
and that from the bispectrum. In Table [8] we list our 
estimates of b src . The raw maps and cleaned maps yield 
somewhat different values, indicating a possible leakage 
from the diffuse foreground to an estimate of b src . Our 
best estimate in the Q band is b src = 4.3 ± 1.3 //K sr 2 
(68% CL). See lNolta etaTI (|2008l ) for the comparison be- 
tween b src , C ps , and the point-source counts. 

Incidentally, we also list b src from the KQ75pl mask, 
whose source mask is exactly the same as we used for the 
first-year analysis. We find b src = 8.7 ± 1.3 /iK 3 sr 2 in 
the Q band, which is in an excellent agreement with the 
first-y ear result, b src = 9.5 ±4.4 [iK 3 sr 2 (jKomatsu et al.1 
[2003^ . 

3.5.4. Results: Minkowski Functionals 



For the analysis of the Mi nkowski functionals, w e fol- 
low the meth o d des cribed in iKomatsu et al.l (|2003l ) and 
ISpergel et"aTl (|2007h . In Fig. U\ we show a ll of the three 
Minkowski functionals dGott et al.1 Il990l: iMecke et al. 
|1994t lSchmalzing fc BuchertHl997l : ISchmalzing fc Gorski 
119981 : IWinitzki &: Kosowskvl[l9981 ) that one can define on 
a two-dimensional sphere: the cumulative surface area 
(bottom), the contour length (middle), and the Euler 
characteristics (which is also known as the genus; top), 
as a function the "threshold," is, which is the number of 
cr's of hot and cold spots, defined by 



AT 

0o ' 



(34) 



where o~o is the standard deviation of the temperature 
data (which includes both signal and noise) at a given 
resolution of the map that one works with. We compare 
the Minkowski functionals measured from the WMAP 
data with the mean and dispersion of Gaussian realiza- 
tions that include CMB signal and noise. We use the 
KQ75 mask and the V+W-band map. 

While Fig. [7] shows the results at resolution 7 (N si( i c = 
128), we have carried out Gaussianity tests using the 
Minkowski functionals at six different resolutions from 
resolution 3 (A^ido = 8) to resolution 8 (N S id c = 256). We 
find no evidence for departures from Gaussianity at any 
resolutions, as summarized in Table [9j in this table we 
list the values of % 2 of the Minkowski functionals relative 
to the Gaussian predictions: 



XWMAP — 



WMAP 



(^m>], 



F 3 — (F J ) 

WMAP \ r sim/ 



(35) 



where F^y MAp and F s l im are the ith Minkowski function- 
als measured from the WMAP data and Gaussian sim- 
ulations, respectively, the angular bracket denotes the 
average over realizations, and EJJ ' is the covariance ma- 
trix estimated from the simulations. We use 15 different 
thresholds from v = —3.5 to v = +3.5, as indicated by 
the symbols in Fig. and thus the number of degrees of 
freedom in the fit is 15 x 3 = 45. We show the values 
of Xwmap an d the degrees of freedom in the second col- 
umn, and the probability of having \ 2 that is larger than 
the measured value, F(> Xwmap)' m the third column. 
The smallest probability is 0.1 (at N S id e = 8), and thus 
we conclude that the Minkowski functionals measured 
from the WMAP 5-year data are fully consistent with 
Gaussianity. 

What do these results imply for /]^£ al ? We find that 
the absence of non-Gaussianity at N S id c = 128 and 64 
gives the 68% limits on f° c L al as /^£ al = -57 ± 60 and 
-68 ± 69, respectively. The 95% limit from N sidc = 128 
is —178 < /j^£ al < 64. The errors are larger than those 
from the bispectrum analysis given in § 13.5.31 by a factor 
of two, which is partly because we have not used the 
Minkowski functional at all six resolutions to constrain 
/jy£ ■ For a combined ana lysis of the WMAP 3-year 
data, see lHikage et al.l (|2008l ). 

It is intriguing that the Minkowski functionals pre- 
fer a negative value of /]v£ al , Inl^ ~ — 60, whereas 
the bispectrum prefers a positive value, /]^£ al ~ 60. In 
the limit that non-Gaussianity is weak, the Minkowski 
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functionals are sensitive to three "skewness param- 
eters:" (1) ((AT) 3 ), (2) ((AT) 2 [d 2 (AT)]), and (3) 
([<9(AT)] 2 [d 2 (AT)]), all of which can be written in 
terms of the weighted sum of the bispectrum; thus, the 
Minkowski functionals are sensit ive to some selected con- 
figurations of the bispectrum (jHikage et al.l l2Q06h . It 
would be important to study where an apparent "ten- 
sion" between the Minkowski functionals and the KSW 
estimator comes from. This example shows how impor- 
tant it is to use different statistical tools to identify the 
origin of non-Gaussian signals on the sky. 

3.6. Adiabaticity of primordial fluctuations 
3.6.1. Motivation 

"Adiabaticity" of primordial fluctuations offers impor- 
tant tests of inflation as well as clues to the origin of mat- 
ter in the universe. The negative correlation between the 
temperature and E-mode polarization (TE) at / ~ 100 
is a ge neric signature of adiabatic super-horizon fluctu- 
ations (|Spergel fe Zaldarriagal 119971 : iPeiris et alJ l2003f ). 
The improved measurement of the TE power spectrum 
as well as the temperature power spectrum from the 
WMAP 5-year data, combined with the distance infor- 
mation from BAO and SN, now provide tight limits on 
deviations of primordial fluctuations from adiabaticity. 

Adiabaticity may be defined loosely as the following re- 
lation between fluctuations in radiation density and those 
in matter density: 



3Sp r 8p Tl 
Ap r p m 



(36) 



This version 38 of the condition guarantees that the en- 
tropy density (dominated by radiation, s r oc pV^) per 
matter particle is unperturbed, i.e., S(s r /n m ) = 0. 

There are two situations in which the adiabatic con- 
dition may be satisfied: (1) there is only one degree of 
freedom in the system, e.g., both radiation and matter 
were created from decay products of a single scalar field 
that was solely responsible for generating fluctuations, 
and (2) matter and radiation were in thermal equilib- 
rium before any non-zero conserving quantum number 
(such as baryon number minus l epton number, B — L) 
was cr eated (e.g.. IWeinbergll2004f ). 

Therefore, detection of any non-adiabatic fluctuations, 
i.e., any deviation from the adiabatic condition (Eq. [31)]). 
would imply that there were multiple scalar fields during 
inflation, and cither matter (baryon or dark matter) was 
never in thermal equilibrium with radiation, or a non- 
zero conserving quantum number associated with matter 
was created well before the era of thermal equilibrium. 
In any case, the detection of non-adiabatic fluctuations 
between matter and radiation has a profound implica- 
tion for the physics of inflation and, perhaps more im- 
portantly, the origin of matter. 

For example, axions, a good candidate for dark matter, 
generate non-adiabatic fluctuations between dark matter 

38 A more general relation is Sp x /p x = $py/py, where x and y 
refer to some energy components. Using the energy conservation 
equation, p x = —3H(1 + w x )p x (where w x is the equation of state 
for the component x), one can recover Eq. 13611 . as w r = 1/3 and 
w m = 0. For a recent discussion on this topic, see, e.g.. IWeinberel 



and photons, as axion density fluctuations could be pro- 
duced during inflation independent of curvature pertur- 
bations (which were generated from inflaton fields, and 
responsible for CMB anisotropies that we observe today) , 
and were not in th ermal equilibrium with radiation in the 
early universe (see lKolb fc TurneriH990HSikividl200cl for 
reviews). We can therefore place stringent limits on the 
properties of axions by looking at a signature of devia- 
tion from the adiabatic relation in the CMB temperature 
and polarization anisotropies. 

In this paper we focus on the non-adiabatic perturba- 
tions between cold dark matter (CDM) and CMB pho- 
tons. Non-adiabatic perturbations between baryons and 
photons are exactly the same as those between CDM and 
photons, up to an overall constant; thus, we shall not con- 
sider them separately in this paper. For neutrinos and 
photons we consider only adiabatic perturbations. In 
other words, we consider only three standard neutrino 
species (i.e., no sterile neutrinos), and assume that the 
neutrinos were in thermal equilibrium before the lepton 
number was generated. 

The basic idea behind this study is not new, and 
adiabaticity has been constrained extensively using the 
WMAP data since the first year release, including gen- 
eral (phenomeno logical) studies w i thout references to 
specific models j£eiris et alJ 120031: ICrottv et all 120031: 
Bucher et alll2004 Moodlcv ct al. 2004; Lazaridcs 'et all 
2004 iKurki-Suonio et all 120051: IBeltran et al 
Dunklev et al.l 120051: iBean et al l l2006t iTrott" 



2005 



2007 



Keskitalo et al.ll2007h . as well as c onstraints on specific 
models such as double infla t ion (|Silk fc Turner! 119871: 
Polarski fe Starobinskvi Il992t I1994D. axio n (IWeinbere 



19781: iWilczekl 119781: iSeckel fc Turnerl [19851 : lLindej[l98^ 
. .Turner fc Wilczeld Il99lf). and cu r vaton 
dLvth fc Wandd 120031: jMoroi fc Takahashil l200l 12001 : 
iBartolo fc Liddla l2002h . all of which can be con- 
strained from the limits on non-adiabatic fluctua- 



iBeltran et al. 2004 


; Lazaridesl 12005: Parkinson et al. 


120051 IBeltran et al. 


120051 120071: Kawasaki & Sckiguchi 



120071) . 

We shall use the WMAP 5-year data, combined with 
the distance information from BAO and SN, to place 
more stringent limits on two types of non-adiabatic 
CDM fluctuations: (i) axion-type, and (ii) curvaton- 
ty pe. Our study given below is similar to the one 
by iKawasaki fc Sekiguchil ([2001 for the WMAP 3 -year 
data. 

3.6.2. Analysis 

We define the non-adiabatic, or entropic, perturbation 
between the CDM and photons, 5 Ci7 , as 

6p c 35py 
Pc 4p 7 ' 

and report on the limits on the ratio of the power 
spectrum of 5 Cj7 , Ps(k), to the curvature perturbation, 
P-n(k), at a give n pivot wavenumber, fco, given by (e.g., 
IBean et al.ll2006D 



^c,7 — 



(37) 



a(k 



Ps(kg) 

l-a(Ad) ~ Pn(k )' 
We shall take k to be 0.002 Mpc -1 . 



(38) 
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Fig. 8. — Constraint on the axion entropy perturbation fraction, ao (§ 13,6,31 . In all panels we show the WMAP-only results in blue 
and WMAP+BAO+SN in red. (Left) One-dimensional marginalized constraint on ao, showing WMAP-only and WMAP+BAO+SN. 
(Middle) Joint two-dimensional marginalized constraint (68% and 95% CL), showing the correlation between ao and n s for WMAP-only 
and WMAP+BAO+SN. (Right) Correlation between n B and Q m h 2 . The BAO and SN data help to reduce this correlation which, in turn, 
reduces correlation between qq and n s , resulting in a factor of 2.2 better limit on ao. 
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Fig. 9. — Constraint on the curvaton entropy perturbation fraction, o_i (§ I3.6,4j l, In all panels we show the WMAP-only results in blue 
and WMAP+BAO+SN in red. (Left) One-dimensional marginalized constraint on a_i, showing WMAP-only and WMAP+BAO+SN. 
(Middle) Joint two-dimensional marginalized constraint (68% and 95% CL), showing the correlation between c«_i and n s for WMAP-only 
and WMAP+BAO+SN. (Right) Correlation between n s and Q m h 2 . The BAO and SN data help to reduce this correlation which, in turn, 
reduces correlation between ot—\ and n s , resulting in a factor of 2.7 better limit on ct—\. These properties are similar to those of the axion 
dark matter presented in Fig. \E\ 

While a parametrizes the ratio of the entropy power 
spectrum to the curvature power spectrum, it may be 
more informative to quantify "how much the adiabatic 
relation (Eq. [35]) can be violated." To quantify this, 
we introduce the adiabaticity deviation parameter, 5 a di, 
given by 

Spc/Pc - 3fy 7 /(4/3 7 ) 



There cou ld be a signi f icant correlation between 
S c .y and K ilLangloisl [l999t iLanglois k Riazueld [2000l : 



iGordon et al.ll2001 ). We take this into account by intro- 
ducing the cross-correlation coefficient, as 



f3(k ) 



Ps,Tz(k ) 



r(c,7) 
adi 



(39) 



(42) 



z[5p c /p c + 3<Sp 7 /(4p 7 )] 

which can be used to say, "the deviation from the adi- 
abatic relation between dark matter and photons must 

be less than 100<5^' J 7 ' ) %-" The numerator is just the def- 
inition of the entropy perturbation, <S C , 7 , whereas the 
denominator is given by 

I ( ^£ + 3( ^°7 



= 3^ + 0(5). 



Pc 4p 7 

Therefore, we find, up to the first order in <S/7?., 
i(c,7) _ S _ V" 



(40) 



where Ps,ii(k) is the cross-correlation power spectrum. 

Here, we have a negative sign on the left hand side 
because of the following reason. In our notation that 
we used for Gaussianity analysis, the sign convention of 
the curvature perturbation is such that it gives temper- 
ature anisotropy on large scales (the Sachs-Wolfe limit) 
as AT/T = -(1/5)K. On the other hand, those who 
investigate correlations between S and 1Z usually use an 
opposite sign convention for the curvature perturbation, 
such that the temperature anisotropy is given by 



adi 



3TZ 



(41) 



AT 1 - 2 „ 

= -K - -S 

T 5 5 



(43) 



for a«l. 



on large angular scales (|Lang lois 1999]), where TZ = —TZ, 
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and define the correlation coefficient by 



0(ko) = 



(44) 



Therefore, in order to use the same sign convention for (3 
as most of the previous work, we shall use Eq. (|42p . and 
call P = +1 "totally correlated," and (3 = — 1 "totally 
anti-correlated," entropy perturbations. 

It is also useful to understand how the correlation or 
anti-correlation affects the CMB power spectrum at low 
multipoles. By squaring equation (|43[) and taking the 
average, we obtain 

((AT 1 ) 2 

J^2 



Therefore, the "correlation," (3 > 0, reduces the temper- 
ature power spectrum on low multipoles, whereas the 
"anti-correlation," (3 < 0, increases the power. This 
point will become important when we interpret our re- 
sults: namely, models with (3 < will result in a positiv e 
correlation between a and n s ([Gordon fc Lewis! 12003). 
Note that this property is similar to that of the tensor 
mode: as the tensor mode adds a significant power only 
to I < 50, the tensor-to-scalar ratio, r, is degenerate with 
n s (see Fig. El . 

Finally, we specify the power spectrum of S as a pure 
power-law, 



P s (k) oc ¥' 



PsMk) oc fc(«'+« s )/2-4 7 (46) 



in analogy to the curvature power spectrum, P-jz(k) oc 
fc™ s_4 . Note that j3 does not depend on k for this choice 
°f Ps,n{k)- With this parametrization, it is straight- 
forward to compute the angular power spectra of the 
temperature and polarization of CMB. 

In this paper we shall pay attention to two limit- 
ing cases that are physically motivated: totally un- 
correl ated (0 = 0) entropy pert urbations - axion- 
tvpe (ISeckel fc Turner! [19851 : iLmdel [1985L ILTndel Il99l 
iTurner fc Wilczekl Il991h . and totally anti-correlated 
(P = —1) entropy pertu r bations - curvaton- 
tvpe dLinde fc Mukhanwl I1997K iLvth fc Wands] 120031 : 
iMoroi fc Takahashill200lL 120021: iBartolo fc Liddld l2002h . 
Then, we shall use ao to denote a for (3 = and a_i for 
= -l. 

3.6.3. Results: Implications for Axion 

First, let us consider the axion case in which S and 1Z 
are totally uncorrelated, i.e., (3 = 0. This case represents 
the entropy perturbation between photons and axions, 
with axions accounting for some fraction of dark matter 
in the universe. For simplicity, we take the axion pertur- 
bations to be scale invariant, i.e., m = 1. In Appendix 151 
we show that this choice corresponds to taking one of the 
slow-roll parameters, e, to be less than 10 -2 , or adding 
a tiny amount of gravitational waves, r <C 0.1, which 
justifies our ignoring gravitational waves in the analysis. 

The left panel of Fig. [8] shows that we do not find any 
evidence for the axion entropy perturbations. The lim- 
its are a < 0.16 (95% CL) and a < 0.072 (95% CL) 
for the WMAP-only analysis and WMAP+BAO+SN, re- 
spectively. The latter limit is the most stringent to date, 
from which we find the adiabaticity deviation parameter 



of S c a ^ < 0.089 (Eq. [39]); thus, we conclude that the ax- 
ion dark matter and photons should obey the adiabatic 
relation (Eq. [35]) to 8.9%, at the 95% CL. 

We find that n s and ao are strongly degenerate (see 
the middle panel of Fig. [8]). It is easy to understand 
the direction of correlation. As the entropy perturba- 
tion with a scale invariant spectrum adds power to the 
temperature anisotropy on large angular scales only, the 
curvature perturbation tries to compensate it by reduc- 
ing power on large scales with a larger tilt, n s . However, 
since a larger n s produces too much power on small an- 
gular scales, the fitting tries to increase f^/i 2 to suppress 
the second peak, and reduce f2 c /i 2 to suppress the third 
peak. Overall, £l m h 2 needs to be reduced to compensate 
an increase in n s , as shown in the right panel of Fig. [8] 

Adding the distance information from the BAO and SN 
helps to break the correlation between fl m h 2 and n s by 
constraining fl m h 2 , independent of n s . Therefore, with 
WMAP+BAO+SN we find an impressive, factor of 2.2 
improvement in the constraint on ag. 

What does this imply for the axions? It has been 
shown that the limit on the axion entropy perturbation 
can be used to place a constraint on the energy scale of 
inflation which, in turn, leads to a stringent constraint 
on the tensor-to-scalar ratio, r jK ain 2006; Be ltran et al.l 
I2007t ISikivi3l2008t iKawasaki fc Sekiguchifl2007h . 

In Appendix [B] we study a particular axion cosmol- 
ogy called the "misalignment angle scenario," in which 
the Pecci-Quinn symmetry breaking occurred during in- 
flation, and was never restored after inflation. In other 
words, we assume that the Pecci-Quinn symmetry break- 
ing scale set by the axion decay constant, f a , which has 
been constrained t o be greater tha n 10 10 GeV from the 
supernova 1987A (jYao et al.ll200l . is at least greater 
than the reheating temperature of the universe after 
inflation. This is a rather reasonable assumption, as 
the reheating temperature is usually taken to be as low 
as 10 8 GeV i n order to avoid overproduction of un- 
wanted relics dPagels fc Primackl Il982t ICoughlan et al.l 
119831 : lEllis et al.lll986h . Such a low reheating tempera- 
ture is natural also because a coupling between inflaton 
and matter had to be weak; otherwise, it would termi- 
nate inflation prematurely. 

There is another constraint. The Hubble parame- 
ter during inflation needs to be smaller than f a , i.e., 
Hi n f ^ fa', otherwise, the Pecci-Quinn symmetry would 
be re stored by quantum fluctuations (|Lyth fc StewartI 
[1992). 

In this scenario axions acquired quantum fluctuations 
during inflation, in the same way that inflaton fields 
would acquire fluctuations. These fluctuations were then 
converted to mass density fluctuations when axions ac- 
quired mass at the QCD phase transition at ~ 200 MeV. 
We observe a signature of the axion mass density fluctu- 
ations via CDM-photon entropy perturbations imprinted 
in the CMB temperature and polarization anisotropies. 

We find that the tensor-to-scalar ratio r, the axion 
density, fi a , CDM density, f2 c , the phase of the Pecci- 
Quinn field within our observable universe, 9 a , and ao, 
are related as (for an alternative expression that has f a 
left instead of 9 a , see Eq. [B7j ) 



4.7 x 10~ 12 ftt r h 



,10/7 



7 



12/7 



fie 



2/7 



1 — ao 



, (47) 



WMAP 5-year Cosmological Interpretation 



25 



< 



(0.99 x 10- 13 )a 
^ 0/7 7 i2/7 



fie 
0„ 



2/7 



(48) 



where 7 < 1 is a "dilution factor" representing the 
amount by which the axion density parameter, fl a h 2 , 
would have been diluted due to a potential late-time en- 
tropy production by, e.g., decay of some (unspecified) 
heavy particles, between 200 MeV and the epoch of 
nucleosynthesis, 1 MeV. Here, we have used the limit 
on the CDM density parameter, fl c h 2 , from the ax- 
ion entropy perturbation model that we consider here, 
Q c h 2 = 0.1052+q qq7q, as well as the observational fact 
that ao -C 1. 

With our limit, a < 0.072 (95% CL), we find a limit 
on r within this scenario as 



r < 



6.6 x 10 



-15 



? l°/ 7 7 12/7 



2/7 



(49) 



Therefore, in order for the axion dark matter scenario 
that we have considered here to be compatible with 
il c ~ il a and the limits on the non-adiabaticity and Q c h 2 , 
the energy scale of inflation should be low, and hence the 
gravitational waves are predicted to be negligible, unless 
the axion density was diluted severely by a late-time en- 
tropy production, 7 ~ 0.8 x 10 -7 (for 9 a ~ 1), the axion 
phase (or the misalignment angle) within our observable 
universe was close to zero, a ~3x 10~ 9 (for 7 ~ 1), or 
both 7 and 9 a were close to zero with lesser degree. All 
of these possibilities would give r ~ 0.01, a value that 
could be barely detectable in the foreseeable future. One 
can also reverse Eq. (|49|) to obtain 



a, 
a 



< 



3.0 x 10 



-39 



Oh" 



0.01 



7/2 



(50) 



Therefore, the axion density would be negligible for the 
detectable r, unless 9 a or 7 or both are tuned to be small. 

Whether such an extreme production of entropy is 
highly unlikely, or such a tiny angle is an undesirable 
fine-tuning, can be debated. In any case, it is clear that 
the cosmological observations, such as the CDM density, 
entropy perturbations, and gravitational waves, can be 
used to place a rather stringent limit on the axion cos- 
mology based upon the misalignment scenario, one of the 
most popular scenarios for axions to become a dominant 
dark matter component in the universe. 

3.6.4. Results: Implications for Curvaton 

Next, let us consider one of the curvaton models in 
which S and TZ are to t ally anti-correlated, i.e., = — 1 
dLvth fe Wandl [200l i Moroi fe Takahashl [200ll [200l 
iBartolo fe Lidda T2002'l . One can also write S as iS = 



BIZ = -BIZ where B > 0; thus, B 2 = a_i/(l - cli). 39 
We take the spectral index of the curvaton entropy per- 
turbation, m, to be the same as that of the adiabatic 
perturbation, n s , i.e., n s — m. 

The left panel of Figure [9] shows that we do not 
find any evidence for the curvaton entropy perturba- 
tions, either. The limits, a_i < 0.011 (95% CL) and 
a_i < 0.0041 (95% CL) for the WMAP-only analysis 



and WMAP+BAO+SN, respectively, are more than a 
factor of 10 better than those for the axion perturba- 
tions. The WMAP-only li mit is better tha n the previ- 
ous limit by a factor of 4 (|Bean et al.ll2006T ). From the 
WMAP+BAO+SN limit, we find the adiabaticity devia- 
tion parameter of 5^?' < 0.021 (Eq. [35]); thus, we con- 
clude that the curvaton dark matter and photons should 
obey the adiabatic relation (Eq. [36]) to 2.1% at the 95% 
CL. 

Once again, adding the distance information from the 
BAO and SN helps to reduce the correlation between n s 
and Q m h 2 (see the right panel of Fig. [9]), and reduces the 
correlation between n s and a_i. The directions in which 
these parameters are degenerate are similar to those for 
the axion case (see Fig. [8]), as the entropy perturba- 
tion with = — 1 also increases the CMB temperature 
power spectrum on large angular scales, as we described 
in § I3~Q 

What is the implication for this type of curvaton sce- 
nario, in which = — 1? This scenario would arise when 
CDM was created from the decay prod ucts of the curva - 
ton field. One then finds a prediction (jLvth et al.ll2003f l 



1 



9 



1 Pcurvaton/ptotal 
Pcurvaton/ Ptotal 



(51) 



where p C urvaton and ptotai are the curvaton density and to- 
tal density at the curvaton decay, respectively. Note that 
there would be no entropy perturbation if curvaton dom- 
inated the energy density of the universe completely at 
the decay. The reason is simple: in such a case all of the 
curvaton perturbation would become the adiabatic per- 
turbation, so would the CDM perturbation. Our limit, 
a_i < 0.0041 (95% CL), indicates that p CU rvaton/ptotai is 
close to unity, which simplifies the relation (Eq. [51]) to 
give 



Pcurvaton 
Ptotai 



1 - 



/Q '-! = 1 _ r(c,7) 
^ adi ' 



(52) 



Note that it is the adiabaticity deviation parameter given 
by Eq. (J39j that gives the deviation of Pcurvaton/ptotai 
from unity. From this result we find 



1 > Pcurvaton > q gg 
Ptotai 



(95% CL). 



(53) 



As we mentioned in § 13.5. 1[ the curvaton scenario is 
capable of producting t he local form of non-Ga ussianity, 
and /^£ al is given by (]Lvth fc Rodriguez! [20051 and ref- 
erences therein 40 ) 



Hocal 
JNL 



5ptotal 



4p c 



5p c 



6ptotal 



(54) 



which gives —1.25 < f£ c L 21 (curvaton) < —1.21, for 
a_i < 0.0041 (95% CL). While we need to add ad- 
ditional contributions from post-inflationary, non-linear 
gravitational perturbations of order unity to this value in 
order to compare with what we measure from CMB, the 
limit from the curvaton entropy perturbation is consis- 
tent with the limit from the measured /]^£ al (see § I3.5.3jl . 

However, should the future data reveal ,f j°r al ^> 1, then 
either this scenario would be ruled out (Beltran 2008; 



39 T his variable, B, is the same as B used in Gordon & Lewis 
(2003), including the sign convention. 



40 Note that the sign convention of /Jj?r al m Lvth & Ro drigue3 



(2005) 



is such that / local 



NL ,WMAP 



INL 
rlocal 

JNL , theirs 
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iLi et al.ll2008b( ), or the curvaton dark matter must have 
been in thermal equilibrium with photons. 

For the other possibiliti es, including pos s ible b aryon 
entropy perturbations, see lGordon fc Lewis] (|2003f ). 



4. PROBING PARITY VIOLATION OF THE UNIVERSE: TB 
AND EB CORRELATION 

4.1. Motivation 

Since the temperature and E-mode polarization are 
parity-even and the B-mode polarization is parity-odd, 
the TB and EB correlation s should vanish in a uni- 
verse that conserves parity dKamionkowski et al.l 
j Kamionkowski et al.l Il997t ISeliak fc Zaldarriagal 
IZaldarriaga fc Sell aid 1 19971 ). For this reason the TB and 



1997 



1997 



EB correlations are usually used to check for systematics, 
and not widely used as a cosmological probe. 

However, parity is violated in th e weak interactions 
(|Lee fc Yangl 119561 IWu et al.lfl957T l. Why can't parity 
be violated at cosmological scales? 

Polarization of photons offers a powerful way of prob- 
ing the cosmologica l parity violation, or the "cosm olog- 
ical birefringence" (|Lue et al.l Il999t ICarrol]|[T99i ). Let 
us consider a parity-violating interaction term in the 
Lagrangian such as the Chern-Simons term, Cqs = 
-{l/2)p a A f} F a ' :t , where F a(3 and Ap are the usual elec- 
tromagnetic tensor and vector potential, respectively, 
F af} = {l/2)e aP ^F lllJ is a du al tensor, and p a is an ar- 
bitrary timelike four- vector. 41 ICarroll et al.l (|1990f ) have 
shown that the Chern-Simons term makes two polariza- 
tion states of photons propagate with different group ve- 
locities, causing the polarization plane to rotate by an 
angle Aa. 

What would p a be? We may take this to be a derivative 
of a light scalar field, p a — 2(d a (f>)/M, where M is some 
unspecified energy scale. In this case, the rotation angle 
is given by Aa = f ^<ft/M = (A(j))/M (|Carroll et all 
[1995 1 Carrolill998l : ILiu et "aT1[2006l IXraeFaTJ[200l) . Such 
a field might have something to do with dark energy, 
for example. We are therefore looking at a potential 
parity-violating interaction between the visible section 
(i.e., photons) and dark sector (i.e., dark energy). 

Such an unusual rotation of polarization vectors has 
been co nstrained by observations of radio galaxies and 
quasars (|Carrollll 19981 ): one of the best data sets available 
today at a single redshift is 3C9 at z = 2.012, which gives 
a limit on the rotation angle, Aa = 2° ± 3° (68% CL). 
There are about 10 measurements between z = 0.425 
and z = 2.012, whose average is Aa = -0.6° ±1.5° (68% 
CL). 

The rotation of polarization plane converts the E-modc 
polarization to the B-mode. As a result, B modes can be 
produced from E modes even if inflation did not produce 
much B modes. This is similar to the gravitational lens- 
ing effect, which also produ ces B modes from E modes 
(jZaldarriaga fc Seljaldll99 8:). but there is an important 
difference: the lensing does not violate parity, but this 
interaction does. As a result, the lensing does not yield 
non-zero TB or EB, but this interaction yields both TB 
and EB. 



41 S ee fLeporal (119981) ; iKlinkhamcr ( 20001; lAdanTfc' Klmkhamcr 
pOOlh for studies on a spacelike p a , including its signatures in 
CMB. 



We shall constrain Aa between the reionization epoch 
z ~ 10 and the present epoch, as well as Aa between 
the decoupling epoch, z ~ 1090, and the present epoch, 
using the TB and EB spectra that we measure from the 
WMAP 5-year data. 

4.2. Analysis 

Before we proceed, we should remember that the mag- 
nitude of polarization rotation angle, Aa, depends on 
the path length over which photons exp erienced a parity- 
violating interaction. As pointed out bv lLiu et al.1 (2006), 
this leads to the polarization angle that depends on I. We 
can divide this ^-dependence in two regimes: 

• Z < 20: the pola rization signal was generated dur- 
ing reionization (|Zaldarriagalll997T ). We are sensi- 
tive only to the polarization rotation between the 
reionization epoch and present epoch. 

• Z > 20: the polarization signal was generated at 
the decoupling epoch. We are sensitive to the po- 
larization rotation between the decoupling epoch 
and present epoch; thus, we have the largest path 
length in this case. 

Below, we shall explore two cases separately. Note that 
we shall use only the polarization spectra: TE, TB, 
EE, BB, and EB, and do not use the temperature spec- 
trum, as the temperature spectrum is not affected by the 
parity-violating interaction. 

Moreover, for the analysis at I < 23 we only vary the 
polarization angle Aa, and the optical depth, r, and fix 
the other parameters at flk = 0, Slbh 2 = 0.02265, il c h 2 = 
0.1143, H = 70.1 km s" 1 Mpc" 1 , and n s = 0.960. At 
each value of r, we re-adjust the overall normalization of 
power spectra such that the first peak of the temperature 
spectrum is held fixed. For the analysis at I > 24 we fix 
r at 0.085, and vary only Aa, as there is no correlation 
between Aa and r at high multipoles. We ignore EE, 
BB and EB at I > 24, as they are much noisier than TE 
and TB and thus do not add much information. 

When the polarization plane is rotated by Aa, the 
intrinsic (primordial) TE, EE, and BB spectra are 
converted into TE, TB, EE, BB , and EB spectra as 
ijLue et al.lll999llFeng et al.ll2005l ) 



TE.obs 



-cT E 



gTB.obs _ (jTE 



; cos(2Aa), 
~sin(2Aa), 



(55) 
(56) 



c EE,obs = C EE cos 2 (2Aa) + C BB ^^Aa), (57) 



C, 



BB.obs 
I 



Cf E sin 2 (2 Aa) + Cf B cos 2 (2 Aa) , (58) 
(Cp* --f'pMnUAni. (r,.) i 



jjEB.obs ^ (-,BB 



where C;'s are the primordial power spectra in the ab- 
sence of parity violation, while C° bs 's are what we would 
observe in the presence of parity violation. To simplify 
the problem and maximize our sensitivity to a potential 
signal of Aa, we ignore the primordial BB, and use only 
a reduced set: 



a 



TE,obs 



TB.obs 



EE,obs 



= C; cos(2Aa), 
-Cj E sin(2Aa), 
■-Cf E cos 2 (2Aa), 



(60) 
(61) 
(62) 
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Fig. 10. — C onstraint on the polarization rotation angle, Aa, due to a parity-violating interaction that rotates the polarization angle 
of CMB (§H3). We have used the polarization spectra (TE/TB/EE/BB/EB at I < 23, and TE/TB at I > 24), and did not use the TT 
power spectrum. (Left) One-dimensional marginalized constraint on Aa in units of degrees. The dark blue, light blue, and red curves show 
the limits from the low-Z (2 < I < 23), high-i (24 < I < 450), and combined (2 < I < 450) analysis of the polarization data, respectively. 
(Right) Joint two-dimensional marginalized constraint on r and Aa (68% and 95% CL). The bigger contours are from the low-i analysis, 
while the smaller ones are from the combined analysis. The vertical dotted line shows the best-fitting optical depth in the absence of parity 
violation (t = 0.086), whereas the horizontal dotted line shows Aq = to guide eyes. 



BB.obs 



■a 



EE , 



2 (2Aa) 



EB,obs 



= ^sin(4Aa) 



(63) 
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Therefore, TB and EB will be produced via the "leak- 
age" from TE and EE. Note that E and B are totally 

l i_ j ■ ,i ■ (r^EB,obs\i r~iEE ,obs ^BB ,obs 

correlated m this case: (C l ) — C t ' L l 

Several groups have constrained Aa from the WMAP 
3-year data as well as from the BOOMERanG data 
dFeng et alJ 120061: iLiu et all 120061: [Kosteleckv k Mewesl 
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200J). All Tut 



. Cabella et al.l l2007t Kia et al , 
Liu et al. I (120061) assume that Aa is constant at all mul- 
tipoles, which is acceptable when they consider the TB 
and EB data at I > 20, i.e., the BOOMERanG data and 
high-/ WMAP data. However, it requires care when one 
considers the low-/ WMAP data. Moreover, all of the 
authors used a Gaussian form of the likelihood function 
for Ci, which is again acceptable at high multipoles, but 
it is inaccurate at low multipoles. 

For the 5-year data release we have added capabilities 
of computing the likelihood of TB and EB spectra at 
low multipoles, 2 < I < 23, exactly, as well as that of 
TB spectrum at high multipoles, 24 < I < 450, using the 
MASTER (pseudo-C;) algorithm. We shall use this code 
to obtain the limit on Aa from the 5-year WMAP polar- 
ization data. For the low-Z polarization we use the Ka, 
Q, and V band data, whereas for the high-/ polarization 
we use the Q and V band data. 

4.3. Results 

Fig. [T0l shows our limit on Aa between (i) the reioniza- 
tion epoch and present epoch from the low-/ polarization 
data (dark blue), (ii) between the decoupling epoch and 
present epoch from the high-/ polarization data (light 
blue), and (iii) combined constraints from the low-/ and 
high-/ data assuming a constant Aa across the entire 
multipole range (red). We find no evidence for parity- 
violating interactions: the 95% CL (68% CL) limits are 
-22.2° < Aa < 7.2° (Aa = -7.5° ± 7.3°) for (i), 
-5.5° < Aa < 3.1° (Aa = -1.2° ± 2.2°) for (ii), and 
-5.9° < Aa < 2.4° (Aa = -1.7° ± 2.1°) for (iii). 



The previous 95% CL (68% CL) limits on Aa 
are largely based upon the high-/ TB and EB data 
from the WMAP 3-year data and/or BOOMERanG: 
-13. 7° < Aa < 1.9° (Aa = -6.0° ± 4.0°) (|Feng et all 
120061). -25° < Aa < 2° (Aa = -12° ± 7°) 
(jKosteleckv k Mewesl 12001 ■ -8.5 ° < Aa < 3.5° (Aa = 
-2.5° ± 3.0°) dCabella et aH 12001. an d -13.8° < Aa < 
1.4° (Aa = -6.2° ± 3.8°) (|Xia et all 12001 . Our lim- 
its from the WMAP 5-year data are tighter than the 
previous ones by a factor of 1.5 to 2, and already compa- 
rable to those from the polarization data of radio galax- 
ies and quasars (see § 14. 1| . Note that the radio galaxies 
and quasars measure the rotation of polarization between 
up to z — 2 and the present epoch, whereas our lim- 
its measure the rotation between the decoupling epoch, 
z ~ 1090, and the present epoch. 

These results show that the TB and EB polarization 
data can provide interesting limits on parity- violating in- 
teraction terms. The future data will be able to place 
more stringent limits (|Xia et al.l 120081 ). In particular, 
adding the Ka and W band data to the high-/ polariza- 
tion should improve our limit significantly. 

5. DARK ENERGY 

5.1. Motivation 

Dark energy is one of the most mysterious obser- 
vations in physics today. The issue is the following: 
when the luminosity d i stances out to Type l a super- 
novae (|Riess et al.lll998t iPerlmutter et alJ ll999) and the 
angular diameter dist ances measured fr om the BAO 
(lEisenstein et alJ l2005h as well as CMB ([Bennett et al.l 
l2003bf l are put together in the context of homogeneous 
and isotropic cosmological models, one cannot fit these 
distances without having an accelerated expansion of the 
universe today. A straightforward interpretation of this 
result is that we need an additional energy component 
in the universe that has a large negative pressure, which 
causes the expansion to accelerate. 

However, we do not know much about dark energy. 
A study of review articles written over the p ast twenty 
years reveals a growing circle of ignorance ( Weinberg 
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1989f ICarroll et all Il992t ISahni fe Starobinskvl 
200a iPadmanabhanl 120031: iPeebles fc Ratral 1200a 
Padmanabhanl l2005t ICopeland et al.l I2006D : physicists 
first struggled to understand why the cosmological 
constant or vacuum energy term was so close to zero, 
then to understand why it was non-zero. Cosmolo- 
gists then explored the possibility that dark energy 
was d ynam i cal, e .g ., in a fo r m of s ome light scalar 
field fcbrdl Il987t iWetterichl [l988t iRatra fc Peebles! 
19881 [Peebles fc Ratral Il988t iFuiTi fc Nishiokal 119901: 
Chiba et al.lll997t ICaldwell et al.l 119981: ICopeland et al l 
19981: iFerreira fc Jovcd Il998t IZlatev et alJ Il999f l. 



Recently, there has been significant interest in modifica- 
tions to General Relativity, in the co ntext of exp l aining 
the acceleration of the universe (jDvali et all 120001: 
IDeffavet et aT1l2002l) . 

Currently, the properties of dark energy are mainly 
constrained by the distance information. There are other 
promising ways of finding dark energy independent of dis- 
tances: the expansion rate of the universe at higher (z > 
0.5) redshifts, the Integrated Sachs- Wolfe (ISW) effect, 
and a slow-down of the growth of large-scale structure in 
the universe due to dark energy. While these tools are 
powerful in principle, the current data are not accurate 
enough to distinguish between the effects of dark energy 
and spatial curvature of the universe, owing to the degen- 
eracy between them fe.g.. iNesseris fc Perivolaropoulosl 
2008; IHo et al.ll2008l: IGiannantonio et al.ll2008T) ~ 

Indeed, the properties of dark energy, such as the 
present-day density and its evolution, e.g., the equation 
of state parameter w, are degenerate with the spatial 
curvature of the universe, £lk- In this section we shall 
explore both flat and curved universes when we report 
on our limits on the dark energy properties. 

In § 15.21 and § 15.31 we explore constraints on time- 
independent (i.e., constant) equation of state, w, assum- 
ing flat (ilfc = 0) and curved (fife + 1 0) geometries, re- 
spectively. In § 15.41 we introduce a set of "WMAP dis- 
tance priors," and use them to explore a wider range of 
model space that has time-dependent equation of state, 
w = w(z). Throughout § 15.41 we use the distance infor- 
mation only to constrain the properties of dark energy. 
We thus assume the standard homogeneous and isotropic 
Friedmann-Lemaitre- Robertson- Walker universe, and do 
not consider modifications of gravity or local inhomo- 
geneity, as the distance information alone cannot dis- 
criminate between these models and the accelerated ex- 
pansion due to dark energy. Finally, in S I5.5l we introduce 
a "WMAP normalization prior." 



5.2. 



Constant equation of state: Flat universe 



What are we doing by assuming a flat universe, when 
we constrain the dark energy equation of state, w? Most 
inflation models in which the inflationary periods last for 
much longer than 60 e-folds predict flk ~ 10 -5 , which is 
three orders of magnitude below the current constraint 
(see § 13. 4p . In this subsection, we use a "strong inflation 
prior," imposing a flatness prior, and explore dark energy 
models in the context of such inflation models. We shall 
explore curved universes in § 15.31 

Figure QT] shows the constraints on w and the present- 
day dark energy density, f^A. The WMAP data alone 
cannot constrain this parameter space very well, as cer- 
tain combinations of w and can produce very similar 



angular diameter distances out to the decoupling epoch. 

The HST prior helps a little bit -0.47 < 1 + w < 
0.42 (95% CL) by constraining Q A : the WMAP data 
measure Q m h 2 , and a flatness prior imposes a constraint, 
Oa = 1 — (Q, m h 2 )/h 2 ; thus, an additional constraint on 
h from the HST Key Project helps determine f^A better. 

The current angular diameter distance measurements 
from the BAO do not quite break the degeneracy between 
w and as they constrain the distances at relatively 
low redshifts, z = 0.2 and 0.35, whereas the transition 
from matter to dark energy domination, which is sensi- 
tive to w, happens at earlier times. (Therefore, the future 
BAO surveys at higher redshifts should be more sensi- 
tive to w.) The WMAP+BAO yields -0.68 < 1 + w < 
0.21 (95% CL). 42 

Finally, the Type la supernova data break the de- 
generacy nicely, as their constraint on this parameter 
space is nearly orthogonal to what is determined by 
the CMB data: WMAP+SN yields -0.12 < 1 + w < 
0.14 (95% CL). 43 

With a flatness prior, the constraint on w from 
SN is so powerful that WMAP+SN is similar to 
WMAP+BAO+SN. We conclude that, when the equa- 
tion of state does not depend on redshifts, dark energy 
is consistent with vacuum energy, with —0.12 < 1 + w < 
0.13 (95% CL) 44 (from WMAP+BAO+SN), in the con- 
text of a flat universe at the level of curvature that is 
predicted by long-lasting inflation models. 

5.3. Constant equation of state: Curved universe 

In this subsection we do not assume a flat universe, 
but do assume a constant equation of state. (For time- 
dependent equation of state, see § 15.4.21 ) As we dis- 
cussed in § 13. 4[ the WMAP data alone are unable to 
place meaningful constraints on the spatial curvature of 
the universe; however, two or more distance or expan- 
sion rate measurements break the degeneracy between 
flk and fi m . As Fig. [6] shows, the combination of the 
WMAP measurement of the distance to the decoupling 
epoch at z ~ 1090, and the distance measurements out 
to z = 0.2 and 0.35 from BAO, strongly constrains the 
curvature, at the level of 1-2%. 

However, when dark energy is dynamical, we need 
three distance indicators that cover a wide range of red- 
shift. As the current SN data cover a wider range in 
redshifts, 0.02 < z < 1.7, than the BAO data, the SN 
data help to constrain the evolution of dark energy, i.e., 
w. 

Figure [12] shows the constraints on w 
and fifc from the WMAP 5-year data alone, 
WMAP+HST, WMAP+BAO, WMAP+SN, as 
well as WMAP+BAO+SN. The middle panel is 
particularly illuminating. The WMAP+BAO com- 
bination fixes f2fc, nearly independent of w. 45 The 
WMAP+SN combination yields a degeneracy line 



+0.21 



(WMAP+BAO; C fe = 0). 



42 The 68% limit is to = -1.15_ .... 

43 The 68% limit is w = -0.977±°;^ (WMAP+SN; Q k = 0). 

44 The 68% limit is w = -0.992±g°H (WMAP+BAO+SN; 

n k = o). 

45 For the WMAP+BAO limit, there is a long degenerate valley 
with a significant volume at w < — 1. Models anywhere in this 
valley are good fits to both data sets. It is dangerous to marginalize 
over these degenerate parameters as conclusions are very sensitive 
to the choice and the form of priors. 
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Fig. 11. — Constraint on the time- inde pendent (constant) dark energy equation of state, w, and the present-day dark energy density, Ha, 
assuming a flat universe, = (§ 15. 211 . Note that we have imposed a prior on w, V) > —2.5. {Left) Joint two-dimensional marginalized 
distribution of w and S7 fc . The contours show the 68% and 95% CL. The WMAP-only constraint (light blue) is compared with WMAP+HST 
(gray), WMAP+BAO (red), WMAP+SN (dark blue), and WMAP+BAO+SN (purple). This figure shows how powerful a combination 
of the WMAP data and the current SN data is for constraining w. (Middle) One-dimensional marginalized constraint on ui for a flat 
universe from WMAP+HST (gray), WMAP+BAO (red), and WMAP+SN (dark blue). The WMAP+BAO+SN result (not shown) is 
essentially the same as WMAP+SN. (Right) One-dimensional marginalized constraints on Q\ for a flat universe from WMAP+HST (gray), 
WMAP+BAO (red), and WMAP+SN (dark blue). The WMAP+BAO+SN result (not shown) is essentially the same as WMAP+SN. 
See Fig. 1121 for the constraints on w for non-flat universes. Note that neither BAO nor SN alone is able to constrain w. they need the 
WMAP data for lifting the degeneracy. Note also that BAO+SN is unable to lift the degeneracy either, as BAO needs the sound horizon 
size measured by the WMAP data. 
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Fig. 12. — Joint two-dime nsio nal marginalized constraint on the time-independent (constant) dark energy equation of state, w, and the 
curvature parameter, Q k (§ 15.31 1. Note that we have imposed a prior on w, w > —2.5. The contours show the 68% and 95% CL. (Left) 
The WMAP-only constraint (light blue; 95% CL) compared with WMAP+BAO+SN (purple; 68% and 95% CL). This figure shows how 
powerful the extra distance information from BAO and SN is for constraining and w simultaneously. (Middle) A blow-up of the left 
panel, showing WMAP+HST (gray), WMAP+BAO (red), WMAP+SN (dark blue), and WMAP+BAO+SN (purple). This figure shows 
that we need both BAO and SN to constrain fl k and w simultaneously: WMAP+BAO fixes O fe , and WMAP+SN fixes w. (Right) The 
same as the middle panel, but with the BAO prior re- weighted by a weaker BAO prior from the SDSS LR G sample (Eiscnstei n et al.l 
120051) . The BAO data used in the other panels combine the SDSS main and LRG, as well as the 2dFGRS data UPercival et al.H2007l ). The 
constraints from these are similar, and thus our results are not sensitive to the exact form of the BAO data sets. Note that neither BAO 
nor SN alone is able to constrain ui or they need the WMAP data for lifting the degeneracy. Note also that BAO+SN is unable to lift 
the degeneracy either, as BAO needs the sound horizon size measured by the WMAP data. 



that is tilted with respect to the WMAP+BAO line. 
The WMAP+BAO and WMAP+SN lines intersect 
at flk ~ and w ~ — 1, and the combined con- 
straints are -0.0179 < Q k < 0.0081 (95% CL) and 
-0.14 < 1 + w < 0.12 (95% CL). 46 It is remarkable 
that the limit on f2fc is as good as that for a vacuum 
energy model, -0.0178 < n k < 0.0066 (95% CL). This 
is because the BAO and SN yield constraints on and 
w that are complementary to each other, breaking the 
degeneracy effectively. 



46 The 68% 

q+0.067 



limits are fli 



-0.0049±°;°064 and 



-1.006"* 



(WMAP+BAO+SN). 



These limits give the lower bounds to the curvature 
radii of the observable universe as i? C urv > 33 /i _1 Gpc 
and i? C urv > 22 /i _1 Gpc for negatively and positively 
curved universes, respectively. 

Is the apparent "tension" between the WMAP+BAO 
limit and the WMAP+SN limit in Fig. [12] the signa- 
ture of new physics? We have checked this by the BAO 
distance scale ou t to z = 0.35 from the S DSS LRG sam- 
ple, obtained bv lEisenstein et all (|2005l ). instead of the 
z = 0.2 and z — 0.35 constraints based on the combi- 
nation of SPSS LRGs with th e SDSS main sample and 
2dFGRS (|Percival et al J 120071 ). While is it not an inde- 
pendent check, it does provide some measurement of the 
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TABLE 10 

WMAP DISTANCE PRIORS OBTAINED FROM THE 
WMAP 5-YEAR FIT TO MODELS WITH SPATIAL 
CURVATURE AND DARK ENERGY. THE 

CORRELATION COEFFICIENTS ARE: 
Tl A ,R = 0.1109, r u>s . = 0.4215, AND 
r R z = 0.6928. 





5-year ML" 


5-year Mean 6 


Error, a 


Ia(z,) 


302.10 


302.45 


0.86 




1.710 


1.721 


0.019 




1090.04 


1091.13 


0.93 



"Maximum likelihood values (recommended) 
'Mean of the likelihood 

"Maximum likelihood values (recommended) 
d Mean of t he l ikelihood 
e Equation II66I I 

sensitivity of the constraints to the details of the BAO 
data set. 

The right panel of Fig. [T^] shows that the results are 
not sensitive to the exact form of the BAO data sets. 47 
Eisenstein et al.'s BAO prior is a bit weaker than Perci- 
val et al.'s, and thus the WMAP+BAO contours extend 
more to w > —1. The important point is that the di- 
rection of degeneracy does not change. Therefore, the 
combined limits from WMAP, SN and Eisenstein et al.'s 
BAO, -0.15 < 1 + w < 0.13 (95% CL), and -0.0241 < 
n k < 0.0094 (95% CL) are similar to those with Per- 
cival et al.'s BAO, -0.14 < l + w < 0.12 (95% CL), 
and -0.0179 < fl k < 0.0081 (95% CL). As expected, a 
weaker BAO prior resulted in a weaker limit on Q k - 

While the above argument suggests that there is no 
serious tension between WMAP+BAO and WMAP+SN 
constraints, would it be possible that the tension, if any, 
could be caused by the WMAP data? As the BAO 
data use the sound horizon size measured by the WMAP 
data, r s (z £ j), some systematic errors causing the mis- 
calculation of r s (z ( i) could lead to a mis-interpretation 
of the BAO data. The current measurement errors in 
r s (z d )/D v (z) from the BAO data are 2.9% at z = 0.2 
and 3.0% at z = 0.35. On the other hand, WMAP mea- 
sures r s (zd) with 1.3% accuracy (see Table [3]). We are 
confident that the systematic error in r s (zd), if any, is 
smaller than the statistical error; thus, it is unlikely that 
WMAP causes a mis-interpretation of the BAO data. 

From these studies, we are able to place rather strin- 
gent, simultaneous limits on Q k (to 1-2% level, depend- 
ing upon the sign), and w (to 14% level). The spatial 
curvature is consistent with zero, and the dark energy is 
consistent with vacuum energy. How does this conclusion 
change when we allow w to vary? 

5.4. WMAP distance priors for testing dark energy 
models 

5.4.1. Motivation 

Dark energy influences the distance scales as well as 
the growth of structure. The CMB power spectra are 

47 To obtain the WMAP+BAO contours in the right panel 
of Fig. [12J we have re-weighted the WMAP+BAO data 
in the middle panel of Fig. 1121 by the likelihood ratio of 
L(Eisenstein's BAO)/L(Percival's BAO). As a result the contours 
do not extend to w ~ 0; however, the contours would extend more 
to w ~ if we ran a Markov Chain Monte Carlo from the beginning 
with Eisenstein et al.'s BAO. 



TABLE 11 
Inverse covariance matrix for the 
wmap distance priors 





l A (z*) 


«(*.) 


z* 


Ia(z*) 


1.800 


27.968 


-1.103 






5667.577 


-92.263 


z, 






2.923 



TABLE 12 

Inverse covariance matrix for the extended WMAP 
is 100O 6 h 2 = 2.2765. 





l A {z«) 


R(z-,) 


z* 


won b h 2 


Ia(z») 


31.001 


-5015.642 


183.903 


2337.977 


R(z«) 




876807.166 


-32046.750 


-403818.837 


z* 






1175.054 


14812.579 


wm b h 2 








187191.186 



sensitive to both, although sensitivity to the growth of 
structure is fairly limited, as it influences the CMB power 
spectrum via the ISW effect at low multipoles (I < 10), 
whose precise measurement is hampered by a large cos- 
mic variance. 

On the other hand, CMB is sensitive to the distance 
to the decoupling epoch via the locations of peaks and 
troughs of the acoustic oscillations, which can be mea- 
sured precisely. More specifically, CMB measures two 
distance ratios: (i) the angular diameter distance to the 
decoupling epoch divided by the sound horizon size at 
the decoupling epoch, D,4(z*)/r s (z*), and (ii) the angu- 
lar diameter distance to the decoupling epoch divided 
by the Hubble horizon size at the decoupling epoch, 
Da{z^)H(z^)/c. This consideration motivates our us- 
ing these two distance ratios to constrain various dark 
energy models, in the presence of t he spatial curvature, 
on the basis of dist ance information (|Wang fc Mukhcrice 
[2TMIWri^[2Tj0l . 

We shall quantify the first distance ratio, 
DA{z*)/r s (z*), by the "acoustic scale," ^a, defined 

by 



Ia = (1 + 2*)- 



(65) 



where a factor of (1 + z*) arises because Da(z*) is the 
proper (physical) angular diameter distance (Eq. [2]), 
whereas r s (z*) is the comoving sound horizon at z* 
(Eq. [5]). Here, we shall use the fi tting function of z* 
proposed bv lHu &: Sugivamal (|1996l ): 



1048 [1 + 0.00124(fVi 2 r - 738 ] 1 + 9i (tt m h 



2\92 



where 



,9i 



92 



0.0783(» b fe 2 )- - 238 

l + 39.5(fW) - 763 ' 

0.560 
l + 21.1(fi 6 /i 2 ) 1 - 81. 



(66) 

(67) 
(68) 



Note that one could also use the peak of the probability 
of last scattering of photons, i.e., the peak of the visi- 
bility function, to define the decoupling epoch, which we 
denote as Zd ec - Both quantities yield similar values. We 
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Fig. 13. — Joint two-dimensional marginalized constraint on the time-independent (cons tant) d ark energy equation of state, w, and the 
curvature parameter, f^, derived solely from the WMAP distance priors (I a, R, z t ) (see § 15.4.11 . combined with either BAO (red) or SN 
(dark blue) or both (purple). The contours show the Axj^ = 2.30 (68. 3% CL) and Ax^ a] = 6.17 (95.4% CL). The left (BAO data from 



Pcr cival ct al. (2007)) and right (BAO data from Eiscnstcin ct al. (2005)) panels should be compared with the middle and right panels of 
Fig. 1121 respectively, which have been derived from the full WMAP data combined with the same BAO and SN data. While the WMAP 
distance priors capture most of the information in this parameter space, the constraint is somewhat weaker than that from the full analysis. 



shall adopt z* here because it is easier to compute, and 
therefore it allows one to implement the WMAP distance 
priors in a straightforward manner. 

The second distance ratio, Da(z*)H(z*) /c, is often 
called the "shift parameter," R, given by (jBond et al.l 
1997) 



i?(z* 



(1 + z*)D A {z*). 



(69) 



This quantity is different from Da(z*)H(z*)/c by a fac- 
tor of \/l + 2*, and also ignores the contributions from 
radiation, curvature, or dark energy to H(z*). Never- 
theless, we shall use R to follow the convention in the 
literature. 

We give the 5-year WMAP constraints on I a, 
R, and z* that we recommend as the WMAP 
distance priors for constraining dark energy mod- 
els. However, w e note an important c aveat . As 
pointed out by [Ekareiv fc Multamakl (|2007l ) and 
ICorasaniti fe Melchiorril (|2008f ). the derivation of the 
WMAP distance priors requires us to assume the un- 
derlying cosmology first, as all of these quantities are 
derived parameters from fitting the CMB power spectra. 
Therefore, one must be careful about which model one 
is testing. Here, we give the WMAP distance priors, 
assuming the following model: 

• The standard Friedmann-Lemaitre-Robertson- 
Walker universe with matter, radiation, dark 
energy, as well as spatial curvature. 

• Neutrinos with the effective number of neutri- 
nos equal to 3.04, and the minimal mass (m„ ~ 
0.05 eV). 

• Nearly power-law primordial power spectrum of 
curvature perturbations, \dn s /dhik\ <C 0.01. 

• Negligible primordial gravitational waves relative 
to the curvature perturbations, r <C 0.1. 



• Negligible entropy fluctuations relative to the cur- 
vature perturbations, a <C 0.1. 

In Fig. [T3] we show the constraints on w and fifc from 
the WMAP distance priors (combined with BAO and 
SN). We find a good agreement with the full MCMC re- 
sults. (Compare the middle and right panels of Fig. [T2] 
with the left and right panels of Fig. [T3] respectively.) 
The constraints from the WMAP distance priors are 
slightly weaker than the full MCMC, as the distance pri- 
ors use only a part of the information contained in the 
WMAP data. 

Of course, the agreement between Fig. [12] and [13] does 
not guarantee that these priors yield good results for the 
other, more complex dark energy models with a time- 
dependent w; however, the previous studies indicate that, 
under the assumptions given above, these priors can be 
used to constrain a wide variety of dark energy models 
(IWang fc M ukhcri ei l2007t lElgarOv fc MultamaB 1200 ' 
ICorasaniti fc Melchiorrill2008f ). See also iLi et all ([2008 
for the latest comparison between the WMAP distance 
priors and the full analysis. 

Here is the prescription for using the WMAP distance 
priors. 



compute z» from 



(1) For a given f2&/i 2 and £l m h 2 , 
Eq. §<o$. 

(2) For a given Hq, Q m h 2 , fl r h 2 (which includes iV e gt = 
3.04), and u>(z), compute the expansion rate, 
H(z), from Eq. ([7]), as well as the comoving sound 
horizon size at z„, r s (z«), from Eq. ([6]). 

(3) For a given f2fc and H{z) from the previous step, 
compute the proper angular diameter distance, 
D A (z), from Eq. ©. 



(4) Use Eq. §5$ and 
respectively. 



to compute Ia{z„) and i?(z*), 
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(5) Form a vector containing Xi 
order. 



(I A, R, 2*) in this 



(6) Use Table [TU] for the data vector, 



(1 



WMAP jjWMAP „WMAP 



). We recommend the 



maximum likelihood (ML) values. 

(7) Use Table [TT] for the inverse covariance matrix, 

(8) Compute the likelihood, L, as Xwmap = ~2 m -k = 
(xi - d l )(C~ 1 ) ij (x j - dj). 

(9) Add this to the favorite combination of the cosmo- 
logical data sets. In this paper we add Xwmap 



to the BAO and SN data, i.e., Xtotai 



Xbao 



Xsn- 



Xwmap 



(10) Marginalize the posterior distribution over f^ft 2 , 
Q, m h 2 , and Hq with uniform priors. Since the 
WMAP distance priors combined with the BAO 
and SN data provide tight constraints on these pa- 
rameters, the posterior distribution of these param- 
eters is close to a Gaussian distribution. There- 
fore, the marginalization is equivalent to minimiz- 
ing y 2 „+,i with r e spect to flhh 2 , f2 m /i 2 , and H (see 



also ICashl H97l IWrightl l200l . We use a down- 
hill simplex method for mi nimization (amoe ba rou- 
tine in Numerical Recipes (| Press et a l. 1992)). The 
marginalization over flbh 2 , f2 m /i 2 , and Hq leaves us 
with the the marginalized posterior distribution of 
the dark energy function, w(a), and the curvature 
parameter, fife. 

Note that this prescription eliminates the need for run- 
ning the Markov Chain Monte Carlo entirely, and thus 
the computational cost for evaluating the posterior dis- 
tribution of w(a) and fife is not demanding. In § !5.4.2l we 
shall apply the WMAP distance priors to constrain the 
dark energy equation of state that depends on redshifts, 
w = w(z). 

For those who wish to include an additional 
prior on fl^h 2 , we give the inverse covariance ma- 
trix for the "extended" WMAP distance priors: 
(Ia{z*), R(z*), z*, lOOrib/i 2 ) , as well as the maximum 
likelihood value of 100fi b /i 2 , in Table HU We note, 
however, that it is sufficient to use the reduced set, 
(Ia(z*), R(z*),z*), as the extended WMAP distance pri- 
ors give ve ry similar constraints on dark energy (see 
IWanal200l . 

5.4.2. Application of the WMAP distance priors: 
Variable equation of state 

In this subsection we explore a time-dependent equa- 
tion of state of dark energy, w(z). We use the following 
parametrized form: 



w(a) 



aw (a) 



Cltr 



(70) 



where w(a) = Wq + (1 — a)w a . We give a motivation, 
derivation, and detailed discussion on this form of w(a) 
in Appendix [Cj This form has a number of desirable 
properties: 



• w(a) approaches to —1 at early times, a < atrans, 
where dtrans = 1/(1 + ztrans) is the "transition 
epoch," and Ztrans is the transition redshift. There- 
fore, the dark energy density tends to a constant 
value at a < a tran s. 

• The dark energy density remains totally sub- 
dominant relative to the matter density at the de- 
coupling epoch. 

• We recover the widely used linear form, w(a) = 
wo + (1 — a)w a (jChevallier fc Polarskl 2001; Linde3 
2003), at late times, a > ^trans- 

• The early-time behaviour is consistent with some of 
scalar field models classif i ed as the "thawing mod- 
els," ( Caldwe ll fc Lind"erl |2005| ) in which a scalar 
field was moving very slowly at early times, and 
then began to move faster at recent times. 

• Since the late-time form of w(a) allows w(a) to go 
below — 1, our form is more general than models 
based upon a single scalar field. 

• The form is simple enough to give a closed, an- 
alytical form of the effective equation of state, 

Weft (a) = (lna)^ 1 J lna din a' w(a') (Eq. lC6j) . which 
determines the evolution of the dark energy density, 
p dc (a) = / 9 dc (0)ar 3[1+tOrff(a)I ; thus, it allows one to 
compute the evolution of the expansion rate and 
cosmological distances easily. 

While this form contains three free parameters, wq, w a 
and Ztrans, we shall give constraints on the present-day 
value of w, wq = w(a — 1), and the first derivative of 



w at present, w 



dw/dz 



instead of wq and 



as they can be compared to the previous results in the 
literature more directly. We find that the results are not 
sensitive to the exact values of Ztrans- 

In Fig. [14] we present the constraint on wo and w' that 
we have derived from the WMAP distance priors (I a, 
R, z*), combined with the BAO and SN data. Note 
that we have assumed a flat universe in this analysis, 
althoug h it is straightforward to inc lude the spat i al cur - 
vature. IWang fc Mukheried (|2007f ) and IWrightl (|2007l ) 
show that the two-dimensional distribution extends more 
towards south-east, i.e., w > — 1 and w' < 0, when the 
spatial curvature is allowed. 

The 95% limit on wo for z trans = 10 is 
0.21 48 . Our results are consistent 



work using the WMAP 3-y e ar data ( see [Z 
Wang fc Mukheried [20071 : iWrightJ [2007 



tng 

20071 for recent work and references therein). The 



-0.33 < 1+wo < 
wit h the previous 
lao et al.ll200 
Lazkoz et al 



WMAP 5-year data help tighten the upper limit on w' 
and the lower limit on Wq, whereas the lower limit on w' 
and the upper limit on wo come mainly from the Type 
la supernova data. As a result, the lower limit on w' 
and the upper limit on wo are sensitive to whether we 
include the systematic errors in the supernova data. For 
this investigation, see Appendix [D] 

Alternatively, one may take the linear form, w(a) = 
wq + (1 — a)w a , literally and extend it to an arbitrarily 



48 The 68% intervals are wo 
(WMAP+BAO+SN; O fc = 0). 



-1.06±0.14 and w' = 0.36±0.62 
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Fig. 14. — Constraint on models of time-dependent dark energy equation of state, w(z) (Eq. [70]), derived from the WMAP distance priors 
(/A; R, and z») combined with the BAO and SN distance data 15. 4. 21 1. There are three parameters: wo is the value of w at the present 
epoch, wq = w(z = 0), w' is the first derivative of w with respect to z at z = 0, w' = dw/dz\ z=f); and ztrans is the transition redshift above 
which w(z) approaches to —1. Here, we assume flatness of the universe, Q k = 0. (Left) Joint two-dimensional marginalized distribution of 
wo and w' for Ztrans = 10. The constraints are similar for the other values of ztrans- The contours show the Axt ota j = 2.30 (68.3% CL) 
and Ax? otal = 6.17 (95.4% CL). [Middle) One-dimensional marginalized distribution of wq for ztrans = 0.5 (dotted), 2 (dashed), and 10 
(solid). (Middle) One-dimensional marginalized distribution of w 1 for ztrans = 0.5 (dotted), 2 (dashed), and 10 (solid). The constraints are 
similar for all ztrans- We do not find evidence for the evolution of dark energy. Note that neither BAO nor SN alone is able to constrain 
wq or w': they need the WMAP data for lifting the degeneracy. Note also that BAO+SN is unable to lift the degeneracy either, as BAO 
needs the sound horizon size measured by the WMAP data. 




Fig. 15. — Constraint on the linear evolution model of dark energy equation of state, w(z) = wo + ui'z/(l + z), derived from the 
WMAP di stan ce priors (I a, R, and z„) combined with the BAO and SN distance data as well as the Big Bang Nucleosynthesis (BBN) 
prior (Eq. |71|). Here, we assume flatness of the universe, f2j- = 0. (Left) Joint two-dimensional marginalized distribution of wo and w' . 
The contours show the Axt otal = 2.30 (68.3% CL) and Axt ota j = 6.17 (95.4% CL). (Middle) One-dimensional marginalized distribution of 
wq. (M iddle) One-dimensional marginalized distribution of w' . We do not find evidence for the evolution of dark energy. Note that ILinderl 
( 2003) defines w' as the derivative of w at z = 1, whereas we define it as the derivative at z = 0. They are related by w' lindcr = 0.5w' WMAp . 

high redshift. This can result in an undesirable situation 
in which the dark energy is as important as the radia- 
tion density at the epoch of the Big Bang Nucleosynthe- 
sis (BBN); however, one can constrain such a scenario 
severely using th e limit on th e expans i on rat e from BBN 
(|Steigmad[2007l ). We follow IWrightl (|2007f ) to adopt a 

Gaussian prior on Q fe = 0and«; = -l 2.211 ± 0.083 

Q fe ^0andw = -1 2.212 ±0.084 

!2 t = and to ^ -1 2.208 ± 0.087 

n fe ^0and«j^-l 2.210 ± 0.084 

a k = 0, w = -1 and m v > 2.212 ± 0.083 

U k = 0, w ^ -1 and m v > 2.218 ± 0.085 



TABLE 13 

Amplitude of curvature perturbations, TZ, 

MEASURED BY WMAP AT k W MAP = 0.02 Mpc" 1 



Model 



10 9 X A^(k W MAp) 



A (1 + Z B BN) 3[1+t0cff(2:BBN)1 



On(l + zbbn) 3 + O r (l + zbbn) 4 + CMl + -zbbn) 2 
= 0.942 ±0.030, (71) 

where we have kept Vl m and for definiteness, but they 
arc entirely negligible compared to the radiation density 
at the redshift of BBN, zbbn = 10 9 . Figure [TBI shows the 
constraint on wo and w' for the linear evolution model 
derived from the WMAP distance priors, the BAO and 
SN data, and the BBN prior. The 95% limit on wq is 



WMAP Normalization Prior 



2.21 ± 0.09 



—0.29 < 1 + wo < 0.21 49 , which is similar to what we 
have obtained above. 



49 The 68% intervals are wo = 
(WMAP+BAO+SN+BBN; Sl k 



-1.04±0.13 and w' ■■ 
0). 



: 0.24±0.55 
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5.5. WMAP normalization prior 

So far, we have been mainly using the distance 
information to constrain the properties of dark energy; 
however, this is not the only information that one 
can use to constrain the properties of dark energy. 
The amplitude of fluctuations is a powerful tool for 
distinguishing between d a rk ene r gy and modifications 
to gra vity (llshak et al.l [20061: iKovama fc Maarteni 



Linder fc Cahnl l2007t lUpadhvel 120071: I Zhang et al l 



2007t iDaniel et al 



20061 lAmarzguioui et"ail l2006t iDore et all 120071: 



2007t I Yamamoto et alj 120071: iChiba fc T akahashi 12007 



Bcrtschingcr & Zukin 



Beanet alJ 120071: |Hu fc Sawickil 120071: ISongetall 



Jain fc Zhand 120071: 
IAminetal.1 l2007M Hul 

2008), as well as for determining the mass of neutrinos 
(|Hu et al.|[l998HLesgourgues fc Pastoril2006h . 

The microwave background observations measure the 
amplitude of fluctuations at the decoupling epoch. By 
combining this measurement with the amplitude mea- 
sured from various low redshift tracers, one can learn 
more about the dark energy properties and the mass of 
neutrinos. 

The overall amplitude of CMB anisotropy is set by the 
amplitude of primordial curvature perturbations, 1Z. For 
example, on very large angular scales where the Sachs- 
Wolfe limit can be used, the temperature anisotropy is 
given by AT/T = —TZ/5 or, in terms of the curvature 
perturbation during the matter era, it is given by 
AT/T = — $/3. On small angular scales where the 
acoustic physics must be taken into account, we have 
the acoustic oscillation whose amplitude is also given by 
TZ. 

This motivates our reporting the "WMAP normal- 
ization," a measurement of the overall normalization 
of the curvature perturbations expressed in terms of 
A^ikwMAp), where A^(fc) = fc 3 P 7? ,(fc)/(27r 2 ) is a con- 
tribution to the total variance of 1Z, (1Z 2 ) , per logarithmic 
interval of fc (see also Eq. [15] and descriptions below it). 

Here, kwMAP is different from fco = 0.002 Mpc -1 that 
we used to define n s , dn s /dhik, r, and A^(fco) reported 
in Table Q] or 2J The goal in this subsection is to give the 
normalization that is as model independent as possible. 

At fco = 0.002 Mpc -1 , for example, we find 
10 9 A|(fc ) = 2.48, 2.41, and 2.46 for a flat ACDM 
model, a curved ACDM model, and a flat ACDM model 
with massive neutrinos. The scatter between these val- 
ues comes solely from the fact that fco = 0.002 Mpc -1 is 
not a right place to define the normalization. In other 
words, this is not the pivot scale of the WMAP data. 

We find that kwMAP = 0.02 Mpc -1 , i.e., a factor of 10 
larger than fco, gives similar values of A^kwMAp) for a 
wide range of models, as summarized in Table [T3l From 
these results, we give the WMAP normalization prior: 

A 2 n {k WM Ap) = (2.21 ± 0.09) x 10" 9 , (72) 
which is valid for models with Qj~ ^ 0, w ^ — 1, or m„ > 
0. However, we find that these normalizations cannot be 
used for the models that have the tensor modes, r > 0, 
or the running index, dn s /d\nk ^ 0. We failed to find a 
universal normalization for these models. Nevertheless, 
our WMAP normalization given by Eq. (f72|) is still valid 
for a wide range of cosmological models. 

How can one use the WMAP normalization? In or- 
der to predict the linear matter density power spectrum, 



P\in(k), one needs to relate the primordial curvature per- 
turbations, IZk, to the linear matter density fluctuations 
at arbitrary redshifts, 8 m ,v.{z). From Ei nstein's equa- 
tions , we find (see e.g., Appendix C of iTakada et al.1 
2006) 

2fc 3 

«m,k(z) = -^-TZ k T(k)D(k, z), (73) 

where D(k, z) and T(k) is the linear growth rate and the 
matter transfer function normalized such that T{k) — > 1 
as fc — * 0, and (1 + z)D(k, z) — ► 1 as fc — * during the 
matter era (e.g., z = 30, where the radiation density is 
less than 1% of the matter density), respectively. Note 
that -D(fc, z) does not depend on fc when neutrinos are 
mass less; however, it depend s on fc when they are massive 
(e.g. JHu fc Eisenstein| [l998). The linear matter density 
power spectrum is given by 

fc 3 P lin (fc,z) _, 001 ^ nnnx r in _ 9 / 2fc 2 J 



2tt 2 



(2.21 ±0.09) x 10" 



xD 2 {k,z)T 2 (k) 



k]VMAP 



(74) 



One application of the WMAP normalization is the 
computation of the present-day normalization of matter 
fluctuations, which is commonly expressed in terms of 
OS, given by 



al = (2.21 ± 0.09) x 10 



/ 2 



f D 2 (k, 



fc 4 



^ = 0)T 2 (fc)-^ 
3sin(fci?) 3cos(fci?) 



kwMAP 



(fci?) 3 

where R = 8 /i _1 Mpc. 



(fci?) 2 



(75) 



Both the dark energy properties 
(or modified gravity) and the mass of neutrinos change 
the value of D(k,z — 0). The transfer function, T(k), 
is much less affected, as long as neutrinos were still rel- 
ativistic at the decoupling epoch, and the dark energy 
or modified gravity effect was unimportant at the decou- 
pling epoch. 

Ignoring the mass of neutrinos and modifications to 
gravity, one can obtain the g rowth rate by solving the 
following differential eq uation ([Wang fc Steinhard3ll998l : 
iLinder fc Jenkinsll2005 



9 



dim 



(Slfc(a) - 3w cS (a)n de (a)) 



dg 

din a 



2tt k (a) + -(1 - w e ff(a))tt d e{a) 



5 (a) = 0, 



where 



9{a) : - 
fife (a) = 

^de(a) ; 
Weff(a); 



D(a) 



(l + z)D(z), 



a 2 H 2 {aY 
VaH 2 

a 3[l+w Q{f (a)] H 2( a y 
An a 



l 

In a 



din a' w(a') 



(76) 

(77) 
(78) 

(79) 

(80) 
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Fig. 16. — Predictions for the present-day amplitude of matter fluctuations, erg, as a function of the (constant) dark energy equation of 
state parameter, w, derived from the full WMAP data (blue) as well as from WMAP+BAO+SN (red). The contours show the 68% and 
95% CL. (Left) Flat universe, f2j. = 0. (Right) Curved universe, ^ 0. 



During the matter era, g(a) does not depend on a; 
thus, the natural choice for the initial conditions are 
ff(ainitiai) = 1 and dg/d\na\ a=a _ tM = 0, where a in itiai 
must be taken during the matter era, e.g., ainitiai = 1/31 
(i.e., z = 30). 

In Fig. [IB] we show the predicted values of as as a 
function of w in a flat universe (left panel) and curved 
universes (right panel). (See the middle panel of Fig. [17] 
for eg as a function of the mass of neutrinos.) Here, 
we have used the full information in the WMAP data. 
The normalization information alone is unable to give 
meaningful predictions for erg, which depends not only 
on A^(kwMAp), but also on the other cosmological pa- 
rameters via T{k) and D(k,z = 0), especially w and 
fl m h 2 . While the predictions from the WMAP data 
alone are still weak, adding the extra distance informa- 
tion from the BAO and SN data helps improve the pre- 
dictions. We find er 8 = 0.807T. ; 044 for a flat universe, 
and as — 0.795 ± 0.046 for curved universes. 

By combining these results with as measured from var- 
ious low redshift tracers, one can reduce the remain- 
ing correlation between w and as to obtain a better 
limit on w. The precision of the current data from 
weak lensing surveys is comparable to these predictions, 
e.g., cT 8 (fi m /0.25) - 64 = 0.785 ± 0.043 (Canad a-France- 
Hawa ii Telescope Legacy Survey (CFHTLS); iFu et al.l 
2008). The weak lensing surveys will soon become pow- 
erful enough to yield smaller uncertainties in as than 
predicted from WMAP+BAO+SN. 

6. NEUTRINOS 

In this section, we shall use the WMAP data, com- 
bined with the distance information from BAO and SN 



observations, to place limits on the total mass of massive 
neutrino species (§ 16. ip . as well as on the effective num- 
ber of neutrino-like species that were still relativistic at 
the decoupling epoch (§ 



6.1. Neutrino Mass 

6.1.1. Motivation 

The existence of non-zero neutrino masses has been 
established firmly b y the experiments detecting at- 
mospheric n eutrinos jHirata et al.l Il992t iFukuda et al.l 
[1991 fl998t lAllison et all 119991 : lAmbrosio et all 120017 
solar neutrinos [ Davis et al.l Il968t ICleveland et all 



19981: lHampel et al I 119991: lAbdurashitov et all 119991: 



Fukuda et al. 20033;lAhSad et al. 2002: I Ahmed et al.l 
20041) . reactor neutrinos (jEguchi et a l. 2003: Araki et al.l 
2005tl and accelerat or beam neutrinos (|Ahn et al.l 120031 : 
Michael et all |2006[ ). These experiments have placed 
stringent limits on the squared mass differences between 
the neutrino mass eigenstates: Am|i — 8 x 10~ 5 eV 2 
from the solar and reactor experiments, and Am| 2 — 
3 x 10 -3 eV 2 from the atmospheric and accelerator beam 
experiments. 

One needs different experiments to measure the ab- 
solute masses. The next-generation tritium /3-decay ex- 
periment, KATRIN 50 , is expected to reach the electron 
neutrino mass of as small as ~ 0.2 eV. Cosmology has 
also been providing useful limits on the mass of neutrinos 
(seelDolgovll2002tlElgarov fc Lahavll2005HTegmarkll2005l; 
Lesgourgues k Pastorl 120061 : iFukugital 120061 : iHannestadl 
2006bl for reviews). Since the determination of the neu- 
trino mass is of fundamental importance in physics, there 

50 |http://www-ik.fzk.de7 Tr katrin 
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is enough motivation to pursue cosmological constraints 
on the neutrino mass. 

How well can CMB constrain the mass of neutri- 
nos? We do not expect massive neutrinos to affect 
the CMB power spectra very much (except through the 
gravitational lensing effect), if they were still relativis- 
tic at the decoupling epoch, z ~ 1090. This means 
that, for massive neutrinos to affect the CMB power 
spectra, at least one of the neutrino masses must be 
greater than the mean energy of relativistic neutrinos 
per particle at z ~ 1090 when the photon temperature 



Goobar et ail 120061 : iFeng et alj[2006t iLesgourgues et all 
20071) . These analyses reached different limits depending 



of the universe was T~ 



3000 K 



0.26 eV. Since 



the mean energy of relativistic neutrinos is given by 
(E) = (77r 4 T„)/(180C(3)) ~ 3.157; = 3.15(4/11) 1 / 3 T 7 , 
we need at least one neutrino species whose mass sat- 
isfies m v > 3.15(4/11) 1 / 3 T 7 ~ 0.58 eV; thus, it would 
not be possible to constrain the neutrino mass using the 
CMB data alone, if the mass of the heaviest neutrino 
species is below this value. 

If the neutrino mass eigenstates are degenerate with 
the effective number of species equal to 3.04, this argu- 
ment suggests that J2 m v ~ 1-8 eV wou ld be the limit 
to wh ich the CMB data are sensitive. Ichikawa ct al. 
(2005) argue that ^ m v ^1.5 eV would be the limit for 
the CMB data alone, which is fairly close to the value 
given above. 

In order to go beyond ~ 1.5 eV, therefore, one needs 
to combine the CMB data with the other cosmological 
probes. We shall combine the WMAP data with the dis- 
tance information from BAO and SN to place a limit on 
the neutrino mass. We shall not use the galaxy power 
spectrum in this paper, and therefore our limit on the 
neutrino mass is free from the uncertainty in the galaxy 
bias. We discuss this in more detail in the Analysis sec- 
tion below. 

6.1.2. Analysis 

We assume that, for defmiteness, the neutrino mass 
eigenstates are degenerate, which means that all of the 
three neutrino species have equal masses. 51 We measure 
the neutrino mass density parameter, fl v h 2 , and convert 
it to the total mass, ^2 m u, via 



m„ = N v m v = 94 eV(Q I/ /i 2 ), 



(81) 



where N v is the number of massive neutrino species. We 
take it to be 3.04. Note that in this case the mass density 
parameter is the sum of baryons, CDM, and neutrinos: 

Since the release of the 1- year ([Spergel et al.ll2003f ) and 
3-year ([Spergel et alj [20071 ) results on the cosmological 
analysis of the WMAP data, there have been a num- 
ber of studie s with regard to the limits on the mass 
of neutrinos dHannestadl [20031: lElgar0v fc Lahavi 120031: 
Allen et al.l l2003t iTegmark et all l2004at iBarger et al l 
2004]; [Hannestad fc Raffeltl 120041: TCrottv et all 120041: 



Seliaket all l2005allbl : llchikawa et alJ 120051: IHannestadl 



20051: iLattanzi et all l2005t lHannestad fc Raffeltl 120061 : 



51 While the current cosmological data are not yet sensitive 
to the mass of individual neutrino species, i.e., the mass hi- 
erarchy, this situation may change in the future, with high-2 
galaxy redshift surveys or weak lensing surveys jTakada et al.l 
2006; Slosar 2006; Hannestad & Wong 2007; Kitching ct al. 2008; 



Abda la & Rawlings 2007). 



upon (1) the choice of data sets and (2) the parameters 
in the cosmological model. 

The strongest limits quoted on neutrino masses come 
from combining CMB measurements with measurements 
of the amplitude of density fluctuations in the recent 
universe. Clustering of galaxies and Lya forest observa- 
tions have been used to obtain some of the strongest lim- 
its on neutrino masses (|Seliak et alJl2005ri ISeTiak et all 
-2006: IViel et~aT]|2006[ ). As the neutrino mass increases, 
the am plitude of mass f l uctuations on small scales de- 
crcascs ([Bond et alf l980: Bo nd fc Szalavlll983l : lMal[l996l 
see also the middle panel of Fig. IT?)) , which can be 
used to weigh neutrinos in the universe (|Hu et aUll998l 
ILesgourgues fc Pastorl[2006l ) . 

These analyses are sensitive to correctly calculating 
the relationship between the level of observed fluctua- 
tions in galaxies (or gas) and the mass fluctuation with 
the strongest limits coming from the smallest scales in 
the analyses. With the new WMAP data, these limits 
are potentially even stronger. There are, however, sev- 
eral potential concerns with this limit: there is already 
"tension" between the high level of fluctuations seen in 
the Lyman alpha forest and th e amplitude of mass fluc - 
tuations inferred from WMAP (jLesgourgues et aLll2007h : 
the relationship between gas temperature and density ap- 
pears to be more complic ated than assumed in the pre- 
vious Lya forest analysis (|Kim et alJl2007t iBolton et al.l 
2008); and additional astrophysics could potentially be 
the source of some of the small scale fluctuations seen in 
the Lya forest data. Given the power of the Lya forest 
data, it is important to address these issues; however, 
they are beyond the scope of this paper. 

In this paper, we take the more conservative approach 
and use only the WMAP data and the distance measures 
to place limits on the n eutrino masses . Our approach is 
more conservative than iGoobar et al.l (|2006l ) , who have 
found a limit of 0.62 eV on the sum of the neutrino mass 
from the WMAP 3-year data, the SDSS-LRG BAO mea- 
surement, the SNLS supernova data, and the shape of 
the galaxy power spectra from the SDSS main sample 
and 2dFGRS. While we use the WMAP 5-year data, the 
SDSS+2dFGRS BAO measurements, and the union su- 
pernova data, we do not use the shape of the galaxy 
power spectra. See § 12.31 in this paper or iDunklev et al.1 
(poM) for more detail on this choice. 

In summary, we do not use the amplitude or shape of 
the matter power spectrum, but rely exclusively on the 
CMB data and the distance measurements. As a result 
our limits are weaker than the strongest limits in the 
literature. 

Next, let us discuss (2), the choice of parame- 
ters. A few correlations between the neutrino mass 
and other cosmol ogical parameters have been identified. 
IHannestadl (|2005l ) has found that the limit on the neu- 
trino mass degrades significantly when the dark energy 
eq uation of state, w, i s allowed to vary (see also Fig. 18 
of ISpergel et al.ll2007l ). This correlation would arise only 
when the amplitude of the galaxy or Lya forest power 
spectrum was included, as the dark energy equation of 
state influences the growth rate of the structure forma- 
tion. Since we do not include them, our limit on the 
neutrino mass is not degenerate with w. We shall come 
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Fig. 17. — Constraint on the total mass of neutrinos, ^Z" 1 ^ (§ I6.1.3J I. In all panels we show the WMAP-only results in blue and 
WMAP+BAO+SN in red. (Left) Joint two-dimensional marginalized distribution of Ho and tn v (68% and 95% CL). The additional 
distance information from BAO helps reduce the correlation between Hq and "}2 tn». (Middle) The WMAP data, combined with the 
distances from BAO and SN, predict the present-day amplitude of matter fluctuations, a$, as a function of ^2 m v An independent 
determination of ag will help determine ^,m v tremendously. (Right) Joint two-dimensional marginalized distribution of w and ^]m„. No 
significant correlation is observed. Note that we have a prior on w, w > —2.5, and thus the WMAP-only lower limit on w in this panel 
cannot be trusted. 



back to this point later in the Results section. Inciden- 
tally, our limit is not degenerate with the running index, 
dn s /dlnk, or the tensor-to-scalar ratio, r. 

6.1.3. Results 

Figure fTTl summarizcs our limits on the sum of neutrino 
masses, m v- 

With the WMAP data alone we find J2 m i> < 
1.3 eV (95% CL) for the ACDM model in which w = -1, 
and J2 m v < 1-5 eV (95% CL) for w ^ -1. (We assume 
a flat universe in both cases.) These constraints are very 
similar, which means that w and ^ m v are not degener- 
ate. We show this more explicitly on the right panel of 
Fig. [13 

When the BAO and SN data are added, our lim- 
its improve significantly, by a factor 2, to ^ m v < 
0.67 eV (95% CL) for w = -1, and £m„ < 
0.80 eV (95% CL) for w ^ -1. Again, we do not ob- 
serve much correlation between w and ^2 m v While the 
distances out to either BAO or SN cannot reduce corre- 
lation between f2 m (or Hq) and w, a combination of the 
two can reduce this correlation effectively, leaving little 
correlation left on the right panel of Fig. [T7] 

What information do BAO and SN add to improve the 
limit on ^m„? It's the Hubble constant, Hq, as shown 
in the left panel of Fig. [T7l This effect has been explained 
by llchikawa et al.l (|2005l ) as follows. 

The massive neutrinos modify the CMB power spec- 
trum by their changing the matter-to-radiation ratio at 
the decoupling epoch. If the sum of degenerate neutrino 
masses is below 1.8 eV, the neutrinos were still relativis- 
tic at the decoupling epoch. However, they are definitely 
non-relativistic at the present epoch, as the neutrino os- 
cillation experiments have shown that at least one neu- 
trino species is heavier than 0.05 eV. This means that the 
fi m that we measure must be the sum of fib, fi c , and Q u ; 
however, at the decoupling epoch, neutrinos were still rel- 
ativistic, and thus the matter density at the decoupling- 
epoch was actually smaller than a naive extrapolation 
from the present value. 

As the matter-to-radiation ratio was smaller than one 
would naively expect, it would accelerate the decay 
of gravitational potential around the decoupling epoch. 



This leads to an enhancement in the so-called early in- 
tegrated Sachs- Wolfe (ISW) effect. The larger J] 
is, the larger early ISW becomes, as long as the neutri- 
nos were still relativistic at the decoupling epoch, i.e., 
Y,m„<1.8eV. 

The large ISW causes the first peak position to shift to 
lower multipoles by adding power at I ~ 200; however, 
this shift can be absorbed by a reduction in the value 
of Hp. 52 This is why V m v and Hq are anti-correlated 
(see llchikawa et al. 2005, for a further discussion on this 
subject). 

It is the BAO distance that provides a better limit on 
Hq, as BAO is an absolute distance indicator. The SN 
is totally insensitive to Hq, as their absolute magnitudes 
have been marginalized over (SN is a relative distance 
indicator); however, the SN data do help reduce the cor- 
relation between w and Hq when w is allowed to vary. As 
a result, we have equally tight limits on ^ m„ regardless 
of w. 

Our limit, ^m„< 0.67 eV (95% CL) (for w = -1), is 
weaker than the best limit quoted in the literature, as we 
have not used the information on the amplitude of fluc- 
tuations traced by the large-scale structure. The middle 
panel of Fig. [T7] shows how the WMAP data combined 
with BAO and SN predict the present-day amplitude of 
matter fluctuations, a^, as a function of ^2 m v From this 
it is clear that an accurate, independent measurement of 
erg will reduce the correlation between eg and J2 m v , and 
provide a significant improvement in the limit on ^2 rn v . 

Improving upon our understanding of non-linear as- 



trophy sical effects such as those raised by Bolton et al 
(2008) for the Lya forest data and iSanchez fc Cole 



(|2008ft for the SDSS and 2dFGRS data is a promising- 
way to improve upon the numerical value of the limit, 
as well as the robustness of the limit, on the mass of 
neutrinos. 

52 This is similar to what happens to the curvature constraint 
from the CMB data alone. A positive curvature model, < 0, 
shifts the acoustic peaks to lower multipoles; however, this shift 
can be absorbed by a reduction in the value of Ho- As a result, a 
closed universe with ~ —0.3 and Q/± ~ is still a good fit, if 
Hubb le's constant is as low as Ho ~ 30 km/s/Mpc (Spergcl et al. 
[20071) . 
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6.2. Effective number of neutrino species 
6.2.1. Motivation 

While the absolute mass of neutrinos is unknown, the 
number of neutrino species is well known: it's three. The 
high precision measurement of the decay width of Z into 
neutrinos (the total decay width minus the decay width 
to quarks and charged leptons) , carried out by LEP using 
the production of Z in e + e~ co llisions, has yielded N v = 
2.984 ± 0.008 (|Yao et alJl2006l ). However, are there any 
other particles that we do not know yet, and that were 
relativistic at the photon decoupling epoch? 

Such extra relativistic particle species can change the 
expansion rate of the universe during the radiation era. 
As a result, they change the predictions from the Big 
Bang Nucleosynthesis (BBN) for the ab undance of light 
eleme nts such as helium and deuterium tStcigma n et al.l 
fl97l . One can use this property to place a tight bound 
on the relativistic degrees of freedom, expressed in terms 
of the "effective number of neutrino species," N e g (see 
Eq. [53] below for the precise definition). As the BBN 
occurred at the energy of ~ 0.1 MeV, which is later than 
the neutrino decoupling epoch immediately followed by 
e + e~ annihilation, the value of N e g for 3 neutrino species 
is slightly larger than 3. With other subtle corrections 
included, the current standard value is jVft andard = 3,04 
(iDicus et al.lll982t iGnedin fc Gnedh3ll998llDoTgov et all 
119991 iMangano et alj I2002D . The~2cr interval for iV eff 
from the observed helium a bundance, Yp = 0.240±0.006, 
is 1.61 < N cS < 3.30 (see ISteignranll2007l for a recent 
summary) . 

Many people have been trying to find evidence for the 
extra relativistic degrees of freedom in the universe, using 
the cosmo logical probes such as CMB and the large-scale 
struct ure (lPierpaolill2003l: lHannestadll2003t ICrottv et al " 



12004 iBareer et al.1 12001 
It n+^-r,^,,; r,i I loonst 



200: 

2005: iLattanzi et al 



Trotta k M elchiorri 
Cirelli k Strumia 2006; 



200.1 . _ 

Hannestadl l2006at llchikawa et alJ l2007t IM angan o et alj 
20071: lHamann et alj l2007t Ide Bernardis et all 12008) . 



There is a strong motivation to seek the answer for the 
following question, "can we detect the cosmic neutrino 
background, and confirm that the signal is consistent 
with the expected number of neutrino species that we 
know?" Although we cannot detect the cosmic neutrino 
background directly yet, there is a possibility that we 
can detect it indirectly by looking for the signatures of 
neutrinos in the CMB power spectrum. 

In this section we shall revisit this classical problem 
by using the WMAP 5-year data as well as the distance 
information from BAO and SN, and Hubble's constant 
measured by HST. 

6.2.2. Analysis 

It is common to write the energy density of neutrinos 
(including anti- neutrinos), when they were still relativis- 
tic, in terms of the effective number of neutrino species, 
N eS , as 

= Nj-^Tt, (82) 



120 



where T v is the temperature of neutrinos. How do we 
measure Neg from CMB? 

The way that we use CMB to determine N e g is rel- 
atively simple. The relativistic particles that stream 



freely influence CMB in two ways: (1) their energy 
density changing the matter-radiation equality epoch, 
and (2) their anisotropic stress acting as an additional 
source for the gravitational potential via Einstein's equa- 
tions. Incidentally, the relativistic particles that do not 
stream freely, but interact with matter frequently, do 
not have a significant anisotropic stress because they 
isotropizc themselves via interactions with matter; thus, 
anisotropic stress of photons before the decoupling epoch 
was very small. Neutrinos, on the other hand, decoupled 
from matter much earlier (~ 2 MeV), and thus their 
anisotropic stress was significant at the decoupling epoch. 

The amount of the early ISW effect changes as the 
equality redshift changes. The earlier the equality epoch 
is, the more the ISW effect CMB photons receive. This 
effect can be measured via the height of the third acoustic 
peak relative to the first peak. Therefore, the equality 
redshift, z cq , is one of the fundamental observables that 
one can extract from the CMB power spectrum. 

One usually uses z cq to determine Q m h 2 from the CMB 
power spectrum, without noticing that it is actually z eq 
that they are measuring. However, the conversion from 
z e(l to £l m h 2 is automatic only when one knows the ra- 
diation content of the universe exactly - in other words, 
when one knows N e g exactly: 

1 , _ ^™ _ Qmh 2 1 , . 

+ cq n r ~ n 7 h 2 i + o.227iiv eff ' [ ' 

where fl^h 2 — 2.469 x 10~ 5 is the present-day pho- 
ton energy density parameter for T C mb = 2.725 K. 
Here, we have used the standard relation between 
the photon temperature and neutrino temperature, 
T v = (4/ll) 1 / 3 T 7 , derived from the entropy conserva- 
tion across the electron-positron ann ihilation (see, e.g., 
lWeinberelH97l iKolb fc Turner! [HfflOh . 

However, if we do not know N e g precisely, it is not 
possible to use z eq to measure Q, m h 2 . In fact, we lose our 
ability to measure fl m h 2 from CMB almost completely 
if we do not know N e g. Likewise, if we do not know 
£l m h 2 precisely, it is not possible to use z cq to measure 
N c g. As a result, N e g and O m /i 2 are linearly correlated 
(degenerate), with the width of the degeneracy line given 
by the uncertainty in our determination of z cq . 

The distance information from BAO and SN provides 
us with an independent constraint on f2 m /i 2 , which helps 
to reduce the degeneracy between z cq and fl m h 2 . 

The anisotropic stress of neutrinos also leaves 
distinct signatures in the CMB p ower spectrum , 
which is not degenerate with VL m h 2 dHu et al.l 11991 
iBashinskv k Seliald 12001. jTrotta k Melchiorrfpioll ) 
(see also iMelchiorri k Serrall2006t) have reported on evi- 
dence for the neutrino anisotropic stress at slightly more 
than 95% CL. They have parametrized the aniso tropic 
stress by the viscosity parameter, c 2 is (Hu 1998), and 
found c 2 is > 0.12 (95% CL). However, they had to com- 
bine the WMAP 1-year data with the SDSS data to see 
the evidence for non-zero c 2 , a . 

In iDunklev et al] (|2008f ) we report on the lower limit 
to N e g solely from the WMAP 5-year data. In this paper 
we shall combine the WMAP data with the distance in- 
formation from BAO and SN as well as Hubble's constant 
from HST to find the best-fitting value of N e g. 

6.2.3. Results 
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Fig. 18. — Constraint on the effective number of neutrino species, JV e fi (§ 16,2,31 1. Note that we have imposed a prior on N B g, < 7V c ff < 10. 
fn all panels we show the WMAP-only results in blue and WMAP+BAO+SN in red. (Left) Joint two-dimensional marginalized distribution 
(68% and 95% CL), showing a strong degeneracy between f2 m /i 2 and N B g . This degeneracy line is given by the equality redshift, 1 + z oq = 
n m /Q r = (4.050 X W 4 )n m h 2 /(1 + 0.2271A r cff ). The thick solid lines show the 68% and 95% limits calculated from the WMAP-only limit 



on z eq : 



3141 



-154 
-157 



(68% CL). The 95% CL contours do not follow the lines below N e g ~ 1.5 but close there, which shows a strong 



evidence for the cosmic neutrino background from its effects on the CMB power spectrum via the neutrino anisotropic stress. The BAO and 
SN provide an independent constraint on fl m h 2 , which helps reduce the degeneracy between N e ff and Cl m h 2 . (Middle) When we transform 
the horizontal axis of the left panel to z cq , we observe no degeneracy. The vertical solid lines show the one-dimensional marginalized 68% 
and 95% distribution calculated from the WMAP-only limit on z cq : z cq = 3141^^* (68% CL). Therefore, the left panel is simply a rotation 
of this panel using a relation between z cq , Q m h 2 , and N c ff. (Right) One-dimensional marginalized distribution of N e ff from WMAP-only 



and WMAP+BAO+SN+HST. Note that a gradual decline of the likelihood toward AT ei 



for the WMAP-only constraint should not be 



trusted, as it is affected by the hard prior, 7V e g < 10. The WMAP+BAO+SN+HST constraint is robust. This figure shows that the lower 
limit on JV cff is coming solely from the WMAP data. The 68% interval from WMAP+BAO+SN+HST, N cft = 4.4+1.5, is consistent with 
the standard value, 3.04, which is shown by the vertical line. 



Figure [T5] shows our constraint on N e g . The contours 
in the left panel lie on the expected linear correlation 
between Q m h 2 and N e g given by 



iV eff = 3.04 + 7.44 



n m h 2 3139 



1 



(84) 



. 0.1308 1 + z Cl 

which follows from equation (|8"3"|) . (Here, Sl m h 2 = 0.1308 
and z e q = 3138 are the maximum likelihood values from 
the simplest ACDM model.) The width of the degener- 
acy line is given by the accuracy of our determination of 
z oq , which is given by z cq = 3141+^ (WMAP-only) for 
this model. Note that the mean value of z eq for the sim- 
plest ACDM model with N cS = 3.04 is z cq = 3176±]jJ, 
which is close. This confirms that z cq is one of the fun- 
damental observables, and N e g is merely a secondary pa- 
rameter that can be derived from z cq . The middle panel 
of Fig. 1181 shows this clearly: z cq is determined indepen- 
dently of Neg. For each value of N e g along a constant z eq 
line, there is a corresponding Q, m h 2 that gives the same 
value of z eq along the line. 

However, the contours do not extend all the way down 
to N e g = 0, although equation (f8"4")l predicts that N c g 
should go to zero when Vt m h 2 is sufficiently small. This 
indicates that we are seeing the effect of the neutrino 
anisotropic stress at a hi gh significance. Whil e we n eed 
to repeat the analysis of iTrotta fc M clchiorril (|2005l ) in 
order to prove that our finding of N e g > comes from 
the neutrino anisotropic stress, we believe that there is a 
strong evidence that we see non-zero N e g via the effect 
of neutrino anisotropic stress, rather than via z cq . 

While the WMAP data al one can give a lower limit 
on N c g (jDunklev et al.ll200ll K they cannot give an up- 
per limit owing to the strong degeneracy with fl m h 2 . 
Therefore, we use the BAO, SN, and HST data to break 
the degeneracy. We find N c g = 4.4 ± 1.5 (68%) from 



WMAP+BAO+SN+HST, which is fully consistent with 
the standard value, 3.04 (see the right panel of Fig. [T5)) . 

7. CONCLUSION 

With 5 years of integration, the WMAP temperature 
and polarization data have improved significantly. An 
improved determination of the third acoustic peak has 
enabled us to reduce the uncertainty in the amplitude 
of matter fluctuation, parametrized by ag, by a factor 
of 1.4 from the WMAP 3- year result. The E-mode po- 
larization is now detected at 5 standard devi ations (c.f., 
3.0 st andard deviations for the 3-year data; iPage et ail 
2007), which rules out an instantaneous reionization at 
Zreion = 6 at the 3.5er level. Overall, the WMAP 5- 
year data conti nue to support the s implest, 6-parameter 
ACDM model dDunklev et al.ll2008h . 

In this paper, we have explored our ability to limit 
deviations from the simplest picture, namely, non- 
Gaussianity, non-adiabatic fluctuations, non-zero grav- 
itational waves, non-power-law spectrum, non-zero cur- 
vature, dynamical dark energy, parity-violating interac- 
tions, non-zero neutrino mass, and non-standard number 
of neutrino species. Detection of any of these items will 
lead us immediately to the new era in cosmology, and a 
better understanding of the physics of our universe. 

From these studies, we conclude that we have not de- 
tected any convincing deviations from the simplest 6- 
parameter ACDM model at the level greater than 99% 
CL. By combining WMAP data with the distance infor- 
mation from BAO and SN, we have improved the ac- 
curacy of the derived cosmological parameters. As the 
distance information provides strong constraints on the 
matter density (both BAO and SN) and Hubble's con- 
stant (BAO), the uncertainties in Q. m h 2 and Hq have 
been reduced by a factor of 1.7 and 2, respectively, from 
the WMAP-only limits. The better determination of Hq 



40 



Komatsu et al. 



0.90 - 



0.85 - 



g> 0.80 h 



0.75 - 



0.70 - 



- 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




WMAP 

WMAP+BAO+SN 



0.80 - 



0.75 - 



< 



0.70 - 



0.65 - 



0.12 0.13 0.14 



0.15 




Fig. 19. — Four representative cosmological parameters that have improved significantly by adding the BAO and SN data. (See also 
TableQ]) The contours show the 68% and 95% CL. The WMAP-only constraint is shown in blue, while WMAP+BAO+SN in red. (Left) 
The distance information from BAO and SN provides a better determination of Cl m h 2 , which results in 40% better determination of ag. 
(Right) The BAO data, being the absolute distance indicator, provides a better determination of Hg, which results in a factor of 2 better 
determination of f2yy, H;,, and f2 c . 



reduces the uncertainty in Q\ (as well as f2j, and tt c ) 
by a factor of 2, and the better determination of fl m h 2 
reduces the uncertainty in as by a factor of 1.4. These 
results are presented visually in Fig.HH See also Table Q] 
for the summary of the cosmological parameters of the 
ACDM model. The addition of BAO and SN does not 
improve determinations of fl^h 2 or r as expected. Since 
n s is mainly degenerate with fl^h 2 and r, with the for- 
mer being more degenerate, the addition of BAO and 
SN does not improve our determination of n s , when we 
consider the simplest 6-parameter ACDM model. 

To find the limits on various deviations from the sim- 
plest model, we have explored the parameter space by 
combining the WMAP 5-year data with the distance 
measurements from the Baryon Acoustic Oscillations 
(BAO) and Type la supernova (SN) observations. Here, 
we summarize significant findings from our analysis (see 
also Tabled]): 

• Gravitational waves and primordial power spec- 
trum: Improved simultaneous constraint on the 
amplitude of primordial gravitational waves and 
the shape of the primordial power spectrum (from 
WMAP+BAO+SN). In terms of the tensor-to- 
scalar ratio, r, we have found r < 0.22 (95% CL), 
which is the tightest bound to date. A blue pri- 
mordial spectrum index, n s > 1, now begins to be 
disfavoured even in the presence of gravitational 
waves. We find no evidence for the running in- 
dex, dn s /d\nk. The parameter space allowed for 
inflation models has shrunk significantly since the 
3-year data release (§ I3.3[) . most notably for models 
that produce significant gravitational waves, such 
as chaotic or power-law inflation models, and mod- 
els that produce n s > 1, such as hybrid inflation 
models. 

• Dark energy and curvature: Improved simultane- 



ous constraint on the dark energy equation of state, 
w, and the spatial curvature of the universe, fife 
(from WMAP+BAO+SN). Wc find -0.0179 < 
n k < 0.0081 (95% CL) and -0.14 < 1 + w < 
0.12 (95% CL). The curvature radius of the uni- 
verse should be greater than R curv > 22 and 
33 /i _1 Gpc for positive and negative curvature, re- 
spectively. The combination of WMAP, BAO, and 
SN is particularly powerful for constraining w and 
fife simultaneously. 

• Time- dependent equation of state: Using the 
WMAP distance priors (I a, R, z*) combined with 
the BAO and SN distance data, we have con- 
strained time-dependent w. The present-day value 
of to, wo, is constrained as —0.33 < 1 + w < 0.21 
(95% CL), for a flat universe. 

• Non-Gaussianity: Improved constraints on primor- 
dial non-Gaussianity parameters, —9 < /]^£ al < 
111 and -151 < f^ n < 253 (95% CL), from the 
WMAP temperature data. The Gaussianity tests 
show that the primordial fluctuations are Gaussian 
to the 0.1% level, which provides the strongest evi- 
dence for the quantum origin of the primordial fluc- 
tuations. 

• Non-adiabaticity: Improved constraints on non- 
adiabatic fluctuations. The photon and matter 
fluctuations are found to obey the adiabatic rela- 
tion to 8.9% and 2.1% for axion- and curvaton-type 
non-adiabatic fluctuations, respectively. 

• Parity violation: WMAP's limits on the TB and 
EB correlations indicate that parity- violating inter- 
actions that couple to photons could not have ro- 
tated the polarization angle by more than —5.9° < 
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Aa < 2.4° between the decoupling epoch and 
present epoch. 

• Neutrino mass: With the WMAP data com- 
bined with the distance information from BAO 
and SN, we find a limit on the neutrino mass, 
J2 m f < 0.67 eV (95% CL), which is better than 
the WMAP-only limit by a factor of 2, owing to 
an additional constraint on Hq provided by BAO. 
The limit does not get worse very much even when 
w is allowed to vary, as the SN data reduce cor- 
relation between Hq and w effectively. Since we 
rely only on the CMB data and distance informa- 
tion, our limit is not sensitive to our understanding 
of non-linear astrophysical effects in the large-scale 
structure data. 

• Number of neutrino species: With the WMAP 
data alone we find evidence for non-zero N c h 
(jDunklev et al.l l2008h . which is likely coming 
from our measurement of the effect of neutrino 
anisotropic stress on the CMB power spectrum. 
With the BAO, SN, and HST added, we break 
the degeneracy between tt m h 2 and N e g, and find 
N cS = 4.4 ± 1.5 (68% CL) and 1.8 < N cfi < 
7.6 (95% CL), which are consistent with the stan- 
dard value, N c g = 3.04; thus, we do not find any 
evidence for the extra relativistic species in the uni- 
verse. 

The limits that we have obtained from our analysis 
in this paper are already quite stringent; however, we 
emphasize that they should still be taken as a proto- 
type of what we can achieve in the future, including more 
integration of the WMAP observations. 

A smaller noise level in the temperature data will re- 
duce the uncertainty in non-Gaussianity parameters. An 
improved determination of the TE spectrum increases 
our sensitivity to non-adiabatic fluctuations as well as to 
the primordial gravitational waves. The E-mode polar- 
ization will be more dominated by the signal, to the point 
where we begin to constrain the detailed history of the 
rcionization of the universe beyond a simple parametriza- 
tion. Our limit on the B-mode polarization continues to 
improve, which will provide us with a better understand- 
ing of the polarized foreground. The improved TB and 
EB correlations will provide better limits on the cosmo- 
logical birefringence. 

While we have chosen to combine the WMAP data 
mainly with the distance indicators, one should be able to 



put even more stringent limits on important parameters 
such as r, n s , dn s /dhi k, w{z), m„, and N e g, by including 
the other data sets that are sensitive to the amplitude of 
density fluctuations, such as the amplitude of the galaxy 
power spectrum, Lya forest, weak lensing, and cluster 
abund ance. With the Lya forest data from lSeliak et al] 
(2006), for example, the limit on the running index im- 
proves from -0.068 < dn s /dlnk < 0.012 (95% CL) to 
-0.034 < dn s /d\nk < 0.011 (95% CL) for r = 0, 
and -0.11484 < dn s /d\nk < -0.00079 (95% CL) to 
-0.0411 < dn s /d\nk < 0.0067 (95% CL) for r ^ 0. 
Better understanding of systematic errors in these meth- 
ods will be crucial in improving our understanding of 
the physics of our universe. Hints for new physics may 
well be hidden in the deviations from the simplest 6- 
parameter ACDM model. 
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APPENDIX 
FAST CUBIC ESTIMATORS 
We use the following estimators for /]^£ al , //vi"\ and 6 src (the amplitude of the point-source bispectrum) 



/ 



local 
NL 



f7f = (F-'US, + (F- 1 ) 22 5 2 + (OasSs, 
bsrc = {F-^stSr + (F- v h 2 S 2 + (F-^asSa, 



where F™ is the Fisher matrix given by 
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Here, -B^) 2 ; are theoretically calculated angular bispectra (given below), where i = 1 is used for /]^£ al , i = 2 for , 

and i = 3 for 6 src . The denominator of contains the total power spectrum including the CMB signal (Cf mb ) and 

noise (N t ), Ci = Cf mb bf + N h and k is the beam transfer function given in lHill et all (|2008l ). 

While this formula allows one to estimate /]v£ al j /atl' 1 ' anc ^ ^src simultaneously, we find that a simultaneous estimation 
does not change the results significantly. Therefore, we use 

f^f = S 1 /F 11 , f N q f = S 2 /F 22l b src = S 3 /F 33 , (A5) 

to simplify the analysis, as well as to make the comparison with the previous work easier. However, we do take into 
account a potential leakage of the point sources into /]^£ al and /^jf by using the Monte Carlo simulation of point 

sources, as described later. These Monte-Carlo estimates of the bias i n /^ al an d /j^ 1 "' 1 due to the source contamination 
agree roughly with contributions from the off-diagonal terms in Eq. (|A1|) and (j A2|) . 

Assuming white noise, which is a good approximation at high multipoles in which noise becomes important, one can 
compute the noise power spectrum analytically as 

f d 2 h a%M(n) , k x 

Ni = n pix / , 7 ° )J , A6 
J 47r/ sky A obs (n) 

where f2 p i X = 4ir/N p i x is the solid angle per pixel, M(n) is the KQ75 mask, / SKy = 0.718 is the fraction of sky retained 
by the KQ75 mask, ao is the rms noise per observation, and 7V b s (n) is the number of observations per pixel. 
The angular bispectra are given by 

lr[a h (r)(3 l2 (r)/3 l3 (r) + (eye.)}, (A7) 

r 2 dr { \p h (r) lh (r)S h (r) + (eye.)] - 2S h (r)S h (r)S h (r)} , (A8) 

(A9) 
where 
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Various functions in j are given by 





[ k 2 dkg T i(k)ji(kr) 7 


(All) 




f k 2 dkP 9 (k)g T i(k)ji(kr), 


(A12) 




f k 2 dkP^ 3 (k)g T i(k)ji(kr), 


(A13) 


w4j 


f k 2 dkP* /3 (k)g T i(k)ji(kr). 


(A14) 



Here, P&(k) cx fc" s_4 is the primordial power spectrum of Bardeen's curvature perturbations, and gTi{k) is the radiation 

transfer function that gives the angular power spectrum as C\ = (2 / tt) [ k 2 dkP$(k)g%,(k). 

Th e skewness parameters, Si, are given by (jKomatsu et alj|2005t iBabichl 120051: ICreminelli et al.ll2006t TYadav et all 

51 = 4vr J r 2 dr J ^ [A(h, r)B 2 (n, r) 

-2B(n,r)(4 im (n,r)B sim (n,r)) M c - A(n,r)(B s 2 im (n,r)) M c] , (A15) 

5 2 = -3S 1 + 24tt / r 2 dr [ — {[B(h, r)C(n, r)D(h, r) - B(h, r)(C sim (h, r)D sim (h)) MC 

J J w 3 

-C(n,r)( J B sim (n, r)D sim (n)) M c - D(h, r)(B sitn (n, r)C sim (n, r)) M c] 
~ [D 3 (n, r) - 3D(h, r)(D 2 sim (h, r)) MC ] } (A16) 



3 

S 3 = \ [ — [E\h) - ZE(n){E 2 sim (n)) MC \ , (A17) 
■j J w 3 

where w 3 is the sum of the weighting function cubed: 

w 3 ee J d 2 n W 3 {n). (A18) 
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For a uniform weighting, the weighting function is simply given by the KQ75 mask, i.e., W^(n) = M(n), which gives 

w 3 = 47rJ 
given by 



W3 = in j ' s ky. For our measurements of /]vx a \ /jvl''; an d b sr c, we shall use a "combination signal-plus- noise weight," 



W(n) 



M(n) 



+ a%/N ohs (h)' 



(A19) 



cmb 



where cr 2 mb = (1/Att) 5^(2Z + l)Cz cmb fr? i s the CMB signal variance, ao is the rms noise per observation, and A/obs(n) 
is the number of observations per pixel. This combination weighting yields a nearly optimal performance for b src and 
/ivz, ' whereas it results in a minor improvement in /J^£ al over the uniform weighting. The bracket, ()mCj denotes 
the average over Monte Carlo realizations, and "sim" denotes that these are the filtered maps of the Monte Carlo 
realizations. 

The filtered temperature maps, A, B, C, D, and E, are given by 



A(h,r 
S(n, r 
C(n,r 
D(n, r 



E 



1=2 m=-l 



^ ai(r)^-ai m Yi m (h), 



1=2 m=-l 
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1=2 m=-l 

/max I 



k 
Ci 



1=2 m=-l 1 



(A20) 
(A21) 
(A22) 
(A23) 
(A24) 



respectively. Here, /max is the maximum multipole that we use in the analysis. We vary i max to see how much the 
results depend on Z max . 

Eq. (|A15[) and (|A16p involve the integrals over the conformal distances, r. We evaluate these integrals as 



J r 2 dr[...}(r) = J2(m)o P t^An{...](n). 



(A25) 



We use the bispectrum optimization algorithm described in ISmith fc Zaldarriaga! (|2006f) to compute the optimal 
weights, (u>i )opt, and decide on which quadrature points, r,, to keep. We choose the number of quadrature points such 
that the bispectrum computed in this way agrees with that from more dense sampling in r to 10 _o , which typically 
gives ~ 5 quadrature points for /j^£ al and ~ 15 points for . 53 
The measurement of these estimators proceeds as follows: 



1. Generate the simulated realizations of CMB signal maps, T5, from the input signal power spectrum, C ; cmb , 
the beam transfer function, 6/. We have generated 300 realizations for the analysis given in this paper. 



and 



2. Add random noise, Tjv, using the rms noise per pixel given by <tq / \J A bs(n) 

3. Add point sources. We use a simplified treatment for the source simulation, 



T src (n) 
2.725 K 



sinh 2 (a;/2) fircAW 
x 4 67.55 MJy 



(A26) 



where ft p i X is the solid angle of pixel, 



hp/(k B T cmh ) = 56.80 GHz (for T cmb = 2.725 K), e is a Poisson 



random variable with the mean of (e) = n STC £l p i X , and n src is the average number of sources per steradians. This 
simplified model assumes that there is only one population of sources with a fixed flux, -Fs rc , and each source's 
flux is independent of frequency. We choose n src = 0.85 sr _1 and F src — 0.5 Jy, which yields the source power 
spectrum in Q band of C ps = 8.7 x 10~ 3 /1K 2 sr and the source bispectrum in Q band of b src — 8.7 x 10~ 5 /1K 3 sr 2 , 
which roughly reproduce the measured values. However, this model does not reproduce the source counts very 
well. (The source density of n src = 0.85 sr -1 at 0.5 Jy is too low.) The main purpose of this phenomenological 
model is to reproduce the power spectrum and bispectrum - we include point sources in the simulations, in order 
to take into account the potential effects of the unresolved sources on primordial non-Gaussianity, /j^£ al and 

/•equil 

Jnl ■ 



53 Note that ISmith fc Z aldarriaga! pOul l use 10 6 as a criterion, which gives more quadrature points to evaluate. We find that 10 
sufficient for the size of statistical and systematic errors in the current measurements. 
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4. Coadd them to create the simulated temperature maps, T(n) = Tg(fi) + Tjv(n) + T src (n). 

5. Mask and weight the temperature maps, T(n) — > T(n) = W(h)T(n), where W-^n) is given by Eq. (|A19[) . 

6. Remove the monopole and dipole from T(n). 

7. Compute the harmonic coefficients as 

a lm = [ d 2 n f (n)Y* n (n). (A27) 



Generate the filtered maps, A s i m , B sim , C S j m , D sim , and -E s i m , and compute the appropriate Monte-Carlo averages 
such as (Asim-B S j m )MC! e ^ c - This is the most time consuming part. 

Compute the filtered maps from the WMAP data. When we coadd the V and W band data, we weight them as 
T v+W = (T V +0.9T W )/1.9. The beam transfer function of the coadded map is given by 6^ +w = (6^+0.96^)/1.9, 
and <7o / ./Vobs of the coadded map is given by 



(a v + w ) 2 _ 1 \ K) 2 , (0.9^ 



<t W (n) 1-9 2 



v^2 /n o^W\2- 

(A28) 



10. Compute the skewness parameters, Si, from the filtered WMAP data, and obtain /]^£ al , /jvl'', and b src , either 
jointly or separately. 

11. Compute these parameters from the simulated realizations as well, and obtain the uncertainties. 

For the computations of gTi{k) and generation of Monte Carlo realizations, we have used the maximum-likelihood 
values of the WMAP 3-year data (the power-law ACDM model fit by the WMAP data alone with the Sunyaev- 
Zel'dovich effect marginalized) : n b = 0.0414, n cdm = 0.1946, n A = 0.7640, H = 73.2 km/s/Mpc, t = 0.091, and 
n s = 0.954 (|Spergel et alj |2007). These parameters yield the conformal distance to t = as cr = 14.61 Gpc. 

AXION 

In this Appendix we derive relations between the tensor-to-scalar ratio r, the axion mass density £l a h 2 ! , the entropy- 
to-curvature perturbation ratio a, the phase of the Pecci-Quinn field a , and the axion decay constant f a . 
Let us write the expectation value of the complex Pecci-Quinn field, V'pQj as 

(^P Q ) = A e w 0) (Bl) 

where f a is the axion decay constant, and 9 a is the phase. Quantum fluctuations during inflation generate fluctuations 
in the phase, 66 a , as 

M* = irr- ( B2 ) 

27T/ a 

As the number density of axions scales as the phase squared, n a cx the mass density fluctuation is given by 

Pa a nOaJa 

As the energy density of axions was negligible during inflation, the axion density perturbation, Sp a / p a , would produce 
the isocurvature perturbation. While radiation (including photons) is generated by decay of inflaton fields, (some 
of) dark matter is in the form of axions whose generation is independent of photons; thus, the entropy perturbation 
between photons and axions would be generated. We assume that axions were not in thermal equilibrium with photons 
in the subsequent evolution of the universe. 
The entropy perturbations and curvature perturbations are given, respectively, by 

k 3 Ps(k) fg Hi PP n (k) _Hi_ Hi 

where Q a < tt c is the axion mass density, H^ is the expansion rate during inflation at which the wavenumber k went 
outside of the horizon, 4>k is the value of inflation at the same time, e = —H/H as (M^/2)(V' /V) 2 is the usual slow-roll 
parameter, V{4>) is the inflaton potential, and M v \ = l/vo7rG is the reduced Planck mass. We have used the slow-roll 
approximation, j> m -V'/(3H) and H 2 ^/(3Af 2 j). 

Here, let us comment on our choice of m = 1, which makes k 3 P$ cx k 771 ^ 1 independent of k. Since k 3 Ps oc H% oc k~ 2e , 
where e = —H/H 2 , this choice corresponds to having a very small slow-roll parameter, e <g. 1. This is consistent with 
our limit on the curvature power spectrum, n s — 1 + 6e — 4t/ ~ 1 — 4?] < 1, where ry is another slow-roll parameter. As 
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7 J \Q a J VlO 12 GeVy l-a Q 
Alternatively, we can eliminate the axion decay constant, f ai to obtain 

4.7 x io- 12 fn r h 2 \ 12/r /n \ 2/7 



r = 



]l 0/7 \ 7 J \^aj 1-ao' 
This is equation ([48]) . 



the current limit is 1 — n s sa 4?y ~ 0.04, our approximation, m = 1, is valid for e < 0.01. It should be straightforward 
to extend our analysis to the case in which ra^l. 
By dividing Ps(fc) by Pn(k), we find the entropy-to-curvature perturbation ratio for axions, ao(k), as 

i - «o(fc) " i^(fc) niei{f a /M vl f 1 J 

At this point, it is clear that one cannot solve this constraint uni quely for any of e, f a , or 9 a . 
In order to break the degeneracy, we use the axion mass density (|Kawasaki fc Sek iguchi 20071. and references therein) 

n a h* = 1.0 x io- 3 7 ^ (toTo^vJ > ( B6 ) 

where 7 is a dilution factor, representing the amount by which the axion density could have been diluted by a late-time 
entropy production bet ween the QCD phase transition at ~ 200 MeV and the epoch of nucleosynthesis at ~ 1 Mpc. 

Combining equation (|B5j) and (jB6|) to eliminate the phase, 6 a , and using the relation between the tensor-to-scalar 
ratio r and the slow-roll parameter e, r — 16e, we find 

^(I,xl0-)^V^VW^V /6 ^. (B7) 



(B8) 



EQUATION OF STATE OF DARK ENERGY: A NEW PARAMETRIZED FORM 

In this Appendix, we describe the models of dark energy that we explore in § 15.41 Our goal is to obtain a sensible 
form of time-dependent dark energy equation o f state, w(a). One of the most commonly used form of w(a) is a linear 
form (jChevallier fe Polarskj|200ll : lLinderll2003f) 

w(a) =wo + (l — a)w a , (CI) 

where wq and w a parametrize the present-day value of w and the first derivative. However, this form cannot be 
adopted as it is when one uses the CMB data to constrain w(a). Since this form is basically the leading-order term 
of a Taylor series expansion, the value of w(a) can become unreasonably too large or too small when extrapolated to 
the decoupling epoch at z* ~ 1090 (or a* ~ 9.17 x 10~ 4 ), and thus one cannot extract meaningful constraints on the 
quantities such as wo and w a that are defined at the present epoch. 

To avoid this problem, yet to keep a close contact with the previous work in the literature, we shall consider 
an alternative parametrized form. Our idea is the following: we wish to keep the form given by Eq. (|C1|) at low 
redshifts, lower than some transition redshift, ztrans- However, we demand w(a) approach to —1 at higher redshifts, 
z > zt rans . This form of w(a) therefore has the following property: at early times, before the transition redshift, 
^trans, dark energy was just like a cosmological constant, and thus the dark energy density was nearly constant, i.e., 
Pdo(z > ztrans) ~ constant. Then, dark energy began to become dynamical at z ~ -Ztransi with the equation of state 
given by the conventional linear form, Eq. (IC1|) . 

Some of the properties of our form of w(a) are similar to those of "thawing models," (Caldwell fc Linderl l2005h 
in which a scalar field was moving very slowly initially, giving w(a) w —1 at early times, and then began to move 
faster towards low redshifts, causing w(a) to deviate more and more from —1 at low redshifts. Our parametrization 
can describe more general class of models than single scalar field models, as it allows for w to go be low —1. On the 
other hand, models that are based upon a single scalar field cannot have w < — 1 (e.g., lHull2005f) . The "Forever 
regular" parametrization explored in IWang fc Tegmarkl (|2004fl also appro aches to a constant density at early times, if 
the la te-time equation of state is w < —1. The "Kink model" explored in lBassett et al.l (|2002l . 120041 ): ICorasaniti et al.l 
(2004) also extrapolates a constant equation of state at early times to a different constant equation of state at late 
times. Our parametrization is more general than theirs, as their form allows only for a constant equation of state at 
late times. 

We wish to find a smooth interpolation between w ear i y = —1 and iui a t c = wq + (1 — a)w a . We begin by writing 

w(a) = w(a)f(a/a tlans ) + (-1) [1 - /(a/Wans)] , (C2) 

where atrans = 1/(1 + -^trans), and the function f{x) goes to zero for x <ti 1, and to unity for x ^> 1. Here, w(a) is the 
form of w at low redshifts. Any function that has this property is adequate for f(x). We choose: 

f(x) = \ [tanh(lnx) + 1] = (C3) 
2 x + 1 

which gives the desired form of the equation of state of dark energy, 
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Fig. CI. — Evolution of dark energy for wo = —1.1 and w' = 1, and various transition redshifts, ztrans = 0.5, 2, and 10, above which 
w(z) approaches to —1. {Top Left) Evolution of the dark energy to matter density ratio as a function of z. Note that the vertical axis has 
been multiplied by 10 6 . {Top Right) Evolution of the dark energy density relative to the dark energy density at present. The dark energy 
density was nearly constant at high redshifts above ztrans; thus, these models can describe the "thawing models," (Caldwell & Linder 2005) 
in which dark energy was nearly constant at early times, and had become dynamical at lower redshifts. {Bottom Left) Evolution of the 
equation of state, w{z) = Pd c { z ) I Pdo( 2 )- By construction of the model, w{z) approaches to —1 beyond ztrans- {Bottom Right) Evolution 
of the effective equation of state, w e ff{z), which determines the evolution of dark energy density as pde(z) = pdc(0)(l + z) 3 ^ 1+w ° il W. 

where 

w(a) = wq + (1 — a)w a . (C5) 

One nice property of this form is that it allows one to obtain a closed, analytical form of the effective equation of state, 
w e g(a), which gives the evolution of dark energy density, yOde(a) = p(0)a^ 1+w <= !! ^ : 

/■In a 

d In a' w(a') 



w cS (a) - 



1 

In a 



-1 



(1 - a)w a 1 + W + (1 + Otr 



In i 



ln< 



s)w a a 
In ■ 



atr 



1 + Otr 



(C6) 



This property allows one to compute the expansion rate, H(a) (Eq. (7j), and hence the distance (Eq. [2]), easily. 



Finally, we use the present-day value of w, wo 
parameters. They are related to wo and w a as 

l + w 



w(z — 0), and the first derivative, 



dw/dz 



z=o- 



1 + w 



1 



^trails 
W a 



1 + Otr 



Qtrans(l + Wq) 
(1 + Otrans) 2 



as free 

(C7) 
(C8) 



The inverse relations are 



1 + Wo = (1 + a tra ns)(l + W ), 

W a = (1 + atrans)^' + a t rans(l + W ) 

In the limit of very early transition, a tran s <C 1, one finds wo ~ wo and w' ' 
description of our form of w(a). 

Figure [Cl] shows the evolution of dark energy density, pde(z) = p(0)(l + z) 3 [ 1+tUoff ( a ^, the equation of state, w(z), and 
the effective equation of state, w e g(z), computed from Eq. (IC4|) and (|C6[) . We choose wo = —1.1 and w' — 1, which 



(C9) 
(CIO) 

w a , as expected. This completes the 
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TABLE Dl 

Comparison of ACDM parameters from WMAP+BAO+SN with various SN compilations 



Class Parameter Union" Union+Sys.Err.' Davis c Alternative 1 * 



Primary 


ioon 6 /i 2 


9 9K7+0-058 
— 0.059 


2.267 ± 0.059 


2.270 ± 0.060 


2.265 ±0.059 




Q c h 2 


0.1131 ± 0.0034 


0.1134±°;°°f 7 


0.1121 ± 0.0035 


0.1143 ±0.0034 




U A 


0.726 ± 0.015 


0.725 ±0.016 


n 7 on + 0.014 

u. /oz_ 015 


0.721 ± 0.015 




n s 


0.960 ± 0.013 


0.960 ±0.013 


0.962 ± 0.013 


0.960±° ;°^ 




T 


0.084 ±0.016 


0.085 ±0.016 


0.085 ±0.016 


0.084 ± 0.016 




A^(V) 


(2.445 ± 0.096) x 10~ 9 


(2.447±°;°g) x 10" 9 


(2.429™) x ID" 9 


(2.457™) x 10- 9 


Derived 




0.812 ± 0.026 


u.oio_ Q27 


0.807 ± 0.027 


0.817 ± 0.026 




Ho 


70.5 ± 1.3 km/s/Mpc 


70.4 ± 1.4 km/s/Mpc 


70.9 ± 1.3 km/s/Mpc 


70.1 ± 1.3 km/s/Mpc 




n b 


0.0456 ± 0.0015 


0.0458 ± 0.0016 


0451+ 0016 


0.0462 ± 0.0015 




n c 


0.228 ±0.013 


229+ 014 


0.223 ±0.013 


0.233 ±0.013 






u.iooo_ 0036 


1361+ 0038 


0.1348 ± 0.0038 


0.1369 ±0.0037 




■°rcion 


10.9 ± 1.4 


10.9 ± 1.4 


10.9 ± 1.4 


10.8 ± 1.4 




to* 


13.72 ± 0.12 Gyr 


13.72 ±0.12 Gyr 


13.71 ± 0.12 Gyr 


13.73 ±0.12 Gyr 



"Compilation by Kowalski et al. (2008) without the systematic errors included 
'Compilation by Kowalski ct al. (2008) with the systematic errors included 
c Compilation bv lDavis et al.l l|200^ i 

"Compilation us ed in the original versi on of this paper (version 1 of arXiv:0803.0547 ) 
Compilation by Kowalski et al. (2008) without the systematic errors included 
^Compilation by Kowalski ct al. (2008) with the systematic errors included 
Compilation by; Davi s et al., 1 20071 ) 
Compilation used in the original version of this paper (version 1 of arXiv:0803.0547 ) 
*fco = 0.002 Mpc" 1 . A^(fe) = k 3 P n (k) I '(2tt 2 ) (Eq. [B]) 

^"Redshift of reionization," if the universe was reionized instantaneously from the neutral state to the fully ionized state 

at ^rcion 

*The present-day age of the universe 

are close to the best-fitting values that we find in § 15.41 (see Fig. [M)) . We show three curves for the transition redshifts 
of Ztrans = 0.5, 2, and 10. We find that the form of w(z) that we have derived achieves our goal: w(z) approaches to 
— 1 and the dark energy density tends to a constant value at high redshifts, giving sensible results at the decoupling 
epoch. The dark energy density is totally sub-dominant compared to the matter density at high redshifts, which is 
also desirable. 

The constraints that we have obtained for wo and w' are not sensitive to the exact values of Zt ran s (see the right 
panel of Fig. [T4|) . This is because all of the curves shown in Fig. ICll are very similar at z < 1, where the BAO and SN 
data are currently available. 

COMPARISON OF SUPERNOVA COMPILATIONS AND EFFECTS OF SYSTEMATIC ERRORS 




w w Q w' 



Fig. Dl. — The same as Fig. 1151 but the systematic errors in the Type la supernova data are included. 

In Table IdTI we show the ACDM paramet ers derived from WM AP+BAO±SN, where we use various SN compilations: 
"Union" for the latest union compilation ([Kowalski et alj [20081 ). "U nion+Sys.Err . " for the union compilation with 
systematic errors included, "Davis" for the previous compilation bv iDavis et alJ (|2007t ). and "Alternative" for the 
compilation that we used in the original version of this paper (version 1 of arXiv:0803.0547). 

For the " Alternative" compilation we have combi ned measurements from the Hubble Space Telescope ()Riess et al.l 
12004120071 ). the SuperNova Legacy Survey (SNLS) (jAstier et al.ll2006l ), and the Equation of State: SupErNovae trace 
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Cosmic Expansion (ESSENCE) survey (|Wood-Vasev et al1l2007f ). as well as some nearby Type la supernovae. In the 
"Davis" and "Alternative" compilations, different light curve fitters were used for the supernova data taken by different 
groups, and thus these compilations were not as optimal as the union compilation, for which the same SALT fitter was 
used for all the supernovae samples. Moreover, the union compilation is the largest of all. For this reason we have 
decided to update all the cosmological parameters using the union compilation. 

Nevertheless, we find that all of these compilations yield similar results: the mean values shift no more than ~ 0.5a. 

The effects of the systematic errors in the Type la supernova data on the ACDM parameters are also very small 
for WM AP+BAO+SN; however, the effects on the dark energy parameters, wo and w' , turn out to be significant. In 
Fig. ID1I we show the two-dimensional joint constraint on wq and w' from WMAP+BAO+SN+BBN (see also § 15.4.2)) 
with the systematic errors included. Comparing this to Fig. 1151 where the systematic errors are ignored, we find that 
the constraints on wo and w' weaken significantly: we find wq = — 1.00±0.19 and w 1 = 0.11±0.70 with the systematic 
errors included, whereas wq = —1.04 ± 0.13 and w' = 0.24 ± 0.55 without the systematic errors. 
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